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Aerodynamic Coefficients of an Oscillating 
Airfoil with Control Surface in Two- 
Dimensional Subsonic Flow’ 


A. 1. vAN DE VOORENT 
National Aeronautical Research Institute, Holland 


SUMMARY 


The method of reference 1, yielding a direct solution in terms of 

) known functions for the aerodynamic forces acting on an oscillat- 

ing wing in two-dimensional subsonic flow, is extended to the case 

Pof an airfoil with a control surface hinging about its leading edge. 

Hinge moment coefficients are calculated for translation and rota- 

| tion of the whole airfoil, while, for rotation of the control surface 

lift, pitching moment and hinge moment coefficients are pre- 

sented. Results are given for Mach Numbers 8 = 0.35, 0.5, 

06,0.7, and 0.8; for chord ratios r = 0.1, 0.2, and 0.3; and fora 
'gumber of w values smaller than 1.0. 


List OF SYMBOLS 


rectangular coordinates measured from the mid-chord 

point of the airfoil (x = chordwise coordinate) 
= x/vV1 — 6? 
=y 

semichord 

ratio of control surface chord to total chord 

air speed 

speed of sound 

Mach Number, v/c 

frequency 

reduced frequency, vl/v 

w/(1 — B?) 

acceleration potential 

density of the undisturbed airflow 

elliptic coordinates (see reference 1) 

amplitude of translation (positive if downward ) 

amplitude of wing rotation about mid-chord (positive 
if trailing edge downward) 

amplitude of control surface rotation (positive in the 
same sense as B) 

K total force on airfoil (positive if downward ) 


Received September 9, 1952. 

* This investigation was sponsored by the Netherlands Aircraft 
Development Board. 

t Chief, Section for Flutter and Theoretical Aerodynamics. 


= total moment about mid-chord point (positive if tail 
heavy ) 

= control surface moment about hinge axis (positive if 
tail heavy) 


INTRODUCTION 


U HAS BEEN DESCRIBED in reference | how exact for- 
mulas for the aerodynamic force and moment acting 
on an oscillating wing in two-dimensional subsonic 
compressible flow can be derived. These expressions 
are obtained from the analytical method, originally de- 
veloped by Timman,’ which solves the boundary-value 
problem of the time-free wave equation with pertaining 
boundary conditions and conditions of regularity. 
Complete analytical solutions in terms of Mathieu and 
modified Mathieu functions, applying to a plain wing 
and to a wing with an aerodynamically balanced con- 
trol surface, have been elaborated by Timman in col- 
laboration with the present author.* The numerical 
computations based on this theory have been performed 
by the Computation Department of the Mathematical 
Centre at Amsterdam under the direction of van Wijn- 
gaarden and with the assistance first of Scheen and later 
of Berghuis. 

Since the coefficients for a plain wing have already 
been published, this paper contains only those coeffi- 
cients that are introduced by the addition of a control 
surface. All tabulated values should be accurate up 
to one unit in the last decimal supplied. 

To facilitate use in flutter calculations for tapered 
wings (strip theory), all coefficients are now being inter- 
polated to 2 in the range 0 < 2 < 1.0 with an interval of 
Q = 0.02. Tabulated results can be applied for at the 
National Aeronautical Research Institute. 
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00338 
00800 
01125 
01651 

02390 


0 
0 
0 
0 
0 
0 
0 
0 


Na’ 


00013 
00236 
00695 
01384 
02304 
03457 
06483 


00039 
00423 
01168 
02273 
03747 
05593 
09042 


00006 
00171 
00549 
O1121 
O1889 
02853 
04658 
06847 


.10173 


00005 
00214 
00863 
02226 
03314 
05004 
06837 


.09558 


00000 
00290 
00850 
02014 
02809 
04041 

05658 


—0.01838 


-0.01017 





Na” ny’ 
0 0.01781 
—0.00141 —0.01296 
0.00250 —0.01117 


0.00998 
0.00953 
0.00992 
0.01140 
0.01381 
0.01794 


0.00406 
0.00668 
0.00889 
0.01292 


0.02906 


0 -~0.02119 
0.00162 0.01437 
0.00296 0.01277 
0.00455 0.01247 
0.00769 0.01386 
-~0. 01567 





—0.01484 0.01921 
0.02231 0.02341 
Na” np’ 

0 0.03837 
0.00376 0.03005 
0.00653 —0.02541 
0.00911 0.02202 
0.01177 0.01888 
0.01466 —0.01564 
0.01790 0.01217 
0.02601 —0.00455 
0 0.04151 
0.00456 0.03042 
0. 00802 
0.01163 0.02222 
0.01588 0.01940 
0.02124 0.01699 
0.02830 0.01522 
0.04365 0.01459 
0 0.04494 
0.00364 0.03432 
0.00637 0.02954 
0.00907 0.02685 
0.01204 0.02503 
—0.01555 0.02381 
0.01986 0.02326 
0.02846 0.02395 
0.04074 0.02716 
0.06462 0.03564 
0 0.05034 
0.00399 0.03665 
0.00705 0.038194 
0.01154 0.02958 
0.01923 0.03036 
0.02578 0.03310 
0.03772 0.03982 
0.05383 0.04955 
0.08492 0.06588 
0 0.05992 
0.00454 0.04063 
0.00836 0.03679 
0.01296 0.03724 
-0.02221 0.04389 
0.02953 0.05090 
0.04323 0.06381 


0.06476 0.07950 


The graphs and tables in this paper will be 
panied by a short review of the complications intro- 


duced by the control surface. 


THE SINGULARITY AT THE 


As has been shown in reference 1, the problem is to 


LEADING 


CONTROL SURFACE 


EDGE 


accom- 


OF 


TABLE 2 ( Concluded ) 


np” k,’ k.” me’ m-” ne" Ne” 
0 1.53979 0 +0. 19937 0 0.03291 0 
0.00337 1 62627 + 0.38387 0.07622 0. 22769 0.03022 0.00148 
0.01040 0.81591 0.35330 0.20376 0.22813 0.02953 0.00504 
0.02004 -~0. 69501 0.29674 0.29421 0.20713 0.02970 0.00961 
0.03253 0.62506 0.36481 0.17390 0.03061 0.01481 
0.04019 0.59628 0.24081 0.39506 0. 14836 0.03134 0.01762 
0.05018 0.55894 0.23085 0.42063 0. 10686 0.03232 0.02098 
0.05935 0.51803 0.22050 0.42531 0.06122 0.03310 0.02401 
0.07069 0.45669 +0. 18922 0.39204 0.00554 0.03364 0.02844 
0 1.83272 0 +0. 23730 0 0.03917 0 
0.00436 1.02022 +0.51461 0.21281 0. 28293 0.03543 0.00195 
0.01322 0.75804 0.44108 0.38230 0.22049 0.03533 0.00629 
0.02222 0.64412 0.37695 0. 46206 0.13816 0.03629 0.01002 
0.03493 0.55952 0.32591 0.48864 0.01550 0.03819 0.01395 
0.04177 0.52173 0.30925 0.46383 +0. 04502 0.03907 0.01560 
0.05067 0.46675 0. 28732 0.38877 0.09445 0.03964 0.01778 
0.06007 0.39880 0.25137 0.28985 +0. 06636 0.03937 0.02150 
TABLE 3 
r=03 
nh be,” ke.” m-.’ me” n-" ne” 
0 1.41072 0 +0. 26935 0 0.05905 0 
0.00887 1.10235 +0. 20347 + 0.11341 0.15912 0.05355 0.00527 
0.02459 0.65413 0.14124 + 0.038775 0. 18334 0.05068 0.01515 
0.04092 0.88017 + 0.05735 0.00004 0.19634 0.04887 0.02539 
0.05740 0. 84006 0.02489 0.02039 0.21040 0.04740 0.03563 
0 07405 0.81785 0. 10268 0.03155 0.22717 0.04601 0.04586 
0.09093 0. 80634 6.17658 0.03733 0.24657 0.04458 0.05609 
0.12568 0.80340 0.31547 0.03961 0.29213 0.04145 0.07672 
0 1.52592 0 +0. 29135 0 0.06387 0 
0.01256 1.09364 +0. 25659 +0.06747 0.20377 -0.05655 0.00752 
0.03354 0.93220 0.17034 0.02279 0.22928 0.05383 0.02067 
0.05509 0.86719 + 0.07787 0.06627 0. 25033 0.05254 0.03391 
0.07715 0.84251 0.00360 0.09233 0.27646 0.05180 0.04707 
0.09988 0. 83928 0.07379 0.11149 0.30767 0.05137 -0.06018 
0. 12332 0.84958 0.13344 0.12857 0.34284 0.05120 0.07325 
0.15945 0. 88136 0. 20269 0.15601 0.40031 0.05142 0.09271 
0 1.65186 0 + 0.31539 0 0.06914 0 
0.00722 1.21649 +0.32838 +0.08740 0.22502 0.06206 0.00410 
0.02307 1.01721 0. 28630 0.02584 0.25472 -0. 05922 0.01406 
0.03959 0.92204 0. 22007 0.08836 0. 26894 0.05812 0.02426 
0.05639 0.87349 0.15956 0.12896 0.28395 0.05783 0.03437 
0.07355 0.84931 0.10845 0.16249 —0.30150 0.05803 -0.04432 
0.09110 0. 83897 0.06709 0.19162 0.32104 ~0.05860 -0.05409 
0.11810 0.83911 + 0.02074 —(0. 23409 0.35065 0.06009 —0.06842 
-0. 14516 0.84575 0.00628 —0. 27930 0.37705 0.06217 — 0.08213 
0.17920 0.85158 0.01722 0.34237 0.39907 0.06549 —0.09907 
0 1.85045 0 +0. 35332 0 0.07746 0 
0.00828 1.24729 +0. 42663 +0. 02943 0.28288 0.06831 0.00468 
0.02697 1.00777 0.37007 0.11794 0.30200 0.06583 0.01639 
0.05264 0.87695 0. 28028 0.22457 0.30679 0.06600 0.03169 
0.08576 0.81218 0.20454 0.31771 0.30843 0.06879 0.04954 
0.10588 0.79120 0.17705 0.36627 0.30253 0.07136 0.05937 
0.13181 0.76605 0.15520 0.42239 0.28111 0.07528 -0.07118 
0.15499 0.73444 0.14076 0.46148 0. 24396 0.07904 0.08145 
0.18221 0.67419 +0.10929 0.47221 0.17612 0. 08289 -~0.09508 
0 2.20248 0 +0. 42053 0 0.09219 0 
0.01057 1.25037 +0.58131 0.10 0.36417 ~0.07941 -0.00604 
0.03403 0.95069 0.48003 0.31453 0.32939 0.07865 -0.02038 
0.05770 0.82328 0.39118 0.42820 0.27136 0.08155 -0.03312 
0.09050 0.73610 0.31207 0.51120 0.16492 —0.08831 —0.04720 
0.10760 0.70204 0.28443 0.51795 0.09807 ~0 09227 —0.05299 
0.12870 0.65273 0.25122 0.47347 0.01874 0.09629 —0.05965 
0. 14894 0.59148 +0. 20789 0.36517 0.00342 0.09758 0.06863 


THE 


determine a function « related to the acceleration po- 




























































































tential y by the formula 

‘ > , . , / > 
u(X, Y) = o(XV/1 — B*, Yoexp (—ivBX/cV/ 1 — B?) 
which satisfies the time-free wave equation with cor- 
responding boundary conditions and conditions of 
regularity. Besides the usual singularity at the leading 
edge of the airfoil, there exists in the case of a control 


surface a second singularity, which is logarithmical for 
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Fic. 1. Aerodynamic coefficient, a. 


linearized flow, at the leading edge of the control sur- 
face. 


The corresponding singular part in « can be obtained 
from the function 


SCIENCES—MARCH, 


; cosh — — cos ( ) 
f(&, n) = tog | ms = = 


cosh € — cos (yn — m) 


1953 


representing the acceleration potential of the steady 
flow about an airfoil with zero angle of incidence and q 
control surface with an infinitely small angle of deflec- 
tion. In the same way as for the singularity at the lead- 
ing edge, correction terms must be added to obtain a 
solution of the wave equation. The result is 


cosh & — cos (y + m) 
cosh — — cos (n — m) 


fos Ne, (€) 
4> seat - ” 
1 N 


‘ Ne,“" (0) 
—mé 


- e sin mn, 
> ® B,,™ | Se, ” | 


m 1 m 


u(é, n) = v Jog | 


The factor by which this singular potential enters into 
the expression for the resulting potential is obtained by 
identifying the vertical velocity of the flow in an arbi- 
trary point of the control surface as calculated from 
the resulting potential with its prescribed value. Ob- 
viously, the factor is proportional to the discontinuity 
in the prescribed normal velocity at the leading edge of 
the control surface. From the complete function 
u(é, n) the acceleration potential and the pressure dis- 
tribution follow, and the aerodynamic forces and mo- 
ments can then be obtained by appropriate integrations. 


NUMERICAL CALCULATIONS 


The complex dimensionless derivatives k,, m., Na, M» 
and n, defined by 
K = mplv*e" (Ak, 4. Bk, + Ck,) 
M mpl v7e""(Am, + Bm, + Cm.) 
N = rpil*ve'"(An, + Bn, + Cn.) 


were calculated for the same 8 and w values for which in 
reference 1 the coefficients k,, ky, m,, and m, have been 
given, with the exception of w values larger than |, 
which were omitted in the calculations referring to a 
wing with control surface. Three values were given to 
the chord ratio; viz., 0.1, 0.2, and 0.3. 

As new basic functions the incomplete Bessel func- 
tions, 


R,(6?2, m) = Si” e265 9 Cos nn dy 


appeared in the calculations for the control surface 
system. The method of computation was by expand- 
ing e’°°S” into its Fourier series and performing 
the integration. This leads to a series of Bessel func- 
tions J,,(8°2), m = 1.2..., with coefficients depending 
onm. In this way the incomplete Bessel functions were 
computed‘ in eight decimal places for n = 0(1)20, 
cos 7 = 0.3(0.1)0.8, and for a varying number (47 for 
n = 0to1 for n = 20) of B?Q values between 0 and 4.5. 

The resulting values for the aerodynamic derivatives 
are given in Tables 1-3. Partially, they can be com- 
pared with results given by Dietze® [w = 0, (0.02), 0.1, 
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(0.1), 0.7, 8 = 0.7, r = 0.15] and recalculated and ex- 
tended by Turner and Rabinowitz,’ using Dietze’s 
method [w = 0, (0.02), 0.1, (0.1), 0.7, 8 = 0.7, 7 = 0.15, 
(0.09), 9. 12]. Though an accurate comparison is 
hampered by the different chord ratios, general agree- 
ment is good, the largest discrepancies being 4 per cent 
5 deg. in the argument, which 
It may, how- 


in the modulus and 
occurs in k, for large values of w and r. 
ever, be repeated that the coefficients given in this 
paper are the more accurate. (See Figs. 1-5.) 
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On the Pure Rolling Motion of Winged and/or 


Finned Missiles in Varying Supersonic Flight’ 


JOHN D. NICOLAIDES?# anp RAY E. BOLZi 
Ballistic Research Laboratories and Case Institute of Technology, Respectively 


ABSTRACT 


Based on the assumption that the aerodynamic force acting at 
any point on a lifting surface is linearly dependent upon the local 
angle of attack at that point and that the aerodynamic coefficients 
vary in a particular manner with Mach Number, the pure rolling 
motion of a winged and/or finned missile with large or small vari- 
ation in flight Mach Number is given. 

This refined theory of pure rolling motion and a suggested re- 
duction technique are investigated in the Aerodynamics Range 
for small changes in Mach Number where it is found that the 
“*fit’’ of the theory to the observed motion is consistent with the 
small errors of measurement. The results obtained from the re- 
fined theory of motion are compared with those obtained from the 
standard theory based on constant aerodynamic coefficients, and 
it is found that for these particular tests the differences are ex- 
tremely small and may be reconciled by a simple empirical rela- 
tion. 

The aerodynamic coefficients, Cip and Cis, and the effective- 
ness parameter, pb/2V, associated with the pure rolling motion 
of an arrow-type missile configuration having rectangular, single 
wedge, 8 per cent thickness fins are determined over a Mach Num- 
ber range of 1.3 < M < 3.0. These experimental values are 
compared with those obtained from previous experiments on the 
same configuration with 16 per cent thickness fins and with theo- 
retical values as determined by the linearized theory for rectangu- 
Qualitatively, the results show 
that the values for the 8 per cent thickness fins lie in between. 


lar fins of 0 per cent thickness. 


The aerodynamic coefficients are also obtained at two Mach 
Numbers on two additional configurations, wing only and wing 
plus fins, to indicate the suitability of the Range for a study of 
wing-fin roll interference. The results indicate that for the two 
models studied the damping effectiveness of the fins was almost 
eliminated, being 14 pér cent at a Mach Number of 1.95 and 0 
per cent at a Mach Number of 1.73. 


SYMBOLS 


A p = area of fin or wing panel 

A,B = arbitrary constants depending on boundary condi- 
tions 

b = wilg span 

i, = damping constant 

Cx = control constant 

Cp = drag coefficient, drag/aS 

Cis = roll moment derivative due to canted surface, 
(OC, /05)5—>0 

Clip = roll moment derivative due to rolling velocity, 
[OC /2( pb/2V) l¢pn/21)—>o 

Ci = roll moment coefficient, L/gAb 

d = diameter of body 
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6L5 = rolling moment due to angle of cant, ( CisdgbA pny 
, ply = rolling moment due to rolling velocity, damping 
moment, Cip(pb/2V )qbA pn 
I = axial moment of inertia, slugs per sq.ft 
M = Mach Number 
m = mass, slugs 
ns = number of canted fins 
n = total number of fins and/or wings 
p = Rolling velocity, ¢, rad. per sec. 
gq = pV?/2 
PE. = probable error 
pb/2V = fin tip helix angle 
psb/2V = steady-state fin tip helix angle 
S = xd?/4 
5 = steady-state rolling velocity 
V = linear velocity of missile along trajectory, ft. per sec 


Z = linear distance along the missile flight trajectory 


a, 8, X = constants 
¥ = —d/[1 + (B/aM)) 
6 = angle of cant 
p <« = free-stream air density, slugs per cu.ft 
o = angle of roll 
Wy, 2) = In(1 + yz)/y2 
Superscripts 


Prime mark is the derivative with respect to 2, distance along 
the trajectory. 


INTRODUCTION 


I REFERENCES | AND 2, a technique was suggested 
for determining the aerodynamic coefficients asso- 
ciated with the pure rolling motion of free-flight missiles. 
This technique was applied to free-flight data obtained 
in the Aerodynamics Range with the result that the 
aerodynamic coefficients were determined with good 
statistical accuracy and were reproducible. 

However, the theory of rolling motion employed was 
derived under the assumption that the aerodynamic 
coefficients were constant during any one test flight. 
For the Aerodynamics Range where the observational 
distance is approximately equal to the significant tran- 
sient rolling distance for the models employed, the total 
change in Mach Number was 7 per cent maximum, for 
which the associated variation of the aerodynamic co- 
efficients from the average value was about 3 per cent. 

Therefore, in view of the change of the coefficients 
associated with the variation in Mach Number both 
in the Aerodynamics Range and also in full-scale free 
flight, it is one of the purposes of this paper to present a 
refined theory of pure rolling motion of free-flight 
missiles which allows for variation of the aerodynamic 
coefficients with the varying flight Mach Number. 


The value of this theory is to: (1) provide a more 
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general reduction technique applicable to other test 
facilities where the change in the coefficients is large, 
and (2) investigate the association between Mach 
Number and the values of coefficients as obtained in the 
Range by the previous theory. 

Also in reference 1, the aerodynamic coefficients in 
roll are presented over a range of test Mach Numbers 
from 1.7 to 2.6 for missiles with 16 per cent thick fins. 
As discussed in the reference, the deviation of the ex- 
perimental results from those calculated from the linear- 
ized theory was thought to be primarily a result of the 
large wing thickness. It is, as a consequence, a second 
purpose of this paper to present data for identical test 
configurations except with 8 per cent thick fins, over 
arange of Mach Numbers from 1.3 to 3.0, in order to 
evaluate the effect of fin thickness on the aerodynamic 
coefficients. 

Finally, it is the third purpose of this paper to ex- 
plore in a preliminary fashion the applicability of the 
test technique to an experimental determination of 
wing-tail interference phenomena. 


SOLUTION OF THE ROLL EQUATION 


The differential equation of motion! of a missile in 
pure roll is given by 


o” + C(z)d’ = Ci(z) (1) 


where 


b'nAp : aa ’ 
. * Ci,(z) — “2 C,(3) (la) 


Ci(z) 
2m 


MOTION OF 
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pbn;A pb 


Cis(2) (lb) 


rn) = roll angle 


and where the primes refer to derivatives with respect 
to s, the distance along the trajectory. In reference 
1, C; and C, were taken as constant; 
they are functions of the Mach Number, ./, because of 


here, however, 


the variation of the aerodynamic coefficients with /, 
itself a function only of z. 

From Eq. (1), one may observe that, for a constant 
rolling velocity, ¢” = 0, the ratio C,/C, is a constant 
and 1s the steady-state rolling velocity of the missile. 
As such, it is directly proportional to the steady-state 
fin tip helix angle, pb/2V [reference 1, Eq. (346)]. For 
many configurations the value of pb/2V has-a small 
variation with Mach Number (e.g., see Fig. 14 of refer- 
ence 1) and thus C,(z)/C,(z) might be treated as a con- 
stant; however, for generality we shall let 


C.(z)/Ci(z) = 5 + s’2 (2) 


where s; is the steady-state rolling velocity correspond- 
ing to the Mach Number at the beginning of the tra- 
jectory. 

The solution of Eq. (1) may thus be written in in- 


, . 
So, {Ci(s) dz 

5 2° + Ay € * —_— 

vf alii “J ef 492 dz (3) 


The variation of C\(z) from Eq. (la), neglecting the 
drag term* is 


tegral form as 


ri) = By + 512 + 


Ci(s) = AiC,,(2) (4) 
where 

Ky = — pb*nA P 4] (5) 
Now C,,(z) may be expressed with good accuracy both 
theoretically and experimentally by a linear relation 


C.* = aM+86 (6) 


(a, 8 = constants) for the configuration of reference 
1‘ and many others investigated (e.g., see Figs. 1 and 
2).4 

This variation of C,,, together with a linear vari- 
ation of M(z), approximating the true variation over 
relatively small Mach Number ranges, or a variation 
of C, ~ 1/M for large Mach Number ranges, yields 


M = M\(1 — dz) (7) 


* The contribution of the drag term to C; may be estimated by 
considering the values in Table 4, reference 1, where it is seen that 
\(Kr/m)/C,|100 = 1.73 per cent maximum (for Model 6) 

t The variation being 2-3 per cent for 1.5 < M < 3.0 for the 
configuration of reference 1. 

t Data from references 6, 7, and 8 were used in preparing Figs 


l and 2. 
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Fic. 2. Variation of C;,~! with Mach Number for three delta 
plan forms 

where 1/4) = Mach Number at the beginning of the 

trajectory and ’ = constant, when substituted into 


Eq. (5) yields* 
Ci(z) = C,,/(1 + yz) (S) 
where 


Cy, = Kk, (8 + aM) \ 


| 9) 
y = —d/[1 + (6/aM)]f ' 


Now, employing this variation of C\(z), the refined roll 
equation becomes 


s’ 5” ¥ 
atid ola i‘ (1-< Jefe + 
C.ty 2 ios 
Ae (Cig — ¥) /y] In (1 * ys)} (10) 
When s’ = y = 0, the equation reduces to 
oH) = B + S12 + Ae C192 (1 ; ' 


the relation of the previous paper. ' 
The aerodynamic coefficients are obtained from the 
constants of the roll equation [Eq. (10) |] by 


Cis ws (27 pbA pn55)51Ci, (12) 
4] no Oe 

C= = (c., pon ) (13) 
pb?A pn 2m 


bb 51) | l | 
-/6 = a : (14) 
2] 2(57.3)8L1 + (pSC,/2mC,,) 


* The authors are indebted to R. A. Turetsky, Ballistician, 


Naval Ordnance Laboratory, for his assistance in the development 
of the theory and, in particular, for his contributions to the re- 


duction technique 
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METHOD OF DATA REDUCTION? 


The method of data reduction is fundamentally the 
same as in references 1 and 2. The roll equation [Eq. 
(10)] is expanded in a Taylor’s series about the initia] 
values of the constants, and the higher order terms are 
neglected. The resultant equation is linear in the dif. 
ferentials of the constants, and thus the method of 
Least Squares may be applied directly to determine 
the values of these differential corrections. The cor. 
rected initial values of the constants, obtained by add- 
ing these differential corrections to the original values of 
the constants, may in turn be used as initial values, and 
the process may be repeated until the sum of the 
squares of the residuals of the fit reaches a minimum, 


EXPERIMENTAL TESTS 


The purpose of the experimental program may be 
stated in four parts: 

(1) To determine the accuracy and consistency of the 
refined form of the theoretical roll equation in repre- 
senting the experimental free-flight data obtained in 
the Aerodynamics Range (for description see references 
1-4). 

(2) To compare the results obtained by applying the 
refined theory and the standard theory to data ob 
tained in the Aerodynamics Range for small variation 
in flight Mach Number. 

(3) To determine the effect of fin thickness on the 
aerodynamic coefficients over a large range of Mach 
Numbers and to compare these values with those ob- 
tained from linearized theory. 
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Fic. 4. Measured values of the missiles’ roll angles at the 


Range station locations plotted with the ‘“‘fitted’”’ theoretical 
curve of rolling motion. 


(4) To investigate the applicability of the Range for 
experimental studies of wing-fin interference phenomena 
for rolling missiles. 

The basic design of the missile used in the program 
is identical to that used in the tests of reference | 
namely, a cone-cylinder fuselage (1/d = 10) with fixed 
rectangular-plan-form single-wedge fins and/or wings 
incruciform arrangement. The three variations in the 
basic design are: (1) fin-only models 8 per cent thick 
at the trailing edge (Fig. 3a); (2) wing-only models 16 
per cent thick at trailing edge (Fig. 3b); and (3) wing- 
fin configuration (Fig. 3c) with 16 per cent thick wings 


and fins. 
RESULTS AND DISCUSSION 


Presentation of Results 


To illustrate the data obtained from the Range, Fig. 
4 presents the angular orientation in roll of three repre- 
sentative tests as measured from the shadowgraphs 
(Photos 1, 2, and 3) taken at each spark station along 
the trajectory. 

In this same figure are plotted the equations of the 
curves that best fit the data as obtained by the meth- 
ods discussed in the previous section. The agreement 
of the data to within a degree of the fitted curve cannot 
be shown in the figure because of the insensitivity of 
the scale. (The difference between the curves should 
not be interpreted as a result of the variation in the con- 
figuration designs only, since the moments of inertia 


are not identical. ) 
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The ‘excellent fit’ of both roll equations to the ex- 
perimental points for the various models is given by 
their probable errors (P.E.), which varied from 0.14° 
to 0.64° and had an average of 0.41°. Since the prob- 
able errors are small and of the same order of magni- 
tude as the estimated error in roll angle measurements, 
it may be concluded that both roll equations represent 





Typical shadowgraphs of three types of 


Puortos 1, 2, AND 3. 
missiles in free flight showing roll measurement pins and super- 
sonic flow phenomena. 





the actual motion of the models tested to within the 
experimental error. The small P.E.’s also indicate 
that rolling motion is unaffected by the yawing motion 
for the mean square yaws encountered in the tests. 


Discussion of Results 


(a) Comparison of the Two Reduction Techniques. 
For five projectiles, the observed motion was ‘“‘fitted”’ 
and the aerodynamic coefficients were determined both 
by means of the Standard Reduction Technique of refer- 
ence 1—assuming the coefficients constant over the 
Mach Number variation for any given flight—and by 
means of the Refined Reduction Technique discussed 
in this paper, which allows for variation of the coeffi- 
cients. From the values obtained from the two tech- 
niques, it was seen (1) that the total change in Mach 
Number for each of the models is small, (2) that both 
theories “‘fit’’ the experimental data equally well and 
considerably within the previously estimated meas- 
urement accuracy of 0.7°, (3) that the constants and 
the aerodynamic coefficients are also equally well de- 
termined by both theories as indicated by their prob- 
able errors, and (4) that a small (from 0.5—2 per cent) 
but consistent trend toward increasing the values of the 
coefficients results from the refined technique when the 
center of the timing data is used as a reference point 
(see Figs. 5 and 6). 

In reference to this last observation it should be noted 
that the values of the aerodynamic coefficients calcu- 
lated from the constants of the Refined Technique are 
those at 2, the center of the timing data which, for 
this discussion, may be considered as the center of the 
observational data; whereas the values of the aero- 
dynamic coefficients obtained from the Standard Tech- 
nique were arbitrarily associated with 2, as in refer- 
ence 1. Since the Refined Technique yields the vari- 
ation of the aerodynamic coefficients with distance along 
the trajectory, it is now possible to determine the dis- 
tance down range at which the values of the coefficients 
from the Standard Technique occurred. It was found 
that this distance could be approximated by the follow- 
ing empirical relation: z at ¢’ = (1/2) (@,’ + dp’) 
where ¢,’ = rolling velocity at last station and ¢,’ = 
rolling velocity at first station, which effectively weights 
the stations at the beginning of the range where the 
transient motion is most pronounced. 

Comparing now the differences between the aerody- 
namic coefficients obtained from the two reduction 
techniques and the use of this empirical relation, it is 
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Fic. 5. Experimental values of C,, as a function 
Number plotted with the ‘‘Weighted Fit’’ curve 
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TABLE | 
(B) (C) (D) 
Standard Refined Per Cent 

(A) Reduction Reduction Difference in 
Round Values of Values at s* Values (C — B) of 

Number ch Cis ch Cs Ci, Cis 
6 0.530 0.819 0.541 0.834 1.94° 1. 88°; 
rg 0.511 0.749 0.519 0.760 1.49° 1.49% 
& 0.541 0.798 0.547 0.806 1.14% 1.10% 
Q 0.506 0.765 0.515 0.779 1. 88° 1.84°, 
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Fic. 6. Experimental values of Cig as a function 
Number plotted with the ‘‘Weighted Fit’’ curve 


(E) 
Refined Reduction 
Values at s When Per Cent 
¢’ = (1/2) X 
(or’ + ox’) Values (E 
C!, Cis cs 
0.539 O.817 0.13 
0.510 0.748 —(). 16° 
0.542 0.798 0. 10° 


0.507 0.767 0.23% 
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Fic. 7. Experimental values of C),~! as a function of Mach 


Number plotted with the ‘‘Weighted Fit’’ curve 
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Experimental values of Cis 
Number plotted with the ‘Weighted Fit”’ curve. 


Fic. 8 ' as a function of Mach 


sen from Table 1 that the agreement is considerably 
improved. Although the empirical relation appears to 
provide a good method for associating the aerodynamic 
coefficients obtained from the Standard Technique with 
Mach Number, it should be noted that the original dif- 
ferences were small (0.5-2 per cent) and of little con- 
sequence in most engineering work. 

Therefore, the comparison of the two reduction tech- 
liques as applied to this particular Aerodynamic Range 
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data reveals that the Standard Technique is certainly 
adequate. However, for models and full-scale missiles 
where the change in the aerodynamic coefficients for a 
particular test run are large, the Refined Theory should 
yield significantly better results. 

(b) Linearity of (Cis~') and (C,,)~! with Mach Num- 
ber; Consistency of the Data.—\t was found |see Eq. 
(6)] both experimentally and theoretically that, for 
all of the many configurations investigated, (Ci) ~' and 
(C,,)~!, when plotted against Mach Number, yielded 
a linear variation. Since the Refined Roll Reduction 
Technique was based on this variation, the values of 
the reciprocal of the coefficients obtained from the 
various firings are plotted in Figs. 7 and 8. These ex- 
perimental values* were “‘fitted’’ by straight lines in a 
manner that weighted the influence of each test value 
by the probable error associated with it. 

In Figs. 5 and 6, the coefficients C,, and Cy are 
plotted in conventional manner against Mach Number. 
Also included in the figures are the weighted lines ob- 
tained from the reciprocal plots. It is noted that the 
maximum deviation of any of the experimental values 
from the weighted line is about 5 per cent (for > 
1.4) in both figures. 

(c) Effectiveness Parameter per Cant Angle (pb, 2V/6). 

In order to obtain the consistency of the effective- 
ness parameter with Mach Number, a weighted straight 
line was ‘‘fitted’’ to the data (.\/ > 1.4) and is plotted 
in Fig. 9, together with the data. The maximum devi- 
ation of the data from the “‘fitted’’ line is 6 per cent 
(Model 16). It should be noted that the values of both 
Cis and pb/2V/6 for Models 15 and 16 differ from the 


* The coefficients obtained for 1J < 1.4 have been excluded 


because of the suspected presence of transonic effects 
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The roll moment coefficient produced by a canted fin 
surface plotted as a function of the cant angle. 
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respective curves by approximately 4.8 and 7.8 per cent, 
respectively. Since the determination of both of these 
coefficients depends on an accurate knowledge of the 
angle of cant and since both models were subject to 
relatively severe launching conditions in order to ob- 
tain the largest Mach Numbers of the program, it is 
conjectured that there may have been distortion of the 
aluminum fins. 

One of the assumptions of the reduction technique 
was that s’ = O and thus 


constant 


C,/C; = s = 


since [from Eq. (14) ] 

5/6 ~ (pb/2V)/6 
This assumption may be investigated by inspection of 
Fig. 9, where it is seen that the maximum departure 
of any value from a constant taken as the average is 
approximately 3 per cent over 1.4 < 1/ < 5.0, excluding 
the results of Model 16. 

(d) Linearity of Roll Moment and pb/2V with Cant 
Angle. The analysis presented in this paper is based 
upon the assumption of a linear relationship between 
the lift coefficient at any point of the lifting surfaces and 
the angle of cant of the surface. This in turn implies a 
linearity of both the roll moment, C), and the steady- 
state effectiveness parameter, pb/2V, with angle of 
cant. Figs. 10 and 11 show the experimental values of 
roll moment and the effectiveness parameter angle as 
The 
resultant straight lines passing through the origin sub- 


functions of cant angle for two Mach Numbers. 


stantiate the original assumption and are an indication 


of the accuracy and reliability of the test results. 
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It should be noted that the values of effectiveness 
parameter extend over an extremely large range; the 
value obtained for 4° cant is equivalent to the spin tha 
would be obtained by firing a shell of a diameter equal 
to the fin span from a 1-30 twist gun, which is the 
standard rifling for many guns. 

(e) Effect of Fin Thickness on the Coefficients.—{y 
Figs. 12 and 13 are presented the roll damping moment 
and the roll moment coefficients as a function of Mach 
Number for the seven missiles with 16 per cent thick 
fins‘ and the 16 missiles with 8 per cent thick fins. [py 
all other respects the geometry of the missiles was 
identical within machining tolerances. Also plotted 
in these figures are the linearized theoretical curves as 
calculated in reference 1. Since the linearized theory 
becomes more applicable with thinner surfaces, the data 
substantiated the expectations of lower coefficients for 
the thinner surfaces. Note that the deviation of the 
coefficients from the linearized values for the thin-fin 
missiles is less than 12 per cent over the complete range 
of Mach Numbers. 

A comparison of the theoretical lift of an infinite. 
aspect-ratio single-wedge wing obtained from the 
linearized theory with that obtained from the exact 
theory using the two-dimensional expansion and shock 
equations results in a value approximately 9 per cerit 
lower for the linearized lift at a Mach Number of 2.5 
for the 8 per cent thick fins and 15 per cent lower. for 


the 16 per cent thick fins. The trend then in the two- 
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dimensional case is qualitatively similar to that in the 
three-dimensional case as would be expected. 

(f) Wing-Fin Interference.- The experimental values 
lor the aerodynamic coefficients, Cis, Cp, and (p,b 
21’) ‘6 obtained from the 16 per cent thickness surfaces 
for (1) the seven models with fins only over 1.6 < MJ < 
3.0, (2) the two models with wings only at J = 1.7 and 
1.95, and (3) the two models with wings and fins 
(aligned) at MJ = 1.7 and 1.95, are plotted in Figs. 14, 
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15, and 16. For the models with both wings and fins, 
only the four fins were canted (6 = 2°); as a result, 
the fis were in a region of downwash behind the wing 
resulting only from the angle of attack induced by the 
rolling velocity of the missile. Since the maximum roll 
rate was about 25° per ft. and since the distance be- 
tween the wing trailing edge and fin leading edge was 
about 3 in., the tail rotates a maximum of about 6° 
for a 3-in. flight path, which is quite insufficient to 
allow it to escape the downwash region of the wing. 
However, it is not located symmetrically with respect 
to this region. 

The data from the wing-fin models, Fig. 14, indicate 
an extremely large downwash effect on the damping 
coefficient for the tail surfaces. At a Mach Number 
of 1.95, the surfaces are only 14 per cent as effective as 
when acting alone, and at a \/J = 1.73 they contribute 
no damping at all. 

The large interference effect of the wings on the tail 
surfaces damping moment seems due to the fact that 
the greatest part of the moment is contributed by the 
surface area near the fin tip, and it is this region that is 
most seriously influenced by the downwash field of the 
wing. The geometry of the single-wedge airfoil results 
in a comparatively large wake, which also tends to de- 
crease the aerodynamic effectiveness of the fin sur- 
faces in this wake. 

The roll moment produced by the canted fin surfaces 
appears to be unaffected by the presence of the wing. 
This result seems reasonable, since the downwash field* 

* The turbulent rotational motion of the fluid also caused by 
the wing could conceivably influence the roll moment; however, 


this effect appears small from the data 
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of the wing arises solely as a result of the rolling velocity 
and as such can influence only the fin damping moment. 

One might interpret the small differences in the co- 
efficients for the wing-only models and for the fin-only 
models as a result of the variation in body-wing inter- 
ference for the different geometrical arrangements. 
Additional tests, however, would be necessary for ex- 
ploitation of this possibility. 

The data are only sufficient to indicate that the range 
technique is successfully adaptable to a study of wing- 
tail interference phenomena under different flight and 
design conditions. 
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Theoretical Pressure Recovery Through a 
Normal Shock in a Duct with Initial 
Boundary Layer 


GEORGE McLAFFERTY* 
United Aircraft Corporation 


ABSTRACT Figs. | and 2 for a '/;-power turbulent boundary-layer 
profile approaching the shock and a uniform profile 


The theoretical pressure recovery through a normal shock in a 
downstream of the shock. The theoretical pressure 


duct having an initial boundary layer has been calculated on the 
basis of continuity and conservation of momentum. The results recovery is plotted as a function of the ratio of displace- 
can be used to determine the maximum theoretical pressure re- ment thickness of the boundary laver to the semiheight 
covery in supersonic compressors, supersonic wind-tunnei dif- J : 

fusers, Supersonic ram-jet diffusers, and in any other case in which 
a wall boundary layer is present upstream of a normal shock. 
This calculation procedure has also been applied to diffusing 


passages having boundary-layer bleed systems to determine the ASSUMPTIONS (SEE TEXT AND FIGS. 3 AND 4) 
|. V7 POWER TURBULENT BOUNDARY LAYER 
2. UNIFORM PROFILE AFTER SHOCK 


and the boundary-layer passage. 
It is shown that the theoretical total pressure recovery through 3. NO WALL SKIN FRICTION THROUGH SHOCK 


anormal shock in a duct is reduced by an amount proportional 
to the displacement thickness of the boundary layer approaching ones ae -¢ om samme 
the shock for a wide range of approach boundary-layer profiles h 


and approach Mach Numbers 


of the passage. 





maximum theoretical pressure recovery in both the main stream 


NOTE: SEE FIG.6 FOR 848 
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1 ‘HE EXPERIMENTAL PERFORMANCE of supersonic at 10 
diffuser passages is usually evaluated on the basis ao 
of the relation of this performance to the isentropic S « pe 
performance with ideal approach conditions. In the = Ns 
application of diffusers to supersonic compressors and 4 Pal: 
° ; ine ; ze 
other cases where the inlet conditions are not ideal be- = i MACH NUMBER 
; pie sain - Ss APPROACHING SHOCK, M 
cause of the presence of initial boundary layers, the ulti- 2 = | _— 
mate theoretical performance must be modified by this 3 °’ Ss 750 
consideration to provide a realistic basis for perform- g _ 
ance evaluation. - 6 “Tel 
WwW 
DISCUSSION 8 5 2.00 
Pressure recovery through a normal shock with a = _ 
“s - if . . —— 
uniform profile approaching the shock is calculated by 2 ’ i ee 
slaccic. > r > ac1¢ > 1" 2eP ACL ” 
classical theory on the basis of the following three as- gf — -_. & | 
sumptions: (1) continuity of flow, (2) conservation of @ -3 
¢ ~ , Ps y 2 een 
momentum, and (3) conservation of energy. a mm 
These same assumptions were used to calculate the P .2 3.00- 
pressure recovery through a normal shock in a duct a 
with an initial boundary layer approaching the shock. - , 
The method of analysis is discussed in the Appendix. S 
Airflow and momentum of the initial boundary layer =< o 
were calculated from tabulated values of boundary- Py A 2 3 4 
layer displacement thickness and momentum thick- 3* -LA 
ae 
ness. The results of the analysis are presented in h MI-HEIGHT OF PASSA 
= OR - LAYER A 
Received March 17, 1952. Revised and received October 20, A OF PASSAGE 
1952 Fic. 1. Theoretical normal shock pressure recovery vs. displace- 
* Acronautical Research Engineer, Research Department ment thickness of initial boundary layer. 
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LAMINAR LAYER 
1.0 5 2.0 2.5 3.0 
MACH NUMBER APPROACHING SHOCK, Mo 
Fic. 8. Theoretical pressure recoveries for various approach 


boundary-layer profiles, calculated for 6*/A = 0.10 


Theoretical pressure recoveries for different assumed 
boundary-layer velocity distributions upstream of the 
shock are plotted in Fig. 3. It is shown that the theo- 
retical pressure recovery is relatively insensitive to 
changes in distribution from the !/;-power profile, even 
for the case of a laminar profile, where the difference in 
pressure recovery is only 2 per cent. The theoretical 
pressure recovery for velocity profiles downstream of 
the shock other than the assumed uniform profile is not 
changed if one type of momentum average total pres- 
sure is used to calculate the pressure recovery. Other 
total pressure averages provide pressure recoveries 
slightly greater or slightly less than the momentum av- 
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erage total pressure recovery. On the basis of these 
considerations, it is not felt that the results of the anal- 
ysis would be radically changed if the velocity pro- 
files upstream and downstream of the shock-boundary- 
layer interaction differed somewhat from the assumed 
profiles. 

It has been shown experimentally that a normal 
shock in the presence of a boundary layer occurs over a 
substantial length of passage. According to references | 
and 2, the length required to produce the complete 
static pressure rise through a normal shock is on the 
order of 6 to 15 hydraulic diameters. It was also noted 
in references 1 and 2 that a constant-area passage in the 
region of the normal shock allows a considerably greater 
pressure recovery than would be obtained if the normal 
shock occurs in a divergent passage. This result is ex- 
pected because a normal shock in a divergent passage 
occurs at an average area greater than the throat 
area and, hence, at an average effective Mach Numbe 
greater than the Mach Number at the throat. 

The effect of wall skin friction on the pressure re- 
covery through a normal shock may also be calculated 
using the methods discussed in the Appendix. The re- 
results of this analysis are given in Fig. 4 for a Mach 
Number of 2.0 approaching the shock. The wall skin 
friction is expressed as a fraction of the total stream 
thrust approaching the shock. An example of the use 
of Fig. 4 is calculated in the Appendix for a round pipe 
10 diameters in length. If the skin-friction force is 
calculated from the dynamic pressure of the flow ap- 
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THEORETICAL TOTAL PRESSURE RECOVERY, Rofo 





























A= 
~<-04 
.4 
3 SE EEE 
-2 
of 
° 





° a s 3 4 


_3* , BOUNDARY-LAYER DISPLACEMENT THICKNESS 

h SEMI-HEIGHT OF PASSA\ 

OR  BOUNDARY-LAYER DISPLACEMENT AREA 
AREA OF PASSAGE 





Fic. 4. Example of effect of wall friction forces on the pres- 
sure recovery through a normal shock occurring over a length of 
passage 
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S.., ARY-LAYER P_HEIGHT 


§ BOUNDARY-LAYER THICKNESS 
Fic. 5. Example of calculated pressure recoveries in inlet system 
with boundary-layer bleed. 


proaching the shock, the ratio of wall friction to total 
inlet momentum in this example is 0.034, which reduces 
the maximum total pressure recovery about 4 per cent. 
The skin-friction coefficient used in this example is high 
for two reasons. First, the average dynamic pressure 
is substantially less than the dynamic pressure of the 
stream approaching the shock. Second, the skin-fric- 
tion coefficient would be reduced because the boundary 
layer is separated in the region of the shock. From this 
reasoning it is believed that the total pressure recovery 
would be reduced on the order of only 1 or 2 per cent 
by wall skin friction in this example. 

It is often undesirable to have a constant-area throat 
of sufficient length to obtain the complete static pres- 
sure rise through a normal shock. The length of the 
required throat can be reduced by removing the bound- 
ary layer approaching the shock. Removing the 
boundary layer also has the effect of increasing the 
pressure recovery through the shock in the main pas 
sage. The maximum pressure recovery in the boundary- 
layer passage and the maximum pressure recovery in 
the main stream can be calculated as a function of the 
amount of boundary layer bled off by the same type of 
analysis used to obtain the results in Figs. 1 and 2. 
The results of a sample calculation are given in Fig. 5. 
In this example, bleeding off the complete boundary 
layer provided a maximum theoretical pressure re- 
covery of 72 per cent in the main passage as compared 
to 62 per cent with no boundary-layer bleed. The 
theoretical pressure recovery in the boundary-layer 
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passage for complete removal of main-passage boundary 
layer was 52 per cent. 

It can be shown theoretically that isolation of. the 
high-velocity main stream and the low-velocity bound- 
ary-layer stream will result in less overall increase in 
entropy. This results from elimination of losses due to 
mixing the main-stream high total pressure air and the 
boundary-layer low total pressure air. 

Calculated values of the ratio of boundary-layer dis- 
placement thickness to boundary-layer total thick- 
ness are shown in Fig. 6 as a function of Mach Number 
for several different boundary-layer profiles. The turbu- 
lent boundary-layer profile values were obtained from 
reference 3, and the laminar boundary-layer profile 
values were obtained by integration of the boundary- 
layer equations. 

The results shown in Figs. 1 through 6 have a large 
range of application to all cases where the pressure 
rise through a normal shock is affected by an initial 
boundary layer, such as in supersonic compressors, in 
supersonic wind-tunnel diffusers, and in most types of 


supersonic ram-jet diffusers. 





TURBULENT BOUNDARY LAYER VALUES FROM REF. 3 













































































a Lannea scunpany Laven 4 
\2 
eee 
5 
TURBULENT 
¢ BOUNDARY LAYER 
&O)" 
4 4 
ait e 
| 
° if 
~ 
Pa Net) 
.2 7 - 
mw 
4 
7 
ie) 


° ! 2 3 


MACH NUMBER, Mo 


Calculated boundary-layer displacement thickness 
vs. Mach Number 





172 JOURNAL OF THE 


REFERENCES 


1 Neumann, E. P., and Lustwerk, F., Supersonic Diffusers for 
Wind Tunnels, M.1.T. Meteor Report No. 13, November, 1947. 

2 Neumann, E. P., and Lustwerk, F., High-Efficiency Super- 
sonic Diffusers, M.1.T. Meteor Report No. 56, June, 1950. 

3’ Tucker, Maurice, Approximate Calculation of Turbulent 
Boundary-Layer Development in Compressible Flow, N.A.C.A. 
T.N. No. 2337, April, 1951 


Appendix 


List OF SYMBOLS 


a = height of boundary-layer scoop (see Fig. 5) 
A = cross-sectional area 
Ay = surface area for determination of skin-friction drag 
A/Ax = isentropic area ratio, a function of Mach Number only 
‘ (y + 1)/( 1) 
Cc = airflowconstant = @/ z. ( : ) ci ca 
\ RTz\7 + 1 
Cr = skin-friction coefficient (see example in Appendix ) 
D = hydraulic passage diameter (see example in Appen- 
dix) 
F = stream thrust [see Eq. (1)] 
h = semiheight of passage 
Kp = stream thrust parameter [see Eq. (2)] 
Ky = airflow parameter [see Eq. (10)] 
l = length of passage required to obtain complete static 
pressure rise through normal shock (see example in 
Appendix) 
M = Mach Number 
N = velocity profile parameter in turbulent boundary- 
layer equation, V/V» = (y/6)! ('N 
Ps = static pressure 
Pr = total pressure 
R = gas constant for air 
Tr = total temperature 
V = velocity 
W = airflow 
y = distance from wall of passage 
¥ = ratio of specific heats 
6 = boundary-layer total thickness 
6* = boundary-layer displacement thickness, 
1 tg 
= (poVo — pV) dy 
pol 0 j 
A = wall-friction force expressed as a fraction of upstream 
stream thrust [see Eq. (5)] 
p = density 
6 = boundary-layer momentum thickness, 
1 
za pV( Vo — V) dv 
po Vo? I ° 
Subscripts 
0 = midstream conditions upstream of shock-wave boundary- 
layer interaction 
1 = average conditions upstream of shock-wave boundary- 


layer interaction 
2 = average conditions downstream of shock-wave boundary- 
layer interaction 


In a nonuniform stream with constant static pres- 
sure in a symmetrical passage with semiheight 4 and 
unit area, the stream thrust is defined as 
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1 h 
PF = Py i (1 + yM?) dy | 
h Jo : 
1 *h 
Kr = I (1 + yM?) dy (9 
h Jo F . 


Eq. (1) becomes 


Let 


F = PskKr (3 
The stream thrust before and after a shock must be 
equal if there is no wall skin friction. 
F, = F, (4 
If the wall skin friction is a fraction, A, of the upstream 
stream thrust, Eq. (4) becomes 
Fil — A) = F (5 
Substituting quantities from Eq. (3), 
Ps, Ke(1 — A) = Ps, Kr (6 
or 
Ps,/Ps, = (Kri/Kr2)(1 — A) (7 
In the same manner, the airflow in the passage of unit 


area is 


; : l *h P, 
W = dy (8 
h. 0 (A A *) : 
or 
W l *h l 
* — Ps, dy (9 
C h 0 ( Ps Pr )( A A * ) a 
where C is a function of 77, R, and y. 
Let 
: l *h l 
Kw = dy (10 
h Jo (Ps/Pr) (A/Ax) ~ 
so that Eq. (9) becomes 
W/C = PsKy (11 


The airflow before and after the shock must be equal, 


W, = W, (12) 


From Eq. (11), this becomes 


Ps,Kw, = Ps Ky (15 


Ps, Ps, ae Kw, Ky, (16 
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Combining Eqs. (7) and (16), 


(Kr, Kr.) (l1 — A) = Ky, Ky, (17) 


Aw./Kep, = [1/0 —- A)] (Ay, Kr,) (1S) 


It is convenient to express Aw, and KA,, in terms of 
Ky, and Kr,, which are evaluated assuming that the 
mid-stream conditions in the nonuniform flow approach- 
ing the shock exist across the complete stream. Eq. 


18) can be written 


= (19) 


an Kr, Kr, 


Kw, l Ky Ky, Kw, 
hw ie 


The quantities Aw,/Aw, and Ap,/Ar, are determined 
by the flow conditions approaching the shock. They 
are expressed in terms of more conventional boundary- 
layer parameters by the following derivation. By 


definition, 


Kp, l “hh e ) 
— = / (1 + y Al*) dy} / (1 + ydAlo?) = (20) 
Kr, h JS0 : 


75 *h 
/ (1 + yM?) dy + / (1 + yM?) dy 
of 0 of © 


Ke es (?] 
K, hl + yM,2) - 
: 6+ foe dy + (h — 6)\(1 + yM,?) 
Ky 0 . 

"= . (22) 
Kp Ah(l + y M0") 


The quantity f°y.\/2 dy can be obtained by proper 
manipulation of the expressions for boundary-layer 
displacement thickness, 6*, and boundary-layer mo- 
mentum thickness, @. 


| 6 
f+ 9= - | (peVo — pV) dy + 
po | 0 0 


“ 


=a pV(Vo — V) dy (23) 


po | o7 J/0 


*5 V 
rtom f (1- "0 ) ay + 
0 po | 0 
*6 V |”? 
| (f ~ rr) dy (24) 
0 po | 0 po | o~ 
“ V V le 
e+6= J (: ee - + - -—— : -,) dy (25) 
0 po | 0 py | 0 po | 0~ 
75 V2 
/ (: ne -,) dy (26) 
0 pil 0" 


. 


te 
s— | ace (27) 
Jo poVo? 


pV? = PsyM?* (28) 


* +0 = 


- 


Itcan be shown that 


Substituting in Eq. (27), 


THROUGH A 


NORMAL SHOCK l 


3 Poy M? 
a 4+0=6- = (29) 
> Ps.yMe 


Since the static pressure is assumed constant across the 
passage 


* 


] 6 
-6* +0=6- M 7: yi? (30) 
ably 0 


J yM? = yM,? [6 — 8* — 0 (31) 
J0 


Substituting Eq. (31) in Eq. (22), 


or 


6 + yAly?(6 — 6* — 0) + (h — 6)(1 + yAo?) 
hil + yMo?) 


(32) 


Kr, 
Kr, . 
h(1 + yMo?) — 6 — byMo? +6 + y¥Mo?(6 — 6* — 8) 


hl + yM,2) 


Kr yMy?(6* + 8) 
Kr, 7 h(l + yM,?) 


Kr 6 yM? (= *) : 
— =< f= + (35) 
Kr, h(1l + yMo?) \6 6 


From the definition of boundary-layer displacement 
thickness, 6*, 


5* /h) (36) 


Kw,/AKw, = 1 —.( 
or 
Kw,/Kw, = 1 — (6/h)(6*/5) (37) 
Substitute Eqs. (35) and (37) in Eq. (19), 
Kw, 


Kr, - 
(39 
( l ) Kw, h 6 
1 — A/ Ky, 6 Mv 6* 0 
“lia ewe a 
h (1 + yMo?) \6 6 
If a uniform profile is assumed downstream of the shock, 
Kw,/Kr, calculated from Eq. (38) determines the 


Mach Number downstream of the shock, Js, which in 
From Eq. (16), 


(3S) 


turn determines Aw, and Ps,/ Pr,,. 


Ps, » Kw, _ Kw, [1 — (6*/h)] 
Ps, Kw, — Kw, 


(39) 


The following tautology exists: 


Pr, Ps, Ps, Ps./Pr 
Oe sh (40) 


Pr, Ps, Ps, Ps,/Pr. 
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Ps,/ Ps, is equal to unity from the original assumption of 
1 0 | ] 
uniform static pressure in the upstream flow. Sub- 
stituting Eq. (39) in Eq. (40), 


Pr, _ Ps,/Pr, Kw, ( = * 


= (41 
Pr Ps, Pr, Kw, h 


The theoretical total pressure recoveries through a 
normal shock presented in the present report were cal- 
culated using Eqs. (38) and (41). The wall friction 
term, A, was assumed equal to zero except for the cal- 
culation of the results shown in Fig. 4. As an example 
of the use of Fig. 4, calculate A for flow in a round pipe 
at Mach Number 2.00. 

_ wall friction force (49) 

~ total inlet stream thrust z 

a C,As(¥/2)Ps,Mo’ (4 
Ps, ACL + yMo?) 


C,, wDI(y/2) Mo? 
A= (44) 
(9/4)D7(1 + yMo?) 


MARCH, 1953 


A = 1.70(1/D)C,, (45 


The wall skin-friction coefficient in the region of the 
shock 1s difficult to determine but would probably be 
low because of the separated flow due to the adverse 
pressure gradient through the shock. Also, the friction 
force should be based on an average dynamic pres. 
sure through the shock instead of the dynamic pressure 
upstream of the shock. However, as an example, as- 
sume //D = 10 and C,, = 0.002. The skin-friction 
coefficient of 0.002 was obtained from a standard chart 
of turbulent flat-plate skin-friction coefficients for a 
Mach Number of 2.0 and a Reynolds Number of 20 x 
10°. From Eq. (45), 


A = 1.70(10) (0.002) = 0.034 (46 


From Fig. 4, this amount of wall skin friction will re- 
duce the total pressure recovery approximately 4 per 
cent. The 4 per cent figure is probably an upper limit 
according to the reasoning stated previously, and the 
actual reduction in total pressure recovery is probably 
more nearly on the order of 1—2 per cent in this example. 
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Forced Convection Heat Transfer from 
Horizontal Cylinders in a Rarefied Gas 


F. M. SAUER* anp R. M. DRAKE, JR.7 
Unwersity of California at Berkeley 


ABSTRACT 


An approximate analytical treatment of the energy equation is 
given which yields a specification for the average heat-transfer 
coeficients from horizontal circular cylinders to gases flowing at 
high Mach Numbers and low Reynolds Numbers. The results 
are given for continuum and slip-flow regions in comparison with 


available experimental data. 


NOMENCLATURE 


a = thermal diffusivity, sq.ft. per sec. 

D = diameter of cylinder, ft. 

h = heat transfer coefficient, local, B.t.u./hr. ft.2 °F. 

Jv = Bessel’s functions of first kind 

Kv = Bessel’s functions of second kind 

k = thermal conductivity of the gas, B.t.u.-ft./hr. ft.2 °F. 

L = reciprocal length proportional to molecular mean 
free path 

M = Mach Number (dimensionless ) 

n = anarbitrary constant 

Nu = Nusselt Number, (D/k), (dimensionless ) 

r,0,5 = spatial coordinates, ft *or rad 

r = radius of cylinder, D/2, ft 

Re = Reynolds Number, )D/v (dimensionless ) 

t = temperature at any point, °F 

t = temperature of the free stream, °F. 

ty = temperature of the cylinder surface, °F 

u,v,’ = velocity components in the r, 6, = directions, ft. per 
sec 

V = any constant velocity in the @ direction, ft. per sec 

V = free-stream velocity, ft. per sec 

a = thermal accommodation coefficient (dimensionless) 


= ratio of specific heats of constant pressure and con 


stant volume (dimensionless ) 


Ml = absolute viscosity, Ibs. sec./ft.? 

\ = molecular mean free path, ft 

p = density, lbs. sec.?/ft 

v = temperature ratio, (¢ — tye)/(to — te) (dimensionless 
o = Prandtl Number (dimensionless ) 

T = time, sec 


INTRODUCTION 


i es SUBJECT OF HEAT TRANSFER from bodies in a 
high-speed rarefied gas stream has become impor- 
tant in several fields of investigation. If a gas is suffi- 
ciently rarefied, itis clear that a continuum treatment fails 
because the mean distance between molecular encoun- 


ters becomes comparable with significant body dimen- 

Received June 17, 1952. Revised and resubmitted November 
17, 1952 

* Formerly Assistant Professor Mechanical Engineering; now 
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sions, and the effect of the molecular structure on the 
flow and heat transfer must be considered. The na- 
ture of the rarefaction effects with regard to the aero- 
dynamic effects has been examined by Tsien! and with 
regard to heat transfer by Drake and Kane.” Although 
the requirements attendant to the design of high-speed 
high-altitude missiles have provided the stimulus for 
most investigations, the results are of interest for other 
systems. Rarefaction effects may occur at fairly ordi- 
nary pressures if the body dimensions are small, as for 
instance in a hot-wire anemometer or thermoelement*® 
to be used in high-speed flows. 

The use of hot wires as turbulence measuring devices 
requires extremely thin wires in order to reduce the 
thermal lag. These wires operate at Reynolds Num- 
bers of the order of 10—1,000. 
the molecular mean free path are comparable, and at 
supersonic speeds, the flow around the wire is in some 


The wire diameter and 


cases slip flow. 

It is the purpose of this paper to introduce approxi- 
mate solutions for the specification of the average heat- 
transfer coefficients for horizontal cylinders appropriate 
to laminar flow in the continuum and slip-flow regions 
in order to estimate the effect of the temperature jump 
associated with slip flow on the heat transfer. 

A comparison between continuum and slip-flow heat 
transfer is possible if considerations are made on the 
basis of a continuum formulation of the energy equation 
wherein the rarefaction effects, as specified by the tem- 
perature jump condition, are relegated to the boundary 
values of the problem.’ In particular, if the mean 
tangential velocity at the cylinder surface is taken as 
an arbitrary constant the energy equation may be re- 
duced to a heat conduction equation which can be 
solved to yield a specification of the average heat trans- 
fer coefficient for continuum flow conditions. The re- 
sults so obtained may be made to correspond to existing 
experimental data for average heat-transfer coeflicients 
from horizontal cylinders® '! in continuum flow by spec- 
ifying the magnitude of the mean tangential velocity. 
Using the value of the mean tangential velocity thus 
determined, another solution of the heat conduction 
equation for temperature jump boundary conditions 
yields an average heat-transfer coefficient for the slip 
flow case which may be compared with the continuum 
results for effect of the temperature jump on the heat 


transfer. 
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Fic. 1. Coordinate system. 


ANALYSIS 


Consider an infinite cylinder placed normal to a uni- 
form gas stream. Neglecting the terms relating to com- 
pressibility and dissipation of mechanical energy, the 
energy equation may be written in cylindrical coordi- 
nates, Fig. 1, 


of of , vot | Of _ 
oa ve” “eS 
joer 1dr 1 O% ON 
. Yar: + ror rae? t dstf (1) 


The solution of the problem requires steady state, 
0/07 = O, two-dimensional flow, w = Of/d0z = 0, and 
the additional assumptions: 

(a) Conduction of heat in @ direction is negligible, 
0*t/00? = 0. 

(b) The radial velocity is everywhere zero, u = 0, 
and the angular velocity is a constant, v = V. 

(c) The surface temperature of the cylinder is a con- 
stant. 

Introducing a 
(t — by)/(b 
sumptions, Eq. (1) becomes 


V 00 ford 
= a- 
r 00 lor? 


dimensionless temperature, # = 


t,,), and making use of the above as- 


l od 


») 
r or f -) 


Boundary conditions applicable to continuum or slip 
flows may be considered. In the continuum-flow case 
the fluid adjacent to the cylinder surface has the same 
temperature as the surface; thus, 3 = 0, while the 
fluid far away from the cylinder has a temperature /, 
and 3d = 1. The point 6 = 0, r = m™ provides some 
Q and 


|. This is best seen if the problem is visualized in the 


difficulty, inasmuch as 3 has two values there 


SCIENCES—MARCH, 


1953 
r-6 plane. However, at distances along @ = 0 somewhat 
greater than 7;, the cylinder radius, the oncoming flow 
has the same characteristics as the free-stream flow: 
thus, for r > r; along 6 = 0; 38 = 1. For the slip-flow 
case, the above arguments apply except at the surface, 
r = nr, where there exists a temperature discontinuity 
between the surface and the adjacent gas which is pro- 
portional to the heat rate. Making this change, the 
boundary conditions may be written as: 


Continuum Flow 


v = 0, r=nNn, <9 <s 
v = il, , > Me, 6=0 
Slip Flow 
Ov 
= 16, r= fr, O<0< fr 
or 
g =i], o> he 6=0 


where 


l 2—- 
= 1.996 ( : ) ( *) 
L y+ 1 a o 


SOLUTION ON THE DIFFERENTIAL EQUATION 


Introducing a new independent variable y = 7, 
Eq. (|) has the form 
ov  1day  1e_ | 7 
oy? yoy K0O#@ 


where A = a/4V and the boundary conditions become: 
Continuum Flow 


¢o=0, y= 
v=1, y> ri, 


Vr, O<0<9 


Slip Flow 


Ow , rae 
Ov =i, y=VnO<0< 8 
i=l y> on 6=0 


where & = 2L Yn. 

Eq. (3) is identical with the equation for the transient 
radial conduction of heat into an infinite region bounded 
The slip-flow bound- 

“radiation”’ bound- 


internally by a circular cylinder. 

ary condition is identical with the 
ary condition. The solutions are well known and have 
been found by 


Carslaw and Jaeger’ among others. 


the Laplace transform method by 
Each case will be 


presented separately. 


CONTINUUM FLOW 


The solution of Eq. (2), transforms from the solution of Eq. (3), is the following: 


v= — 
T 


where uw is an arbitrary parameter of integration. 


2 ee Jy(un r) Vo(urn i) = Yo(uy r) Jun r;) du (4) 
a J(u Vr) + Yoru Vr) u 
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The local heat-transfer coefficient at the surface is given by 


h = k (O8/Or),-,, (5) 
Performing the necessary differentiation, the local Nusselt Number becomes 
hD t ” aid du (6) 
= ; 
k wr Jo Jur) + Yoru Yr) u 
Let U =u Vr, and x = Aé 7, = a6/2VD; Eq. (6) becomes 
i ed dU 4 100,1:x)) (6a) 
Sa. ee ars - = — (0,1 5x) 6a 
k wT Jo Jer(l )+ Vyr(V) l 1T- 


Jaeger and Clark® give an excellent tabulation of the function /(0,1;x). The same reference gives a series valid 
for small values of x, (vx << 1). 
oe? a. t te. tis, Pi as act _ 
1(0,1;") =~ | (wx) .~ + x — O.1469x" + 0.263x? — 0.315x" + 0.53607 +... (7) 
~ qT < 
The average heat-transfer coefficient will be calculated over the region 0 < 6 < x. 


hD l In 
( ) “s { T(0,1;x) dx (&) 


k 


“re 
where x, corresponds to the value of x at 6 = 7, 
For small values of x,, using Eq. (7), 


1 


(AD /k)ay = 1.128x,>'? + 0.500 — 0.0941x,'" + 0.0625x, — 0.0588x,°" + 0.0677x,2 . . (9) 


For the region x < 0.1, (kD /k),, was calculated from Eq. (9); for x > 0.1, the average heat-transfer coefficient 
was found by numerical integration of Eq. (8) by means of the tabulated values of Jaeger and Clark. 


The value of |’ has been left arbitrary up to this 


‘ : . ee ge : TABLE | 
pont. Let = I./nand o = 0.70 for air; 
; hD 
nwa nT 2.24n ’ ( k ). me 
x. = = = a 
4 2V,D 2 Reo Re 0.01 12.8 672.0 
0.10 4.04 67.2 
, . . 0 1.57 6.72 
The calculated solution may be compared with the z 0 i 4 3 36 
recommended curve of \[cAdains‘ for forced convection 10.0 0.740 0.672 
of aj . eae eae BS eres 100.0 0.431 0.0672 
fair flowing normal to cylinders. If = 3, the curves 200 0 0 379 0 0336 


are nearly coincident in the region of Re less than 100 


as seen on Fig. 2. The calculated values of Re with _ 7 ; 
Eq. (9a) may be compared with the empirical Eq. (4) 


n = 3 are shown in Table 1. 7 . . 
of McAdams based on the experiments of many investi- 


For Re > 300 and o representative of air, Eq. (9) re- 
duces to gators, 


(hD/k),,. = 0.500 + 0.435Re"” (9a) (hD/k) ay. = 0.32 + 0.43Re°”” 


Sire Flow 
The solution of Eq. (2) with the temperature jump boundary condition is 


[ut (ur/ry) +E Ivu/ ry)? + [uVilurr) + EVo(uY ny) |? u 


YE [ Ji(uvir) [UV (ur) + & Volur/r)| — VYolurr) [usi(uvn) + ESo(uv 4) Je Ku’ du (10) 
7 0 


Note that Eq. (4) is a special case of Eq. (10)—-+1.e. > &. 
Differentiating Eq. (10) and substituting into Eq. (5), the local Nusselt Number is found to be 
hD 1¢2 P isiadit du 1 
= (11) 
k i J luJi(uv/r,) + E Sour) |? + [uVilurr,) + EVo(uv 7) |? u 


Asbefore, let U = uy/r,. After some rearrangement Eq. (11) can be written in the dimensionless form, 
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TABLE 2 
DL 0.2 0.5 1 2.5 5.0 10 
Re (hD/R )ax 
0.672 0.160 0.304 0.485 0.583 0.651 0.679 
1.344 0.167 0.329 0.488 0.690 0.780 0.824 
6.72 0.180 0.392 0.650 1.03 1.24 1.35 
67.2 0.192 0.453 0.837 
672 0.932 2.12 
6,720 4.38 7.15 
67,200 4.81 9.30 
x 0.2 0.5 1.0 2.5 5.0 10.0 
> x2 


l 
DL{I(DL,x)] (11a) 


hD _ DL / e | 
k J0 | -_ F 1 , | dl 


"(DL) ‘AUT (U) + SDE)" JU) | 4 Hu UYU) + "(DL)" Vol | - i 


where x = AO/r, = a0/2VD and DL = const. (D/X). 
For small values of x, the integral 


~ rl? idee 
, il 
I(p, q, x) = - Es — -—— —— * - (12) 
Jo (pus (U) + qJo(u) + [pPUVWU) + @¥(U))? l 
isapproximated by the series 
I(p, gq, x) = (w*/2pq) [1 — 2p-qa- x? + (1/2)p-? q(p + 2g)u +... |] 


For p = (mw 2) (DL) and q = 2/2 (DL)’’, Eq. (12) becomes 
I(DL, x) = [1 — (2/r'*)x"* + (1/2) (DL) (1 + 2DL)x +... | (13) 
For larger values of x, /(p,g;x) is plotted in reference 7, p. 226, and also in reference 8, p. 283, from which curve 
values were read and replotted to larger scale. 
As before, the average heat-transfer coefficient is found in the region 0 <@< 7. For small values of x,, 


hD | is 5: es 
( F ) = 22,112 — —- (DL) x, °>+ (1+ 2DL)x, +... (14) 
4 »T 


av . 
. 


Forx, > 1, the average value was found by graphical and numerical integration of the curves of Eq. (8). Values for 
= 1/3 Vy) are shown in Table 2. 


Assigning values of a = 0.80, y = 1.4, and o = COMPARISON WITH EXPERIMENT 


70 for air, - , ; , ' 
rhe results of the analysis as shown graphically in 


DT. = D/25 > Figs. 2 and 3 may be compared with experimental data 

of heat transfer from cylinders in cross flow as reported 

For conditions peculiar to rarefied gases, there exists a py Kovasznay and Térmarck* and Stalder, Goodwin 
relationship between the molecular mean free path to wail Creager.” The data of reference 3 were taken 
dy dimension ratio and the Mach and Reynolds under conditions of supersonic flow at normal pressures 
Numbers, which is for air, with very small wires, and the data in reference 10 were 
F taken in supersonic flow under very low pressures with 

A/D = 1.5 (M/Re) : - , 

somewhat larger cylinders. The data so obtained cover 
thus, a range 0.07 < MW/+~/Re < 5.5. On the basis of the 
presently defined flow region limits, these data are con- 

DL = 0.267 Re/M sidered to extend from the near slip region into the free- 
molecule flow region.! The experimental data are 

A plot of (AD /k),, 
tter is shown on Fig. 2. Curves with Mach Number as_ wherein the velocity and density are taken as the free- 


versus (I)D/v) with DL asa param- shown in terms of Nusselt and Reynolds Numbers, 


aparameter were obtained by cross plotting from Fig. stream values while the viscosity and thermal conduc- 
-,and these are shown on Fig. 3 as the desired solutions tivity are evaluated at the stagnation temperature. 
lor the average heat-transfer coefficients. The data are considered herein on that same basis, al- 
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2.0 | l] a; * Terror | 
la Nu 20.267 tRe 
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0.06} © 0.05! 1.65 NS 
— @ 0.080 1.60 
0.04} 9 0.126 1.60 to 2.65 aN 
| @ 0.126 1.90 to 2.30 \ 
NOTE (1) SOLID POINTS REPRESENT DATA TAKEN 1 
0.02 IN UNIVERSITY OF CALIFORNIA. 
(2) SHADED AREA IN LOCALE OF KOVASZNAY 
AND TORMARK DATA. 
0.0l ! hediebi tA RS) ! LaediidiA AAS | 
0.0! 0.02 0.04 0.080! 0.2.04 0.81.0 2 4 6 810 
MACH 
YREYNOLDS 
Fic. 4. Comparison of analysis and experiment for heat transfer from horizontal circular cylinders to rarefied gas (air). 


though it is realized that this particular selection of 
property values tends to suppress the Mach Number 
effect. 

In order to compare the analysis with the experi- 
mental data, Fig. 3 was replotted as shown in Fig. 4. 
In this plot the intercept for 17/1/Re = 
the Hilpert continuum experimental data as represented 
by Eq. (4) in reference (6) and by Eq. (9a) herein. 
As the /, ~ Re ratio increases for greater rarefaction, 
this effect is shown on the Nu/+/ Re ratio. The limit- 
ing solution is for the case of infinite 1//+/ Re ratio. 
It may be seen that the analysis predicts correctly the 
trend of the experimental data, although some question 
remains with regard to the representative magnitudes, 
since for the most part the data of Kovasznay and 
Tormarck are higher than the values predicted by the 
theory. The the fluid property values 
would greatly influence the correspondence between 
the experimental data and the analytical results in any 


0 represents 


selection of 


case. 


SUMMARY 


An analytical treatment has been presented which 


yields a specification of the average heat-transfer co- 


efficients for right circular cylinders arranged _hori- 
zontally in a gas stream for flow conditions character- 
istic of both continuum and slip flow conditions. This 
solution of the energy equation was effected for bound- 
ary conditions incorporating a temperature discon- 
tinuity between the cylinder surface and the gas by 
virtue of a simplifying assumption concerning the na- 
ture of the free-stream velocity distribution around the 
cylinder. The solution is approximate in this respect. 
The solution obtained closely represents the con- 
tinuum heat-transfer results for horizontal cylinders to 
air; and, when compared with experimental data for 
horizontal cylinders in the slip-flow region, it was seen 
to predict the trend of the data satisfactorily, though 
some question remains regarding the representative 
magnitudes. The disagreement here may be in part a 
result of the selection of the fluid property values as 
were used in the reduction of the experimental data. 
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Inelastic Instability and Incremental Theories 
of Plasticity’ 


E. T. ONATF ann D. C. DRUCKER? 


Brown University 


SUMMARY 


A most troublesome paradox has existed for a number of years 
with respect to buckling in the plastic range. Theoretical con 
siderations and all direct experimental evidence show conclu 
sively that an incremental or flow type of mathematical theory 
of plasticity is valid. However, the results of plastic buckling 
tests are well correlated by a simple total or deformation theory! 
and bear no resemblance to published predictions of incremental 
theory .* 6 

The suggestion was made?’ that initial imperfections of shape 
or loading might well explain this most peculiar result. How 
ever, subsequent investigations by several authors seem to have 
given the impression that excessively large imperfections would 
be needed and that the answer would be overly sensitive to the 
magnitude of such imperfections. ® * 9 

It is the purpose of this paper to demonstrate that extremely 
small, and therefore unavoidable, imperfections of shape do ac 
count for the paradox in a simple manner. The buckling load is 
shown to be extremely insensitive to the amount of imperfection 
The example chosen is a simplified version of the long rectangular 
plate hinged along one edge and free on the other under uni 
form compressive stress at the ends. This is the equivalent of 
the case of the cruciform column, which has been so disturbing 
in the past because incremental theory applied to a perfect cruci 


form column did lead to an entirely incorrect result.° 


STRESS-STRAIN RELATION 


THE incremental stress-strain 


7 ‘WE EXACT FORM OF 
relation used is not especially significant. Any 
reasonable isotropic loading function will lead to 


essentially the same result, so that the simplest one, 


the Mises function J2, will be chosen. For plastic 
straining, d/, > 0 and 
: (1 + pv) 
dey, = de;;° + de? do 
E 
v é . 
do,,6;; + F(J2)s;dJo (1) 
where 
J (1 /2is.4 
Si; = oi; — (1/3) 04,6 


vis Poisson's ratio, /2 is Young’s modulus, and 4,; is 
the Kronecker delta. 


indicates summation. 


Repetition of an index in a term 
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ISI 


If the only stresses present are a normal stress o and 
a Shearing stress r on the same plane plus its cross- 
shear, the two important components of the relation 


(1) may be written 


do ’ 4 oda 
de = E + F(J2) , 9 — a rdr (2) 
dy Qdéne. = 
2(1 + v) : oda 
kK dr+ F(Js) 47 + rdr (3) 
where de is the normal strain increment in the o 


direction, dy is the (engineering) shear strain incre- 
ment corresponding to 7, do and dr are the increments 


of normal stress and shearing stress, and /2(1 + vy) 
is the shear modulus G, 

J. = (1/3)e? + 7° 

dJ»y = 2[(eda/3) + rd} (4) 


F(J,) can be determined from a simple tension test 


for which 


de = (da/E) + F(J2)(4/9) o*do 
de 
sh es 
or 
F( J.) = (3/4/2)[(U/E,) — C1, E)] (5) 


where /, is the ordinary tangent modulus and JJ» is 
given by Eq. (4) in the special case of combined normal 


and shearing stress only. 


PREVIOUS DIFFICULTY 


Although the paradox described in the Summary 
occurs for all problems studied, it appears at its worst 
in the torsional buckling of a cruciform in the plastic 
If a cruciform, Fig. 1, is loaded uniformly in 
The kine- 
matics is particularly simple if the ends are not re- 
strained (Fig. 2). Each straight line parallel to the 
The state of 


range. 
compression, failure by twisting is likely. 


axis becomes a helix. stress and strain 
at 


straight column at the instant of buckling, 


is the same each cross section. For a_ perfectly 


shearing 
stress and strain are added to the existing state of simple 
compression, point B of the curve / constant, Fig. 


3. Whether or not the Shanley concept of increasing 








182 JOURNAL OF THE AERONAUTICAL SCIENCES—MARCH, 1953 
































—— J 
I ZD_ 
y/ 
7 ' 
' 
\ 
| 
' 
L 

1 
1 
\ 

ee 




















FIG. 2 























INELASTIC INSTABILITY AND TFTHEORIECS OF PLASTICiItTyYy 183 








de B ale 
hex! 
































Fig. J 
t 
#| } 
| } 
a connenenity ghsomms 
| 
i 
| 
FIG. 4 
o 
| 
o| E,* & 
E 
1 
—e€ 
rts. § 











han! y OR 6h 


FIG. 6 


Twist 


OF .THE CRUCIFORM WITHOUT AXIAL STRESS 








the compression is followed and no matter what smooth 
symmetric loading function is assumed, the normal to 
the loading surface and therefore the plastic strain 
increment vector at B points along the negative o- 
axis. This means that all incremental theories of 
plasticity which do not have a corner at B predict 
a purely elastic response to a shearing stress incre- 
ment from point B regardless of how large a normal 
stress increment is applied simultaneously. The criti- 
cal compressive stress a, predicted in this way is there- 


fore given by 

o1,0 = GCO (6) 
or 

ao. = Gh?/b (7) 


where @ is the angle of twist per unit length, G is the 
elastic shear modulus, GC is the torsional rigidity, /, 
is the polar moment of inertia, / is the thickness of 
each leg, and 6 is the width, Fig. 1. Simple J. de- 
formation theory,*® on the other hand, predicts G,h?/b’, 
where G, is the secant modulus in shear. This latter 
result is in excellent agreement with experiment.‘ 

It is well recognized now that a physically valid 
theory must be of the incremental type. Furthermore, 
direct tests on circular tubes under compression and 
torsion verify the elastic response to torsion.'”!! 
Therefore, when the agreement of buckling tests with 
deformation theory was so prominent in technical 
discussion, the suggestion was advanced that initial 
imperfection might well provide the explanation of 
the paradox.’ More confusion than enlightenment 
seems to have resulted from several subsequent in- 
vestigations that were interpreted as proving that 
excessively large imperfections would be needed and 
that the answer would be overly sensitive to the 
magnitude of such imperfections. 


A SIMPLIFIED MODEL 


Considerably clearer understanding of the physical 
behavior involved is obtained by a simplified model 
and a simplified stress-strain diagram. The cruciform 
cross section (kh < 6) will be taken as a thin shell of 
constant thickness ¢, Fig. 4, instead of a solid plate; 
but the cross-sectional shape will still be supposedly 
maintained under twist. As a consequence, the shear 
stress on the cross section follows the contour and is 
constant in magnitude, and the shear strain y is given 
by 

y = hd (8) 


where @ is the angle of twist per unit length. 
The elastic buckling load, Eq. {6), is 
GC S8Gh*bt ’ 
a” = = —— (approximately ) 
le 8tb3/3 
3Gh? /b? (9) 


II 
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under the assumption 4 < 6. Deformation theory 


would substitute the secant modulus for the elastic 
modulus. 

As the paradox and its explanation are not linked 
to any particular stress-strain relation, the simplified 
two sloping lines of Fig. 5 will be chosen for clarity 
of exposition. 


A QUALITATIVE PICTURE OF THE BUCKLING PROCESS 


The basic physical picture depends upon the shear 
stress-shear strain relation for the test piece, which 
may be translated directly into the moment-angle of 
twist per unit length relation. If there is no axial 
force acting on the shell cruciform, the 1/-@ relation 
is simply the stress-strain curve to a diiferent scale, 


Fig. 6. The moment is 
M, Srhbt (10) 
and y 6h if the cruciform is perfect to start. If the 


cruciform is initially twisted by a small amount, 4) pet 
unit length, the lines are just shifted to the right as 
shown. The presence of appreciable axial compressive 
stress changes the picture completely. 
the magnitude of the compressive stress by o and 


Designating 


employing Wagner's concept of the stress following the 
fiber, the axial force produces a twisting moment 


Mos ol 0 (11) 


The external twisting moment needed to maintain the 
twist is now 


M MM, VV. (12) 


Should 1/. equal or exceed 1/;, no external moment 
is needed, and the column will fail by twisting. Equat 
ing ./, to M/. for an initially untwisted column gives 
the elastic result of Eq. (6) for the critical load. It 
can also be readily seen that the maximum load for a 
completely elastic column is not appreciably influenced 
by initial twist (location of point A, Fig. 7). Ai, 
increases more than /), for a given increase in 6, until 
o reaches the elastic buckling stress. 

If the shell column is loaded in axial compression 
only, so that there is no external moment, .1/; M/s and 
(approximately ) 


Srhbt = AM; 


a(S 3)tb°6 


a) (") M, . 
oc = = (1.9) 
(. fa (S/3)th°¢ 


or 


The important physical point here is that the shear 


ing stress and the twisting moment increase as oa 
increases and also as @ increases. From Eq. (13), 


3h (dr rdé 
ado 3 ( —_ : 
fh 4 4- 
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INITIAL IMPERFECTION HAS LITTLE EFFECT ON THE 
ULTIMATE LOAD OF A COMPLETELY ELASTIC COLUMN WHICH 
FAILS BY TWISTING 


It is somewhat more convenient to work in + and 
notation rather than J1/, and 6, which are proportional, 
respectively, Eqs. (10) and (8). 

Substituting Eqs. (5), (13), and (14) in the shear 
strain relation, Eq. (3), which applies in the plastic 
range, and rearranging 


(1/n) |(71/Gy) 


dr 1 fi 
| (1D) 


Gdy nN 1+ 3(7/0)? 


where 7 (3G /f,) + 1 (3G) is approximately 
the ratio of the elastic shear modulus to the tangent 


shear modulus. In terms of J/, and @, 


ar Gy AM / SGh*bt6 


Consider again the idealized stress-strain diagram 
of Fig. 5. Eg. (15) is valid as soon as the yield point 
is reached, which for all practical purposes requires 
a = o* because 7 is zero or very small. The enormous 


importance of extremely small imperfections in shap 


is immediately apparent when the plastic line is flat 
n is infinite. Upon reaching yield, the slope of the 
r-y or j/,-6 curve immediately goes from the elastic 


/ 


slope to the value + y or 1/,/6. No matter how small 


the initial imperfection, no further increase in @ 1S 
required to continue twisting. In fact, the 7/o tern 


makes the process 


in the denominator of 


Eq. (15 
the column will twist to complete collapse 


instable;: 
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ELASTIC INSTABILITY AN 


INE 


=o*. In Fig. 7, if of = (1/2)o,, o will increase 


A to point P and then stay almost con 


at 0 
from point 
stant, as Shown by the dashed line. 

If nm is reasonably large but not infinite, say 20, the 
slope of 7 G vs. y or M,/SGhbt vs. 60h, which 1s unity 
in the elastic range, drops to a value that is larger 
than the slope of the radial line from the origin by 
bout 1 40 in Fig. 8. Therefore, some increase in o 
is required to continue the twisting of the column 
Until the + o term in the denominator of Eq. (15) 
becomes important, the + G vs. y curve is fairly close 
to a logarithmic spiral centered at the origin. How 


ever, the tangent to the spiral is so nearly a radial 
line that a very sinall increase in o over the value o* 
is sufficient to produce a large enough twist for ro 
to take over and produce collapse. 

\s an example, if o* is one-half of the elastic buckling 


stress 


g (3. 2)(Gh?/b?) 


o*1,0 1Ghhte 


the equation of the spiral is approximately 

} rye’? (16) 
where ¢ is the angular coordinate. This means that 
before o exceeds o* by 20 per cent, the radius will be 
e' or 50 times the radius when plastic action began. 
Initial twists that are extremely small compared with 
maximum possible elastic twists therefore require 
stresses only a little higher than o* to cause collapse. 
Furthermore, large differences in initial imperfection 
change the maximum or buckling stress only a little. 
LOAD 


QUANTITATIVE SOLUTION FOR THE MAXIMUM 


AND THE LOAD Twist RELATION 


The preceding discussion gives a physical picture 
of plastic buckling and explains why the buckling 
stress cannot exceed the yield strength of structural 
metals by a large amount. Detailed calculations are 
obtained in more useful form, however, from an equa 


tion in o rather than 7. In terms of y Oh, 


l ' 3 (4) f 
do n ” \3n7] 3h? G : 


(17) 


y T 0 


dy 3 ( =) E b> on 
ke  \ BR, , 


when the material is in the plastic range, or 


6 — 4% k ») + "| = (18 
dw l+ W nN 


Where bh? 3h?)0? and a,‘ (3h? /b7)G. 


The equation for the elastic range is obtained by 


setting 7 |. <Any initial twist will be magnified 


is the compressive stress is increased from 0 to @ 


For convenience however, the twist at o* will be 
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FIG. 8 


INITIAL IMPERFECTION IS THE ANSWER 





designated by @ and the corresponding value of 


by W. The smaller true initial twist can be com 
puted if desired, but it is not especially important 
here. 


The solution of Eq. (17) may be written as 
GA" s1 + &Y 
o o* t 
WV | 4 7 
‘- 
(1 4 c) n 
0, J Vv ; 


2n errr 4 


1)/2np(1 —2n)/2n 7. 
dk 


(n—1)/2n 
V) 


Because V is small compared with unity and Wy is 
extremely small, a good approximation to Eq. (19) 


1S 
1/2n 1/2n 
Wo WV, m 
. o* ) a ae 
Vv Vv (20) 
g¢= 
1+ [(n 1)/2n|(V — W) 
If %) = 0, o = o,'; but, for n = 20 if VW) + O, then 
(%/W)'’*" is close to unity for extremely large ratios 


of V/W%. Large twists therefore occur for small 
increase of o over o*. The denominator of Eq. (20 
then becomes appreciable, and o decreases as the 
twist continues to increase. 

Fig. 9 is a plot of ¢/ o* vs. W for various Wy) under the 
same assumption used for illustration in Figs. 7 and 8, 


as " = 1/2. 20 has been chosen 


oO 0, 
as reasonable for an aluminum alloy. 


The value of n 
Interpretation 
of the results is aided by noting that W = b? 12/°(2b8 
of the relative angular rotation of two cross sections 
a distance apart equal to the approximate overall 
depth of section, 20, is 


ee , 1/2 
260 3.46(h byw rad. 


20W' * deg. 
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FIG. 9 
2b@ 1S THE ANGLE OF TWIST IN A LENGTH 2b ATO: o* 


THE ACTUAL INITIAL TWIST IS THEREFORE ONLY HALF AS 
LARGE 





for h/b close to 1/10. An initial W) = 10°° or 266) = 
3.8 min. is small indeed and certainly will exist. The 
actual initial twist is only half of the magnified value 
of 3.8 min. As shown in Fig. 9, a change in @ by a 
factor of 10 produces a change of less than 10 per cent 
in the critical load. Comparison with deformation 
theory is also included to show that initial imperfection 
can, and in all likelihood does, account for the experi- 
mental data. 


CONCLUDING REMARKS 


Many practical problems involving maximum com 
pressive loading of actual columns and plates in the 
linear elastic range can be treated as static buckling 
problems for perfect elements. <A_ static stability 
criterion may not be useful, however, for nonlinear 
problems. This seems to be especially true when’ the 


system is nonconservative as well. Apparently a 
solution of the boundary value problem for an initially 
imperfect member is often required to obtain the 
physically proper result. In some cases, however, a 
dynamic type of criterion may be needed to analyze 
the effect of lateral or other disturbance. 

In this paper it has been shown bv use of a simplified 


cruciform and a simplified stress-strain relation that 


1953 


SCIENCES 


MARCH, 


incremental plasticity 


results when unavoidable initial imperfection of shape 


theory does lead to proper 
is considered. Detailed consideration of the actya| 
cruciform and a more realistic curved stress-strain 
relation is left for a later paper. 
that failure by 
strength of the material is appreciably exceeded. 


However, it is clear 
twisting will occur before the yield 


Although the predictions of deformation and incre. 
mental theory now agree fairly well, it should be 
noted that it is a somewhat involved function of the 
tangent modulus which is significant rather than the 
secant modulus of deformation theory. Special and 
possibly unrealizable stress-strain curves therefore do 
exist for which predictions of the two theories would 
differ widely. 
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Thermal Distributions in Jeffery-Hamel 
Flows Between Nonparallel Plane Walls’ 


KNOX MILLSAPS? ann KARL POHLHAUSEN? 
flight Research Laboratory, U.S.AF. 


ABSTRACT ° O(ru) Ov 
: Sa = = 0 (1.1) 
The authors give the exact solution for the thermal distribu or ov 
tions for the steady laminar flow of a viscous incompressible fluid 
between nonparallel plane walls held at a constant temperature momentum, 
The velocity profiles are determined with the aid of Jacobian id 
i ; , Ou v Ou v 1 Op 
elliptic functions by using the Jeffery-Hamel solution of the = = 
hydrodynamic problem. It is shown that in this special case or r OW r por 
the energy equation giving the temperature profiles can be re 2 dv " 
‘ - . : . . ; - = Ji 
duced to an ordinary linear differential equation with variable y (vs — — ) (1.2) 
coefficients. After the introduction of dimensionless parameters, r° OW oi 
numerical solutions are given for diverging and converging chan- 
; ; ‘ Switec and momentum, 
nels with total openings of 10° for the possible combinations of 
-_ Yy Y . ‘ Fae De- ’ . 
three Reynolds Numbers and five Prandtl Numbers. ov P v Ov 5 uv 1 dp 
u =— 
or r Ov r pr OW 
INTRODUCTION 
v2 2 Ou v 13 
iii aati ’ ae os v Vv — ) 
I’ IS WELL KNOWN that the flow of a viscous incom- v+ r2 Oe y2 (19 
pressible fluid between plane parallel plates, the 
so-called plane Poiseuille flow, is governed by a parab- where 
olic velocity distribution. This flow is easily shown 22 12 1 2? 
to be an exact solution of the Navier-Stokes equations rave a > sae 
| or? r or r? Ow? 


because of the fact that the nonlinear terms vanish 
identically. For the case of nonparallel plane walls 
where such fortuitous circumstances do not occur, the 
exact solution has been found by Jeffery! and, independ- 


and p and v denote the pressure and kinematic viscosity, 
respectively. 
If the flow is purely radial, then one may put 


ently, by Hamel.’ . . u = vF(d)/r, —- (14) 
It is the purpose of this paper to give the exact solu- 

where F(#) is termed the dimensionless velocity pro- 

file, and one then notes that the terms in Eq. (1.1) 


vanish identically and that Eqs. (1.2) and (1.3) be- 


tion and a numerical discussion of the thermal dis- 
tribution for the Jeffery-Hamel flow under the assump- 


tions of constant fluid density and constant wall tem- 
come, respectively, 


perature. 
, y* FF? | Op yp?" 
(1) FORMULATION OF THE or oe dir (1.5) 
- “i p Or ‘ia 
HYDRODYNAMICAL PROBLEM 
; ; ‘ 1 Op 2v?F’ 
In polar coordinates (7, 3) the equations for the ve- QO= — t a (1.6) 
i ; 3 . 
locity components (u, v) for the steady (0/of = 0) pr ov “ 
flow of an incompressible (density p = constant) fluid yf Eq. (1.5) is multiplied by r*y~*, a rearrangement 
are: continuity, gives 
Paper presented before the Eighth International Congress of (r*?/pv?) (Op/Or) = F? + F" (1.7) 


Theoretical and Applied Mechanics held at Istanbul, Turkey, 
in August, 1952. Received July, 1952. 
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A partial integration of Eq. (1.6) with respect to 3 gives 


r gives 


the entire period required for this study. ‘ 2 
: : : Sei r® Op or 

+ Consultant in Applied Mechanics = —4F + (1.9) 
v’p Or av" 
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One sees that p as given in Eq. (1.8) is a proper func- 
tional relation, provided that S(r) is chosen so that 
S’r5 is a constant; hence, 


S(r) = Kpv?/2r? (1.10) 


where K is a constant. 

It is now clear from Eqs. (1.9) and (1.10) that Eqs. 
(1.5) and (1.6) will be simultaneously satisfied and 
that a solution to the flow between nonparallel plane 
walls, } = +a, will be obtained, provided that 


F?+4F + F°+ KK =0° (1.11) 
with the boundary conditions, 
F(+a) = 0 (1.12) 


expressing the imperviousness of the walls. 

Thus the partial differential equations that formulate 
the conservation of mass and momentum have been 
reduced to an ordinary nonlinear differential equation 
of the second order whose solutions must satisfy in the 
physical problem under consideration two boundary 
conditions of vanishing velocity at the walls. 


(2) SOLUTION OF THE HYDRODYNAMIC PROBLEM 
After multiplying Eq. (1.11) by 6F’, an obvious in- 
tegration yields . 


2F* + 12F? + 3(F’)? + 6KF + 2H = 0 (2.1) 


where 2// denotes the integration constant. Solving 
Eq. (2.1) for (F’)*, one gets 


(F’)? = (2/3) (H7 — 3KF — 6F? — F*) (2.2) 


()" P dF 
v= = s ee es) 
2 / (7 — 3KF — 6F° — F*) 


where the integration limits must be determined later 
from the boundary condition. Thus, it is clear that # is 
given by an elliptic integral. An evaluation is given 
by Hamel? in terms of the Weierstrassian functions, 
but, if numerical results are desired, it seems to be 
more convenient to use Jacobian elliptic functions be- 
cause of the present availability of a suitable table.'* 
The unavailability of an extensive table of elliptic 
functions appears to be the principal reason that 

although the incomplete theoretical solution has been 
known since 1915—numerically calculated velocity 
profiles for assigned Reynolds Numbers have been pub- 
lished only for singularly selected values.*'°''  Fur- 
thermore, Goldstein'’ remarks, ‘“The discussion is 
hardly yet complete,’ but, since completing this paper, 
the authors’ attention has been called to a penetrating 
theoretical analysis of the problem by Rosenhead.!! 
However, a direct comparison between his results and 


or 


those of the present paper is rather difficult because of 
the different definitions of the Reynolds Number, the 
use of different notations for elliptic functions, and his 
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SOURCE 


SINK 


Schematic diagram of flows in converging and diverging 
channels. 


Fic. 1. 


use of the convergence and divergence of the channel 
as a parameter. 

Assuming that the trinomial in Eq. (2.3) has real 
roots e¢; 2 e@. 2 e; (the remaining case of one real root 
is discussed later), one can factor the denominator so 
that 


(’)" : dF 
les i =a, (24 
27 J (la — F) (e — F) (es — F)] 


where 


—5 (2.5 


€) + €2 + 6s 


and the roots are related to the integration constants 


Hand K by 
€\€2 + C3 + Qe; = 3K (2.6 


eee; = IT (2.7 

It is now necessary to distinguish two different types 
of flow that are described by Eq. (2.4). The first type 
is the flow in diverging channels in which it is assumed 
that a line source with a positive efflux velocity is pres 
ent at the origin of the diverging channel walls; the 
second type is the flow in converging channels in which 
it is assumed that a line sink with a negative influx 
velocity is located at the intersection of the converging 
Schematic diagrams of these flows are 

It is convenient to follow the treat- 


channel walls. 

given in Fig. 1. 
ment by Hamel,’ and its excellent review by Gold 
stein,'’ in the subsequent discussion of the flows in 


diverging and converging channels. 


Case (A) Flows in Diverging Channels 
It is clear that in this case F 2 0, e; 2 0,e; < e& < 
II > 0,and0 < F & e, and that 


3\ "7 e1 dF 
gg = : (2.8) 
2 ak F we F) (F — eo) (F — es) | 


The reduction of this integral to a Jacobian elliptic 


function is given by 


where 
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THERMAL DISTRIBUTIONS 


k? = (€; — e2)/(e; — e3) (2.10) 


6 th a= F ws 
j= sn! ,& 2.11) 
€1 — @3 a — 


by m, an inversion yields 


Thus, 


After denoting |(e: — é3)/6] 


F = e, — (€; — es) sn*(md, k) (2.12) 


The roots may be expressed in terms of m and k, and 


these expressions are 


é, = 2[m*(1 + k*) — 1] ] 
é. = 2[m*(1 — 2k?) — 1] (2.13) 
e; = 2[m*(k? — 2) — 1] 


Now, if the velocity along the axial stream line, 3 
= (0), is prescribed to be ms and if the Reynolds Number, 
Ry = Fy = wor/v, is introduced (the case of asymmetric 
back-flow patterns requires a slight modification), then 
the fact that sz(0, k) = 0 enables one to conclude from 


Egs. (2.12) and (2.13) that 


Ro = e = 2[m?(1 + Rk?) — 1] (2.14) 
or 

m = {{1 + (Ro/2)]/(1 + R*%)}” (2.15) 
Thus, the velocity profile is given by 

F(3) = Ry — 6m?k?sn2(md, k) (2.16) 


and, since the boundary condition (1.12) must be satis; 
fied, it follows that 


Ry = 6m*k?sn*(ma, k) (2.17) 


It is clear that the elimination of m from Eqs. (2.15) 
and (2.17) yields an equation determining ; this 


transcendental equation is 


1+ k Ro/2)] 
: ~ == 5g? At =! | Qa, et (2.18) 
3k[1 + (2/ Ro) ] Lite f 


which is most easily solved by graphical means. 


Thus, if the divergence a and the Reynolds Number 
Ry are specified, k may be found as a root of the trans- 
cendental Eq. (2.18), and m may be obtained from 
Eq. (2.15). The velocity profile may then be calcu- 
lated from Eq. (2.16). 

As was previously noted in Eq. (2.4), it is pos- 
sible that the trinomial in the denominator on the right 
side of Eq. (2.3) has one real root e; and that the re 
maining two roots, é: and e;, are conjugate complex 
numbers, 


@e=art di, é=a — bhi 


In this case Eq. (2.4) becomes 


Bi (*) i dF 
Z / “= F) (F? — 2aF + a” i b?) | 


(2.19) 
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Now, since e; + é: + e; = —6, it follows that 
a= —[3 + (e,/2)] (2.20) 
If one introduces the following notations, 
é,;> — 2ae, + a? + 5? = gq? (2.21) 
Gqta-—e@ 3+ q+ (8e, 2) F 
= =k? (2.22) 
2q 2q 


then a familiar transformation of elliptic integrals gives 


*e dF 
I [(e. — F) (F? — 2aF + a? + 6%]? 


| = F 
7, cn 1 (! é: + ‘ e) (2.23) 
q”* g+e—F 


The combination of Eqs. (2.19) and (2.23) yields 
g-atF 2q\"” 
- = CN ( : vk 
q + = F 3 
If, as in the case with real roots, one prescribes that 
F(O) = Ro, then, since cn(0, k) = 1, it follows that 


(2.24) 


qQ = Ry (2.25) 
Since the boundary conditions require F(a) = 0, one 
has 
gq — Ro 2q\ "7 
t “= cn ( ') a, k (2.26) 
q + Ry 3 


The substitution of Eq. (2.25) into Eq. (2.22) gives 
g = 3[1 + (Ro/2)]/(2k? — 1) (2.27) 

The determination of k can now be made by introduc 

ing g from Eq. (2.27) into Eq. (2.26) and solving the 


resulting transcendental] equation, 
6 4 Rod —_ $k?) 2 of Ry : 
: — = cn . a,k (2.28) 
6 + Ro(1 + 4k?) 2k? — 1 
If a solution for k is known, then g may be calculated 
from Eq. (2.27). Subsequently, the velocity profile 
may be calculated from Eq. (2.24), and it is found to be 
Ry — gt (Ro + q) cn [(2q 3)' v, k| 


F(d) = . 
1 + cn [(2q/3) ” 8, k| 


(2.29) 


At this point it is interesting to consider several par- 
ticular examples of diverging flows. 

Example (1): Flows with Low Reynolds Number in 
Channels of Small Divergence.—In this instance the 
elliptic function sv in Eq. (2.16) reduces to its argu 
ment, and the velocity profile is given by 

F(d) = Ro — 6m'tk?3? 
The boundary condition yields 


Ry — 


6m'k? a? 


hence, 


F(d) = Rol{l — (8?/a?)] 
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which is the familiar parabolic velocity distribution for 
plane Poiseuille flow. 

Example (2): Flows Having Unitary Elliptic Modu- 
If k = 1, it follows from Eq. (2.15) that 


? + =] * & 
m = = 
4 2 


Moreover, it is well known that for unitary modulus 
the elliptic function su(md, 1) degenerates into tanh 


larity. 


md so that 


Ry “9 
r tanh? 
; ss tanh? md? = 2 
F(d) = Ro l —_ ie ee ~ I =~ i/, 
tanh? ma F " 
tanh” 


where from Eq. (2.18) 


4 2\""7]2? «5.232 
= | tanh ; ( ) | = = 
a” 3 a” 


Flows with Vanishing Inclination at 
It is easily seen from Eq. (2.2) that 


Ry = 


Example (3): 
the Channel Walls. 


of (*) ; (*e**) : 
F’'(a) = = 
3 3 


and, if this is to vanish, then at least one of the roots 
must be zero. It is convenient to take es = 0, and with 
this choice Eq. (2.13) shows that 


m = [1/(1 — 2k?)]'” 


Since 0 < k < 1, a critical comparison with Eq. (2.15) 


enables one to conclude that 


Re? wx '/s, m = (Ro/3)'” 


In this example it is clear that Eq. (2.12) reduces to 
sn?(ma, k) = 1 


Since sv is a periodic function, this equation is satisfied 


when 
ma = nK(k) ~ a(Ro/3)'” 


where A(k) denotes the complete elliptic integral of 
the first kind and » is an odd integer. When k? = 
1/5, K(k) = 1.85407, and therefore for n = 1 one has 


a(R)? = 3?K = 3.211 


From these considerations it follows that the first Reyn- 
olds Number having zero slope at the walls is given by 


Ro = (3.211/a)? = 2Ry 


where Ry denotes the Reynolds Number of the flow 
having unitary elliptic modularity. The velocity pro- 
file is easily calculated by appropriate substitution in 
Eq. (2.16). : 


Example (4): Flows with Large Reynolds Numbers. 


As the Reynolds Number of the diverging flow becomes 
larger than the Reynolds Number corresponding to 
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zero slope at the walls, the computation of the corre. 
sponding velocity profiles for the larger Reynolds Num. 
bers shows that a region of back flow appears along 
each wall. As the Reynolds Number increases to stil] 
larger values, positive jets, in addition to the back-flow 
This process con- 
tinues with an increasing number of back-flow regions 


regions, will appear along each wall. 


and positive jets appearing in a symmetrical distriby. 
tion about the main flow. However, it should be noted 
that, since between each forward jet and each back. 
flow region there is a point in the profile with zero ve- 
locity, at this point a wall could be placed leading to 
asymmetrical velocity profiles with the possible num- 
ber of asymmetrical configurations increasing in the 
directly obvious manner with increasing Reynolds 
Number. 

Moreover, at sufficiently large Reynolds Numbers it 
is found that the transcendental Eq. (2.18) has mul- 
tiple roots. Each member of the set of solutions cor- 
responding to a definite Reynolds Number leads to a 
possible symmetrical profile with back-flow regions and 
positive jets. As the Reynolds Numbers increase, so 
do the number of solutions of Eq. (2.18), although it 
should be noted that the number remains finite. The 
asymmetrical profiles may be constructed by the pre- 
viously outlined procedure. 


Case (B) Flows in Converging Channels 


In this case it is easily seen that F < 0 and that the 
roots are real with e, > Oand0 2 e. > e3. Further, one 
notes that 0 2 F 2 e, and, consequently, from Eq. 
(2.4) it follows that 


i (3) or dF 
in 2 / [(e — F) (F — €s) (F cs es) | - 


(2.30 


An application of a standard reduction formula for 


. | (2.31 


F = e3 + (e. — e3) dn~*(md, k) (2:3 


elliptic integrals gives 


l Co — €3 
o = —dn : 
m F — @;, 


The reduction for F gives 


or, in terms of m and k, one has 


F = 2[m*(k? — 2) —1] + 


6m?(1 — k?) dn-?(md, k) (2. 


2.33 

If one introduces the Reynolds Number R> as the ve 
locity along the axial streamline, one sees that 

Ro = Fo = & (2.34 

m? = [1 + (Ro/2)]/(1 — 2k?) (2.39 

The imposition of the boundary condition F(+a) = 9 


yields a transcendental determination of k, 


dn?(ma, k) = 1 — Ro/es 
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TABLE 1 


Diverging Flows 


Ro k? m? e} 
684 1.000 171.5 684 
1,342 0.4989 448.3 1,342 
5,000 0.745 1,393.3 5,000 
or, since dn®x = 1 — k*sn*x, a substitution for e; from 


Eq. (2.13) gives 


k Ro(l _s 2 2”) pe 
sn’'(ma, Rk) = PRET ze 2.37) 
2k?[3k2 — 3 + (Ro/2)(R? — 2)] 





For extremely large Reynolds Numbers, Eq. (2.37) re- 


duces to 
sn*(ma, k) = (2.38) 


If the value of k which is determined by solving either 
Eq. (2.37) or (2.38) is inserted into Eq. (2.35) and if 
this value of m and the value of k are substituted into 
Eq. (2.33), the velocity profile for a prescribed Reyn- 


olds Number may easily be calculated. 


It is interesting to consider several particular ex- 
amples of flow in converging channels. 


Example (1): Flows at Low Reynolds Numbers in 
Channels of Smail Convergence. 


the series development for dn(x, k) shows that 


For small Ry and a, 


(2.39) 


dn(x, k) = 1 


dn-?(x, k) = 1 + k?x? (2.40) 


The substitution of approximation (2.39) into Eq. 
(2.36) gives 


e; = Ro/m?a*k? (2.41) 


and the further substitutions of Eq. (2.34) and approxi- 
mations (2.40) and (2.41) into Eq. (2.32) give 


F = R)[1 (2.42) 


— (3°/a?)] 


which is the familiar parabolic velocity distribution of 
plane Poiseuille flow. 


Example (2): Flows with Large Reynolds Numbers. 


For large Ro, a critical study of Eq. (2.35) and approxi 


mation (2.38) shows that k ~ | and m? = R»/2; hence 
é = é;, and the integral in Eq. (2.30) can then be evalu- 
ated in terms of hyperbolic functions. If this evalu- 
ation is made, one obtains a familiar result due to Pohl- 


hausen :4 


Ro\"” 
F= R, 3 tanh? | (*) (a — v) + 
2\' ) 
tanh (3) | _ 2¢ (2.43) 


Profiles with Back Flows. 
cendental Eq. (2.37) that determines the elliptic modu- 


Example (3): The trans- 


larity of the converging flows has multiple roots—in 


€2 C3 F'(a) 
—345 —345 —7,367 
0 —1,348 0 
— 1,646 —3,360 135.8 X 108 














A * 


Fic. 2. Graphical solution of the transcendental equation for k. 


Sis 





The 


velocity profiles 
The convex profiles are 
1,342, and 
the concave profiles are diverging flows with 
R = 5,000, (B): R = 1,342, and(C): R = 684. 


Fic. 3. Numerical values for the 
broken line is the Poiseuille parabola. 
converging flows with (A): R = 5,000,(B): R = 
(C): R = 684; 
(A): 


fact, a denumerable infinitude. The largest value of k 
leads to the normal symmetrical profile, and each of the 
remaining values leads to a symmetrical profile with 
jets and back-flow regions. 
that asymmetrical profiles may be constructed for con- 


It should also be noted 


verging in exactly the same manner as in the case of 
diverging flows. This infinitude of possible profiles 


was first noted by Rosenhead."! 
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° TABLE 2 
Converging Flows 
R Rk? m? e} Co e3 F"( a) 
—684 0. 9457 382.5 1,486.5 —684 —808.5 2,341 X 10 
—1,342 0.9835 692.8 2,746 —1,342 —1,410 5,887 X 10 
—5,000 0.9998 2,500 9,997 —5,000 —5,003 1,083 & 102 
transcendental determination of k—.e., Eq. (2.18 
One sees that the sn*(ma, k) curve intersects the curve 
1 + k? 
R 3k7[1 + (2/Ro) | 
1o*4 E , 
at k = | when Ry = 684. From this value, k decreases 
0 (as Ry increases) until k*? reaches approximately 0.5 
4ot4 when Ry = 1,342. For still larger Reynolds Numbers, 
-5° o* 5° k again increases. The pertinent data for the three 
vw 





Fic. 4. A numerical example for a diverging flow with R = 


15,000 showing the determination of the three possible values for 


k and the complete set of associated profiles. 











Fic. 5. A numerical example for a converging flow with R = 
684 showing the determination of k and the associated profiles 
for the first three of the infinite number of possible solutions. 


Numerical Examples 


In addition to the general examples previously con- 
sidered in this section, it is instructive to consider a 
few specific examples of the procedure for obtaining 
numerical values for the velocity profiles. In the fol- 
lowing numerical examples, a total angle of 10° has 
been chosen for the divergence or convergence of the 
fixed walls. 

With the prescribed divergence and for the diverging 
solution, the Reynolds Number that corresponds to 
unitary elliptic modularity is 684. One also finds that 
the Reynolds Number corresponding to a profile with 
zero slope at the walls is 1,342. As a third case, a ve 
locity profile that has back-flow phenomena is given; 
a sufficiently large Reynolds Number to illustrate these 
2 is a graph showing the 


Fig. 


phenomena 1s 5,000. 





chosen examples are contained in Table 1. 


Fig. 3 gives the velocity profiles F/ Ry or, equivalently, 
u/uy. An obvious inspection shows that each of these 
profiles is more sharply curved in the middle of the 
channel than the Poiseuille parabola, which is valid for 
parallel flows. It is to be noted that all of the profiles 
have inflection points. 
highest Reynolds Number is particularly interesting, 
Further, one 


The velocity profile for the 


since it exhibits two back-flow regions. 
sees that F vanishes at four points across the profile 
and that walls may be placed at any of these points 
giving flows with two symmetrical back-flow regions 
with a total divergence of 10° or an asymmetrical flow 
(or its mirror image) with one back-flow region with a 
total divergence of 6.9°. 

In order that comparisons may be made, velocity 
profiles for flows in converging channels have been 
calculated for the same set of Reynolds Numbers and 
the same amount of convergence as divergence in the 
numerical examples that were previously discussed. 
The pertinent data for these converging flows are given 
in Table 2. The corresponding velocity profiles are 
given in Fig. 3. It should be noted that the profiles 
are flatter in the middle than the Poiseuille parabola and 
that the profiles have no inflection points. 

As numerical examples of flows with a multiplicity 
of possible profiles, the three possible profiles for a di- 
verging flow with Reynolds Number of 15,000 are given 
in Fig. 4, together with a graph showing the transcen- 
dental determination of k in this case, and the first three 
profiles of the possible infinitude for the converging 
flow with a Reynolds Number of 684 are given in Fig. 9, 
together with the associated determination of k for this 
case. It should be noted that in some cases of con 
verging channels considered here the total flow 1s an 
efflux. 

In concluding this section on the hydrodynamic prob 
lem of the flow between nonparallel walls, many 1™ 
portant conjectures can be made and many interesting 
questions can be asked on the bearing of the existence 


of jets and back-flow regions on the problen of hydro 
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dynamic stability and on the possibilities of studying 
turbulence by examining in greater detail the multi- 
plicity of pe ssible solutions of the Navier-Stokes equa- 


tions. 


2) FORMULATION OF THE THERMODYNAMIC PROBLEM 


The continuity Eq. (1.1) and the two Eqs. (1.2) and 
1.3) for the momenta formulate the hydrodynamic 
problem of the flow through diverging or converging 
channels. To determine the thermal distribution 
within the channels it is necessary to solve the equa- 
tion resulting from the substitution of the solution of 
the hydrodynamic problem into the energy equation. 
Expressing the differential equation in terms of the pre- 
viously used symbolism and denoting the temperature 
of the fluid by /, the energy equation is known to be 


Of 4 
(u >» 


oo) s (e+ | Of 1 l ~~) 4 
roo) p \Or* r or ¢* OF 


\ 


v (26, + e557 + eva) (a.1) 


where c, is the specific heat expressed in mechanical 
mass units, x is the thermal conductivity, and the com- 


ponents of the strain tensor are given by 


€;r = Ou/Or 

— 1 Ov u 

nasil rower (3.2) 
1 Ou Ov v 

er5 = _ 


r Ov or r 


In a closely analogous manner to the analysis used 
by Jeffery and Hamel in obtaining the hydrodynamic 
solution, it is found that the partial differential equa- 
tion, Eq. (3.1), can be reduced to an ordinary differ 
ential equation by using the Jeffery-Hamel transfor- 
mation, Eq. (1.4), together with a transformation 

t = [g(8)/r?] + ty (3.3) 
where /,, is the constant temperature of the walls. This 


reduction yields 


K » ” vO » 
—2vc,F¢ (4g + 2”) + v(4F? + FF) (84) 
p 
md, upon introducing the Prandtl Number, o 
pve, kK, and the dimensionless function, G te/r; 2 
lurther reduction gives 
G” + (4 4 2a Fk )G o(4F? +4 F’?) 0 (3.9 


lf one wishes to determine the thermal distributions 


within a channel of total opening 2a, then the solution 


must satisfy the boundary conditions 


if Eq. BS 


G = () (3.6 
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and it is clear that the dimensionless temperature, 
T = (¢,/uo”)t, is given by 


T = (G/R,*) + Ty (3.7) 


where Ry is the previously defined Reynolds Number 
and 7\, is the dimensionless wall temperature. In the 
case of symmetrical profiles, it is advantageous in per- 
forming numerical calculations to replace the boundary 
conditions (3.6) by 


G'(0) = 0, G(a) = 0 (3.8) 


(4) SOLUTIONS OF THE THERMODYNAMIC PROBLEM 


Eq. (3.5) is an ordinary linear differential equation 
of the second order in normal form; it is‘an example of 
the generalized Lamé equation that is usually treated 
in studies on ellipsoidal harmonics. However, since 
the homogeneous part has a variable coefficient con- 
taining the square of an elliptic function and since 
the excitation is a complicated combination of elliptic 
functions, the analytical integration is not simple even 
in the exceptionally solvable instances due principally 
to Lamé, Hermite, Halphen.'? Although the 
analytical solution is known only for certain cases, 
treated to illustrate 


and 


some special instances are 
nontrivial solutions of the energy equation in finite 
form. 

In accordance with the general theory of linear dif- 
ferential equations, it is possible to obtain the solution 
for the homogeneous equation and then to treat the 
problem of the inhomogeneous equation by the method 
of the variation of parameters. Following this pro- 
cedure, it is convenient to consider the case of diverging 
flows in some detail and then to indicate the result for 
converging flows. It is easily seen from Eqs. (2.16) 
and (3.5) that the homogeneous Lamé equation for 


diverging flows is 


G” + [4 + 120m?*k?sn*(ma, k) — 


l2om*k?sn*(md, k) |G = 0 (4.1 


The most general instance of this equation which has 
been treated is the case when 
(4.2 


12am? = n(n + 1) 


where 7 is a positive integer. In this case the solution 
given by Hermite!’ is 
" H(i + V,) : | 
G A ll J exp IZ (d,) ; + 
= | Of 
jHI(d v 


6 IL | 


AG, + Ba, 


IZ (3 


exp 
Od 


provided that the solutions are not doubly periodic, 
that H(#), O(9), and Z(d) are the principal Jacobian 
functions, and that the 3; are determined from the ” + | 


independent equations comprised ia the system 
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Fic. 6. The thermal distributions in a converging channel 
of total opening of 10° for a flow with R = 684 for parametric 
values of the Prandtl Number. Only the normal solution is 
shown. 











Fic. 7. The thermal distributions in a converging channel of 
total opening of 10° for a flow with R = 1,342 for parametric 
values of the Prandtl Number. Only the normal solution is 


shown. 
ore = 1S. ...5 
"" snd, cnd, dnd + snd, cnd; dnd; ») 
aoe tl sn?3, — snd; | 
n 2 
p cnd, dsb, -- (4.4) 
r= 1 | 


n 


i ns?d, + 4+ n(n + 1)R?sn?(ma, k) = 0 { 
1 


7 


where 5,’ f(7) denotes summation excluding the term 


r 
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r = 1. Having the solution for the homogeneoys 
equation, one uses the method of the variation oj 


parameters to obtain the general solution, 


G\(4F? + F’) 
GiGs' — Gi'Ge 
G(4F? + F’) 


——- 6) 4+ AG + BG. ad 
Es! - GC’. ; 2 


dd — 


An almost completely similar analysis may be mack 
for the case of converging flows after a standard trans 
formation from dn to sn functions and a change oj 
modulus have been made in Eq. (2.33) so that the as 
sociated energy equation has been reduced to an ordi- 
nary Lamé equation. 

The case when [—2 + m(3a)’*]/2 is integral may 
also be treated for certain specific Reynolds Numbers." 

Although a general analytical solution has been given 
for certain specific parametric values of the Reynolds 
and Prandtl Numbers, the result is so complicated and 
involves so much additional computation in the solution 
of the auxiliary equations that it appears advisable to 
use numerical methods ab initio. 

One of the more desirable procedures appears to bi 
the method of finite differences. In using this method, 
Eq. (3.5) is rewritten in the abbreviated form, 


—y" + P(d)y = Q(d) 1.6 
One then approximates this equation with the follow 
ing finite difference scheme :'' 
run - 2F + Tia — 
h? 
12 


(Y,..” + 10Y,” + Y..,") =0 E 


where 
V ee . 
"= 


PY, — Q, (4.8 


h is the interval width, and 7 is a free index varying over 














Fic. 8. The thermal distribution in a converging channel ol 
total opening of 10° for a flow with R = 5,000 for parametrit 
values of the Prandtl Number. Only the normal solution 1S 


shown. 
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THERMAL DISTRIBUTIONS 


the range 1,2,...,” —1,anda subscript denotes the 
value of the indicated functions at ® = ai/2n for 
symmetrical profiles (obvious adjustments can be made 
in the case of asymmetrical profiles). This difference 
scheme is not the usual one, but it is a more refined 
method where the first neglected term in the Taylor 
expansion has a magnitude h*Y,"'/20. The substitu- 


tion of Eq. (4.8) into Eq. (4.7) gives 


h? } dh? y 
ra) sais 12 Pas| Po - v, [2 5 6 P| = ia Y; 1 Xx 
= ) 


h? h? 
[ -" P| =~ 75 lO + 100 + Qual 49) 
[ 2 2 
fori = 1,...,m — 1, and the boundary conditions 


3.6) become 

(4.10) 
Since each ordinate is a function of the ordinates ad- 
acent to it, the successive numerical evaluations of 
Bq. (4.9) lead to the following set: 


Vint = Gini + Oui; (4.11) 
for 1 = 1,2, .,n — 1, and the imposition of Eq. 
4.10) on the member of set (4.11) for which 7 = n — 1 
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The thermal distributions in a diverging channel of 
1,342 for parametric 


Fic. 10. 
total opening of 10° for a flow with R = 
values of the Prandtl Number. 


leads to the complete determination of the Y,; and the 
numerical solution of the energy equation. 

Eq. (3.5) has been solved by the outlined method of 
finite differences for five Prandtl Numbers (0.5, 1, 2, 5, 
and 10) and for the normally symmetrical profiles for 
three Reynolds Numbers (684, 1,342, and 5,000). The 
numerical results are given in Figs. 6-11. 

The thermal profiles for the normal solutions of con- 
verging flows in Figs. 6-8 have similiar properties——a 
definite thermal boundary-layer behavior with rapid in 
creases to a maximum temperature occurring near the 
walls and an almost constant distribution across the 
remainder of the channel with a minimum temperature 
along the channel axis. 

The characteristic behavior in diverging channels, 
Figs. 9-11, is quite different. The thermal distribu 
tions contain oscillations with maxima and minima 
above and below the wall temperature and with their 
number depending on the Prandtl Number of the fluid 
and the number of fundamental flow regions in the 
hydrodynamic profiles. For certain critical Prandtl 
Numbers, there exist points within the channel where 
the temperature becomes infinite with the sign of the 
temperature changing upon passing through the critical 
Prandtl Numbers. The variation of the profiles with 
Prandtl Number is similiar to the stability phenomena 


in solutions of the Mathieu equation. The reader 
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Fic. 11. The thermal distribution in a diverging channel of 
total opening of 10° for a flow with R = 5,000 for parametric 
values of the Prandtl Number. 


should note the incompatibility of the assumed incom- 
pressibility of the fluid and the extreme values of tem- 
perature in the profiles. 
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In view of these theoretical conclusions, a compari- 
son between the theoretical predictions and experi- 
mental data obtained from interferograms of the real 
behavior of the laminar flow in channels with a constant 
wall temperature should prove to be highly interest. 
ing. 
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uadruple Velocity Correlations and Pressure 


Fluctuations in Isotropic Turbulence 


1 


ss | 


MAHINDER S. UBEROI? 
The Johns Hopkins University 


SUMMARY 


It is shown that the correlation of fluctuating static pressure 
in an incompressible and homogeneous turbulence) with any 
fluctuating quantity in the flow field can be expressed in terms 
of the correlation of the same quantity with two or more com 
ponents of the velocity. 

The correlations of pressure with itself are investigated in de 
tail for the case of isotropic turbulence. These correlations can 
be expressed in terms of correlations involving two velocity com- 
ponents at a point and two velocity components at another point 
A postulated relation between the fourth-order and second-order 
correlations is investigated. This relation is satisfied, for ex 
ample, if the joint probability density of the four components of 
velocity is Gaussian. The consequences of this relation are com 
pared with the measurements of the fourth-order correiations 

The root mean square (r.m.s.) pressure and pressure gradients 
are computed from second-order correlation for a range of turbu 
lence Reynolds Numbers. Since the pressure gradient is related 
to diffusion of marked particles from a source, the computed pres 
sure gradient level is compared with that calculated from a set of 


diffusion measurements. 


INTRODUCTION 


ren in isotropic turbulence the dynamical 
(correlation) equation does not contain any pres- 
sure correlations, the correlation of pressure with itself 
is of interest in turbulent diffusion and in the study of 
sound generation by turbulence (for very low Mach 
Numbers of turbulence velocities when there is little 
effect of sound on the turbulent field—i-e., only a slight 
amount of energy is drained away from the turbulent 
energy in the form of sound waves). The static pres- 
sure fluctuations are also intimately connected with the 
onset of cavitation in turbulent flow of liquids. 

The pressure gradient fluctuation, in connection with 
diffusion in isotropic turbulence, was first considered 
by Taylor.' He postulated that (1/p?) (Op Ox,)? is 
of the order U7(0L/0x,)*, independent of Reynolds 
Number of turbulence. Heisenberg? later derived, 
from detailed spectral considerations, an expression for 


(Op Ox)”, 
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l Op 2 ae ol 1 . 
p Ox; Ox; / 


where Ry, = l?/vV (OU/0x,)? is the Reynolds Num- 
ber of turbulence and the symbols have the usual mean- 
ing. 

Obukhov* derived an expression for the pressure cor- 
relation for distances between two points whose spac 
ing is within the Kolmogoroff’s “inertial subrange."’ 
Batchelor’ considered the pressure correlation in iso- 
tropic turbulence. These investigators have, from 
the very beginning, assumed a simple relation between 
the fourth-order and the second-order correlations. 
In order to compare the primary result with experi- 
ment, it is necessary to have the results in terms of the 
fourth-order correlation. 

Before taking up pressure correlation, we shall con- 
sider correlations involving static pressure in general 
and show that these correlations can be expressed in 
terms of correlations involving more than two velocity 
components. This is convenient from the theoretical 
point of view, since it may then be possible to express 
higher order correlations in terms of second-order cor 
relations on the basis of some plausible hypothesis, 
such as those used by other investigators. It may be 
also noted that Heisenberg’s* expression for the spec 
tral transfer term is a hypothesis for the triple corre 
lation in terms of the second-order correlations. The 
fact that correlations involving static pressure can be 
expressed in terms of correlations involving more 
than two components of velocity is also convenient 
from an experimental point of view, since at present 
there exists no technique for the measurement of static 
pressure fluctuations. On the other hand, the stand 
ard hot-wire technique can be extended to correlations 


involving more than two velocity components. 


GENERAL RELATIONS 


Consider a statistically homogeneous and incom 
pressible turbulence with no mean velocity (i.e., a 
‘“‘box turbulence’’). The equations of motion for an 


arbitrary point X are: momentum, 


OU; OU; 1oO ; 
U; =—- P vAU (1) 
rey) Ox; p Ox; 
continuity, 
dU, /d«, = 0 (la) 





where the symbols have the usual meaning and re- 
peated index means summation over the index. Tak- 
ing the divergence of Eq. (1) and making use of Eq. 
(la), we get the well-known relation, 
r r 2 / 
Ou; ae. 
“~ (2) ie = 
Ox OX; 


I 1 { ®UUq'(y) dV 
pq'(é') = / ay) ony) (4 
p Ar 


Oy OY; y—é 


where dV(y) isa volume element. Let 
ou,’ U; 


1 op | . 
Ox;,/Ox,' 


p Ox Ox; p Ox,;’Ox,’ 
Substituting this into Eq. (3), interchanging the order 


The fact that the pressure in an incompressible viscous 
of differentiation, and averaging, we get 


fluid satisfies Poisson’s equation means that it is not a 
yrimary variable, since by using the well-known solu- Sess 
| es , 6 1 — d*pp’ O1U,U,U,' Uy 


tion for Poisson’s equation we can express the pressure = (5 


in terms of the velocity derivatives. p* OF 0§,0E,0E, OF ,0E EOE 
Multiply Eq. (2) with q’, any quantity at another For the isotropic case 
point x’, and taking a statistical average, we get a4) 1 d4[r()] 
1, &p eu U, OEE DE, 7 Or! 
p q O01 ;OX; °_ Ox ,Ox; 


and the bilaplacian equation, Eq. (5), simplifies to 
Making use of the fact that the turbulence is statisti- 
cally homogeneous and interchanging the order of dif- 


L O'rpp'(r) — 8U UU, Uy’ 
1, i - 


ferentiation and averaging, we have p’r or‘ OE ,0E,0F,.0E; 


1 O°pq’(E) 0?U;,U;q' (&) After four repeated integrations, for any fixed /, we get 


(3) 
p OF OE; OE OE; ° ror? . 
sstit : i yg QUiUiU;' UY 
where x" = x + € and p stile 6r Jr mee Oy,OV,OV,OY; ? 
p(x)q'(x’) = pq'(&) (7) 


An overbar denotes statistical average. Eq. (3) is just) U/,U,U,’U;/ is a fourth-order isotropic tensor and its 


the Poisson equation for pg’, and its solution can be form can be derived using some simple results of the 


written immediately as invariant theory.® Thus, 


ee Rt) — &) ~ 2") - 4)... 
l il il k l i { €€)ExE1 a 
r 


I nn ee 7 
Ran (1) 6 ij8x1 + 9 [Run™ (Y) — Reg tr) [64552 + 0; 10 jy + 


l l 
9 Roa (#) a Reg) [EE 6x1 + ££ 6;; | + yy? (2R,,:°"(r) aaa Ran") = Rar) [é E,6, + 
2r? 


Ser 
> 
Str 
as 
> 
x 


wre 


F911 = £5 


The various fourth-order correlations are shown in Fig. 
l. 

Some investigators*® 4 have assumed that joint prob- 
ability density of U, U;, U,’, and U,’ is Gaussian; 
under this assumption the quadruple correlations are 
related to the double correlations by the following: 


U,U,U,'U! = U,U,': UU + UU UU,y! + 
U,U;-U,'U,' (9) 


and, in particular, 


Rull(r) = 2RVr) + Rit 
Ko (7) = 2K, (r) + £0" 

Ror). = RY @)R,"(r) (10) 
Kr) = Rik 

Rae) |, a 


These relations and the expression for U;U;U,'U/' for 
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the case k = / = i were first given by Millionshtchikov.* 
The second-order correlations R;' and R,,” [essentially 
the von Karman-Howarth g(r) and f(r) | are connected 
by the continuity equation of incompressible flow.’ 


R,"(r) = R,'(r) + ih R,|(r) (11 
2 Or 
The hypothesis of Gaussian joint probability density 
makes all odd moments (e.g., triple correlations) zero. 
However, the above relations are only integral condi- 
tions on the joint probability density, and the hypothe- 
sis of Gaussian probability density is not strictly neces- 
sary for their validity. Without making any such 
hypothesis about the probability density, we simply 
take Eq. (9) as a plausible hypothesis. Substituting 
Eq. (8) in Eq. (7) and carrying out some algebraic 
simplifications and integrations by parts, we get 
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order 








2 sie ! ‘ig dy 

| pp'(r) = Ri'(r) a Ri? oo if | 2Ras + Rin” = Sa ta + “ea + R,,"' — ky = UR") = 

Ms ? i y 
(12) 


The pressure gradient is of special interest since it is related to the diffusion of marked fluid particles from a fixed 
surce (see next section on “Pressure Gradient and Diffusion Measurements’). Now, 


(Op/O.x1)? = lim [Op(x,)/Ox;]- [Op(x1’) /Ox,"] 
r—0O 


where x;' = x; + r. Interchanging the order of integration and averaging, we have 
(Op/O0x,)? = —(O*pp’(r)/Or?), <¢ 


Substituting for pp’ from Eq. (12), we get 


1 /Op \? d* ‘ 1\ 
: P) cain | ~ 2Ru'(r) — Raxt*(r) 0 + sf (Ran™® + Rul! — 2Rya"! — 4Ryi") = (13) 
p” \Ox, dr? 0 y3 





In order for the above integrals to converge, it is neces- 


ary that the series expansion of 


n 
(Ren + R,,"* — 2R,,'' — 4Rni") m 
at the origin begin with terms of order r*; therefore, f 


7° 4 
, ii.” — Rn” — 2. = Ra") | = 0 
r=0 


dr? 
We have made use of this fact in deriving the expres- tf 
ee ee —= = Retr) —_ 
son (13). Under the assumption of Eq. (9), Eqs. nbcsaiilay —* 





12) and (13) become 


] a ¥* d . 
- pp '‘(r) = 2 / (v — )| R,4( y | dy 
p- af y/tdy 


(12a) iT nn IT 

l : ( 2 dy R (r) 

(Op/Ox,)? = 4 / | ; Ri) | I) (13a) an 
p J0 dy y 


For vanishing turbulence Reynolds Number, 


R, > O(R, = V U)? X/v) nn 
Roy (r) IT 


\? = U;? (ol ! Ox)" 





1 


where 





R,'(r) can be calculated neglecting the inertia terms of 


the equations of motion,* iy 
| n 
R,'(r) = U,? exp [—r?/2a?] Ry, (r) @ 


Substituting this result in Eqs. (12a) and (13a) we get 





pp'(r) pl? exp [—r?/A?7] = p*(R,')? (14) 


(Op/Ox,)? = 2p7l 2 * (15) IT mm 
On the basis of Kolmogoroff’s hypothesis of local isot a (r) ff 


topy, for large Reynolds Numbers and r with the 








inertial subrange’’ we have® Le r ~ 
y , ’ ’ 1 r . oa 
( l l )* = 2[l 2 = R,*( Yr) ( i( ve) Fic. 1 he fourth order correlations 
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(16) 
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where L’ is the velocity along the displacement vector, 


z= x’ — x’, 
x ls ( ee) 
-- _ ists = 
= : 2 dt 


is the energy dissipation per unit mass of the fluid, the 
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~ 3/4 Vf se 
length ¢ = v’‘e’*, and C, is an absolute constant, 
e and \ are connected by Taylor's decay equa. 


tion,! 
(1/2) (dU,U;,/dt) = —15v U,?/ 


Hence, for values of r within the inertial subrange, 


l ms 1] 2 « 1 2 1 ’ 
(Pp? — pp’) = / | RX) ydy+r | R,'(v | = (17 
2p” Jo Ldy ; * dy : y ' 


The integrals in Eq. (17) can be evaluated by using 


Eq. (16). The result is 


l - : C; a * 
are (‘) 
o 


or 

(1/p)? (p — p’)? = (U — U')® (18) 
This relation was first given by Obukhov,’* and Eqs. 
(12a) and (18a) were given by Batchelor.‘ 
PRESSURE GRADIENT AND DIFFUSION MEASUREMENTS 


In reference 10, the ‘‘Lagrangian’’ microscale \, was 
related to the acceleration in isotropic turbulence with- 


out the neglect of viscosity. A brief résumé follows: 


2 Ley (zy) 
2 V'*L\de dt 


is the r.m.s. velocity in the y-direction. 


(19) 


where |’ 
Squaring and averaging the v-component equation of 
motion, we get 


(>) l (oY) + AV)? 
= - )- 
dt p” \Oy : 


The term (Op/oy AV) = 
terms in the expansion of (AV)? can be expressed in 


(20) 
0 because of isotropy. All 


terms of the derivatives of R,'(r) by making use of the 


relations given by von Karman and Howarth.’ Thus, 


AV) 33 ee 
}* = 
fs or? oS a 


Townsend!! has made measurements of 


[O'R,'(r) or*], 0 


his result is 


(AV) 33 Bea 
2 = = 
2 OF 4h a® 


33 V2? [30 5 P 
2 he 7 + 0.2R, (21) 


For the isotropic turbulence 


dV'?/dt = — 10v (V"?/X?*) 
(the ‘‘decay equation,’’ due to Taylor'). Thus, 
(A2/V’ 4) (dV'/dt)? = 25/R, (22 


Substituting Eqs. (20), (21), and (22) in Eq. (19), we 


get* 
p?V'4 E 20.7 =] ! 
= _ — (23) 
\*(Op/dy)? A” RP RR 


Heisenberg” has calculated an expression for pres- 
sure gradient using the hypothesis that various Fourier 
components of the velocity field are statistically inde 
pendent. As far as the fourth-order correlations are 
concerned, the consequences of this hypothesis can be 
shown to be equivalent to Eq. (9). His result is 
2p7/4 

at = 0.038«R, 

A*(Op/ Oy)? 


(24) 


Instead of Heisenberg’s original extrapolation formula 
for the stationary velocity spectrum, Chandrasekhar’s 
solution'* was used to calculate Eq. (24). x is a “uni 
versal’ constant that must be determined from experi 
mental data. Heisenberg found from turbulence de- 
cay that x = 0.85. Lee" gives its value as 0.13, based 
on skewness of the probability density of (Ol’; Ox) as 

* This is the first step toward the establishment of a possible 
If we 


consider nondecaying isotropic turbulence, then the Lagrangian 


relation between Lagrangian and Eulerian correlations 
correlation 


: - aes "T).\2 
UNUM + 7) = R(x) = UF + D(-1y'(4 ) d 


n=0 dt (2n 


We can express each term of this series in terms of the Eulerian 
correlations by using equations of motion. The first nonconstant 
term involves pressure or fourth-order velocity correlations. The 
next term will involve still higher order correlations. It may be 
possible to relate higher order correlations with second-order 
correlations by making a further hypothesis about the joint prob 
at two points. We have just 


The next step will 


ability density of the velocities 
carried out the first step of this computation. 
involve many more calculations. This is probably not the best 
method of attack of the problem of relating the Lagrangian and 
Frenkiel'? has related Lagrangian 


the Eulerian correlations. 


and Eulerian microscales under some assumptions 
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Proudman" has calculated 


weasured by Townsend."! 
double and triple correlations from Chandrasekhar’s 


lution of Heisenberg’s self-preserving spectrum. He 
has estimated its value as 0.45 from comparison of theo 
retical correlations with measured correlation. Heisen- 
perg’s results are only approximate, and the value of x 
depends on the range of spectrum (or correlation) 
which is made to fit the experimental data. We have 


taken k = 0.45 as a compromise value. 


EXPERIMENTAL OBSERVATIONS 


The measurement of various correlations was made 
in isotropic turbulence at R, = 60. The correlations 


Ru(r), Ran""(1), Ran” (r), and two independent com- 
binations 


(Rk, +- | ae + i). i tR.,**) 
(Rit! + Ran” + 2Rent + 4Rpi™) 


were measured. The measurements are compared 


with the hypothesis of Eq. (9) in Figs. 2, 3, 4, 5, and 6. 
While measuring the fourth-order correlation, the cor- 
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responding second-order correlation was also measured 
simultaneously. For sufficiently large distance, the 
correlation between l’;? and l’,;? disappears and l'?U’ 
tends to a constant value L’,*- l’;*, which is equal to U2 
for isotropic turbulence. If we assume that the prob- 
ability density of U, is Gaussian, then U?° l\4 = 
1 3. This is indicated by a dashed line in Figs. 2, 3, 
and 6. 
hypothesis of Eq. (9) and is fairly well satisfied in prac- 


This assumption is less restrictive than the 
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tice. The correlation R,,"" (Fig. 7) was also meas- 


ured; this is not independent of the above five corre- 
lations. 
stant value, it is of interest to see whether the meas- 
Within the 
experimental accuracy, all the fourth-order correlations 


Since the hypothesis of Eq. (9) gives it a con- 
ured value deviates from this constant. 


are close to those computed from corresponding second- 
order correlations using Eq. (9). 

Batchelor‘ has presented R. W. Stewart’s measure- 
ment of 


in support of the hypothesis of Gaussian joint prob- 
ability density of the four velocity components enter- 
ing in Eq. (5). 


if this hypothesis is satisfied. For the purpose of re- 
lating pressure correlation with second-order velocity 


correlation, it is only necessary to assume the hypothe- 


1 2 > tr? 
; pp'(r) = Ri(r) _ U,? aa 2? / 3R, + R, ~ LR, —_ 2 (R, + Rs oe 2R3) 
p- J) _ 


R», was measured for both positive and negative values 
of r (see Fig. 3) and the point r = 0 was determined 
from the fact that it is an even function of 7-R; and R; 
were measured for positive values of 7, and there is 
some uncertainty (= 0.025 in.) in the determination of 
the point y = O for these two correlations. This un- 
certainty, the experimental scatter, and the fact that 
small differences between relatively large quantities 
are required to compute the pressure correlation have 
made the results extremely uncertain. 


Two widely different sets of curves were drawn for 


each measured correlation, taking into account the ex- 
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sis of Eq. (9) for the fourth-order correlations involy- 
ing two velocity components at one point and two ye. 
locity components at another point. The quantity 

) 

(U, — U,')4/(U, — U,')? 

includes fourth-order correlations involving three ye. 
locity components at one point and one velocity com. 
ponent at another point, and these correlations do not 
enter in the pressure correlation Eq. (5). In passing. 
we mention that, if the fourth-order correlations enter. 
ing in Eq. (5) are related to second-order correlation 
by Eq. (9) and if these new fourth-order correlations 
are related to the second-order correlations by the rela- 
tion 


U,U,U.Uy = UU: UU! + UU Uj, + 
UU; U,l 


then 


This is similar to the hypothesis of Eq. (9) and is less 
restrictive than the hypothesis of Gaussian joint dis- 
tribution of the four velocity components. The fourth 
order correlations U;U,;U,.U,' have been measured. 
The results and further discussion on this point are 
presented in reference 5. 

Let R;, Ro, and R; denote the correlations R,,"", R,, 


and the combination 


(R,"* + a + eg + Lz.,"* 


respectively. This set of correlations, among other 
possible sets, suffices to determine the pressure corre 
lation. In terms of the above notation, the Eq. (12 


becomes 


dy 


3 F 


perimental scatter and the uncertainty in the deter- 
mination of the point r = 0. Two pressure correla- 
tions were computed from these two different sets of 
The result of these com- 


putations and those from the second-order correlation 


fourth-order correlations. 


are compared in Table | and Fig. 8. 

It appears that for the computation of pressure cor- 
relation one should measure directly the differences 
(3R, + R; — 4R;) and (R, + R, — 2R;) instead of the 
individual This the difficult 


problem of measuring small correlation between rela- 


correlations. involves 


tively large quantities. 
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In view of the fact that hypothesis of Eq. (9) is ap 
proximately satisfied we may use double-correlation to 
compute r.m.s. pressure and pressure gradient. How- 
ever, we must not overlook one unfortunate circum- 
stance. It is possible that hypothesis of Eq. (9) is 
approximately satisfied but the r.m.s. pressure and pres 
ure gradient computed from using this hypothesis are 
sill in error, because the r.m.s. pressure and pressure 
gradient depend on the difference of various fourth-order 
correlations. 

Eqs. (12a) and (13a) express the desired quantities 
in terms of the longitudinal second-order correlation 
R,'R;' and R,,” are connected by Eqs. (11). We need 


toexpress (0/Or)R,' in terms of R,.”. Since 
O l _ Oo , 
x '= A ye R, 
or Zr? Or Or 
therefore, 
Oo 2 


. *, ra) 
3 / R,,"(y)y? dy 
9 Oy 


r3 


using this result, Eqs. (12a) and (13a) become 


p° : ~ a S vi weg 

( *) U;* = ° / | 
p- U,?* Jo 0 dy 
(*) 2 72 16 “ a R.% 2d 2dr 
2¢ = rye ’ = 
dx; p 1 J / | rs y)yr ay r 
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a > dr 
R,,"(y)y? dy 


Root mean square pressure and pressure gradient were 
computed from measured R,,"(r) for various Reynolds 
The results of the computation are shown 
The data for the two highest Reyn- 
taken 16. The 


Numbers 
in Figs. 9 and 10. 
reference 


olds Numbers were from 








data in reference 16 are given in terms of grid Reyn- 
lds Numbers. Turbulence Reynolds Numbers, R, 
and \/.\J (where M/ is the mesh spacing) were estimated 
irom grid Reynolds Numbers by assuming that l*, U° 
~f-'in the initial period. lp is the mean velocity. 
Eq. (23) relates the ‘“‘Lagrangian micro-scale,”’ \,, 
with the r.m.s. pressure gradients. Measurements of 
\, have been made by Simmons" and Collis,’ and an 
extensive set of measurements is given in reference 10. 
The experimental data (Table 2) are used to compute 
the pressure gradient [Eq. (23)], and the results are 
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TABLE | 
Op \? / , ae Pom | (OP \7 1+? 

9 9 r 92 ( p / p~ l “wae 2p? / Pp 

p?/p? Ui? |\ 0x / / / Ox] _ 
Two extreme 3.3 85 0.28 
values com- 
puted from 
fourth-order 
correlation 0.45 11 0.20 

0.50 12 0.20 
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TABLE 2 
Op\2 | p* I 
9\2 91) 7 ) dr” ( : ) | , 
r, hin, 2X 20.7 a= l oy Vv. (2 2 
Ry In. In An? Rx? Rx p*V'4 ‘ ed 
Simmons, 1935 37.0 0.16 0.08 8.00 0.072 7.930 0.356 
Collis, 1948 23.0 0.14 0.20 1.00 0.130 0.870 1.070 
Collis, 1948 38.0 0.14 0.20 1.00 0.070 0.930 1.030 
Reference 10, 1951 29.0 0.17 0.11 1.78 0.097 1. 680 0.463 
Reference 10, 1951 32.0 0.12 0.12 2.00 0.086 1.914 0.725 
Reference 10, 1951 36.0 0.41 0.36 2.60 0.074 2.530 0.629 
Reference 10, 1951 42..5 0.23 0.235 1.92 0.061 1.860 0.735 
Reference 10, 1951 $3.5 0.165 0.245 0.910 0.059 0.851 1.084 
Reference 10, 1951 49.0 0.35 0.22 5.04 0.052 $1999) 0.448 
Reference 10, 1951 61.0 0). 23 0.3 0.914 0.040 0.874 1.070 
Reference 10, 1951 67.0 0.34 0.22 1.76 0.036 41.720 0.461 
Reference 10, 1951 74.0 0.19 0.27 0.990 0.0382 0.958 1.022 
a a * Loitsiansky, L. G., Seme Basic Laws of Isotropic Turbulence. 
© Reference 10 C.A.H.T., Moscow Report, No. 440, Translated as N.A.CA 
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SUMMARY 


A field mapping method that is an extension of the usual steady 
state transfer function methods of analyzing closed-loop systems 
for such a system, the open-loop 
as a potential flow field in 


ispresented. It is shown that, 
transfer function may be represented 
the complex plane (either the s-plane or, better, log s-plane) in 
which sources and sinks are placed at the locations of the poles 
ind zeros of this open-loop transfer function. The values of the 


open-loop transfer function along the imaginary axis of this 
plane represent the steady-state (frequency) response of the open 
loop, and from their knowledge the flow field in the vicinity of 
closed-loop roots may be determined graphically by means of 
walytic continuation. The closed-loop roots occur at those 
points s for which the open-loop transfer function achieves the 
value —1. The method allows a quick determination of the 
closed-loop roots for any value of loop gain and further shows 


the dependency of the roots on additionally introduced phase 


shifts 
List OF SYMBOLS 
Laplace transform operator 
6 = direction of fuselage centerline in space 
6 = deflection angle of control surface 
# = complex potential function (¢ + jy) 
¢ = potential function 
¥ = stream function 
E = angle of radial line in s-plane 
db. = decibels, measure of amplitude ratio 
r = magnitude of radial line in s-plane; undamped fre 
quency 
® = angular frequency, rad. per sec.; w = rsin E 


as a proportion of critical 


It 


damping factor, expressed 
damping 


INTRODUCTION 


_ CLOSED-LOOP DYNAMICAL SYSTEM is one in which 

the output of the system is fed back into the sys- 
tem as an additional input. Any closed-loop system 
may be pictured schematically as in Fig. 1. If the ele- 
ments represented by the symbols yu and 6 in Fig. | can 
be represented by a set of linear differential equations 
with constant coefficients, then the analysis of these 
elements by the method of the Laplace transformation 
yields general rational algebraic fractions as the repre- 
sentation terms of the Laplace 


of these elements in 
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transform operator, s.' These functions are termed 


“transfer functions,’’ and their multiplication by any 
general input (in operational form) yields the response 
to that input (also in operational form) 

The transfer function for the closed loop may be seen 


from Fig. | to be given by Eq. (1 


0,6; = w/(1 + wB) (1) 


The poles of 4/9; (i.e., the roots of the equation ob- 
tained by setting the denominator equal to zero) give the 
decay time constants and frequencies of the various 
modes of the system. The poles of %/6; from Eq. (1) 
are those values of s for which y@ is — 1. 
are general rational algebraic fractions, it follows that 
values of s either will be real or will occur in com- 
For the system to be stable, 


Since uw and 8 


these 
plex conjugate pairs. 
these values of s must have negative real parts; a 
pair of pure imaginary roots results in a neutrally stable 
system. 

The complex value of a transfer function for an imag- 
inary value of s is the ratio (amplitude and phase) of 
the output to the input of the corresponding element 
for a steady-state sinusoidal input. A compilation of 
the ratios for various input frequencies is termed the 
‘frequency response’ of the element. In the usual 
methods of frequency analysis, the degree of stability 
of the closed-loop system is inferred from the margins 
with which the complex values of the open loop 8 miss 
being —1 for the imaginary values of s, since this cor- 
responds to neutral stability. In the 
method, presented herein, the complex value of s for 
which yw is actually —1 is obtained from the frequency 
There 


field-mapping 


response of u@ by means of analytic continuation. 
have recently been published a number of papers on this 
subject,*'” and a summary of these and other recent 
developments in the field of automatic control may be 
found in reference 11. 


THE POTENTIAL ANALOGY 


The field-mapping method is based on the consider- 
ation of the transfer function representing any linear 


element as a two-dimensional potential field. This 
relationship is known as the potential analogy,** and 
the manner in which it comes about may be seen by 


considering the general rational algebraic fraction F(s), 
which arises in the solution of the equations represent- 
ing the dynamics of this element by operational means, 
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If the roots of P(s) = 0 are —a; and those of QO(s) = 0 
are —b,, then the fraction may be written (where a 
and 0; are, in general, complex numbers) 

p 

II (s + a;) 

F(s) = K = 


II (s + b;) 


1 
It follows that 


log [F(s)] = log K + > log (s + a,) — 
1 


1 


q 
> log (s + b;) (4) 
t=1 


Consider now the complex potential function ®(s) 
representing radial source flow at the points —c; (where 
c; is, in general, a complex number).* 


¢g(s) + 7W(s) = 


Comparing Eqs. (4) and (5), it is seen that log [F(s) | 
may be considered to be the complex potential func- 


@(s) = —C log (s + ¢;) (5) 


tion for a flow field in which the poles of F(s) act as 
sources and the zeros of F(s) act as sinks (1.e., negative 
sources). Since F(s) may be considered to have an am- 
plitude A and phase B, then 
. 3( 
F(s) = A(s)e” 
and so 


(6) 


log [F(s)] = log[A(s)] + jB(s) 
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Comparing Eqs. (5) and (6), it is seen that the po 
tential function ¢ corresponds to the logarithmic ampli- 
tude of the transfer function /(s) and that the stream 
function y corresponds to the phase of F(s). 


METHOD OF MAPPING 


The (logarithmic) amplitude and phase functions oj 
F(s) were shown to be analogous to the potential and 
stream functions of a flow field determined by sources 
and sinks placed at the locations of the poles and zeros 
of F(s). If the values of the complex potential fune- 
tion of such a field are known along a line of finite length 
within this field, then the flow field may be uniquely de 
termined except in small regions enclosing the singular 
ties within the field.’ This process is known as analytic 
continuation. In its application to the field-mapping 
method, it shall be assumed that the frequency r 
sponse of F(s) is known |i.e., the phase and amplituc 
of F(s) are known along the imaginary axis of the 
plane], and the flow field shall be constructed by flux 
plotting,’ which is the graphical equivalent to analytic 
continuation. 

In order to increase the convenience in application 
of the method in the frequency range of interest, the 
field shall be developed in the plane log s.’ 
contplex variable, expressible in the vector form 


Since s is a 


re” it is seen that 
logs =logr+j jE i 


The log s-plane is thus one in which the abscissa 1s the 
logarithm of the length of the radial line drawn from th 
origin to any point in the s-plane, and the ordinate 1s 
the counterclockwise angle between this line and the 
positive real axis. Since this angle repeats itself ever 
2a rad. it is seen that the log s-plane will consist of a 
series of horizontal strips 27 high and that the flow field 
will be identical in each strip. Thus only one such strip 
need be considered. In addition, it may be noted that 
all singularities either lie along the real axis or occur 11 
conjugate pairs. Thus the flow field in the s-plane 1s 
symmetrical about the real axis, and in the log s-plane 
The 


positive real axis transforms into the base of this strip, 


only a strip between 0 and 7 need be considered. 


the negative real axis transforms into the top of the 
strip, and radial lines from the origin transform into 
horizontal lines in the strip, the vertical coordinate / 
being the angle that the radial line makes with the real 
axis, as shown in Fig. 2. 

In the application of the method to closed-loop sys- 
tems, the information available as a result of the usual 
frequency analysis is assumed to be known. This 
means that the amplitude and phase of F(jw) are known 
[where F(jw) represents the steady-state open-loop 
transfer function], and, hence, the potential and stream 
functions of the flow field are known along the imag! 
nary s-axis (i.e., along the line E = 90° in the log : 


plane). From this information the flow field in the log 
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Fic. 3. Block diagram for numerical example. 
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Fic. 4. Representation of open-loop transfer function in the 
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Fic. 5. Frequency response of open loop 


stip can be constructed by making use of the fact 
that the potential and stream functions are orthogonal 
functions, provided the same scale is used for ordinate 
and abscissa. From Eq. (7) it may be seen that the 
horizontal length of one decade must be log 10, which is 
2.303. The total height of the strip is 7, or 3.1416. 
Three-cycle semilogarithmic graph paper of these di- 
mensions may be obtained from the Raytheon Manu- 
facturing Company, Oxnard, Calif. The size of the 
mesh is determined from similar considerations, and it 


may be shown that 20° of phase corresponds in size to 
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3.03 db. The procedure followed in constructing the 
fieid is as follows. From the known amplitude and 
phase curves of the steady-state open-loop transfer 
function, spacings of, e.g., 3 db. in amplitude and 20° 
of phase are noted along the line E = 90° in the log s- 
plane. Through these points are sketched lines with 
such a slope that the lines through the amplitude and 
phase curves form curvilinear squares. The resulting 
mesh is uniquely determined from the values along the 
line & = 90°, and the sketching may be done by eye with 
reasonable accuracy. The size of the mesh is arbitrary 
so long as the relative lengths along the potential and 
stream lines are maintained; a mesh size of 10° versus 
1.5 db. usually provides sufficient accuracy combined 
with ease of plotting. 


NUMERICAL EXAMPLE 


As an example of the use of the method, the case of a 
small high-speed aircraft stabilized to fly a level course 
by means of a free gyro shall be considered. The block 
diagram for such a system is shown in Fig. 3. The 
aerodynamics of the aircraft may be represented by a 
transfer function of the following type: 


A 10 (s + 10) 
- . (8) 
6 s(s*? + 4s + 104) 


The servo transfer function is assumed to be given by 
Eq. (9). 


[s] = K,/(1 + 0.02s) (9) 


The gyro is assumed to have unity response. The 
open-loop transfer function may thus be obtained from 
Eqs. (8) and (9) as 

10K,(s + 10) 


s(s? + 4s + 104) (s + 50) 


This function may be represented in the s-plane by a 
zero at s = —10and poles at s = 0, s = —50,ands = 
—2 + j10, as shown in Fig. 4. To obtain the fre- 
quency response of the open loop, the substitution 
s = jw is made in Eq. (10). The resulting amplitude 
and phase curves may be seen in Fig. 5 for an assumed 
value of K, = 89. Fig. 6 shows the complete flow field 
in the log s-plane, constructed as discussed in the pre- 
ceding section from the known amplitude and phase 
values along the imaginary s = axis (determined from 
Fig. 5) and from the knowledge of symmetry along the 
top and bottom of the strip. In this particular ex 
ample, the locations of the singularities are known, 
and thus additional aid in drawing the field may be 
had. However, such information is not necessary ex- 
cept to determine the field in a small area around the 
singularity. From Fig. 6, the following characteristics 
of the system may be easily determined: 

(1) The closed-loop roots—namely, the locations of 
the intersections of the 180° phase stream lines and the 
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unity gain potential lines. These in turn determine 
(a) the damping factors of the oscillatory modes [if the 
vertical coordinate on which a root lies has the value 
Fj, then ¢ = sin (F; — .90°)]; (b) the damped natural 
frequency w of the oscillation, from w = r sin /;. 

(2) The effects on the locations of the roots of small 
changes in gain or phase shift of the open loop. 

It should be noted that it is not necessary to develop 
the entire flow field to determine roots corresponding 
to less than critically damped modes. Rather, just 
the region in the vicinity of the 180° phase lines and 
0-db. amplitude lines need be developed. For poorly 
damped roots, this was seen to be near the imaginary 
axis, and thus the location of the root may be deter- 
mined quickly with good accuracy. For well-damped 
roots, the root lies near the real axis, and it becomes 
more difficult to determine its location accurately. 
However, in such cases great accuracy is not usually 
required. For the problem investigated above, only 
the region shown in Fig. 7 would have had to be de- 
veloped, a job requiring only a few minutes. The con- 
struction of the field may be made with any desired 
accuracy, the only restrictions being the effort put in by 
the operator and the accuracy with which the basic data 
are known. 

Another advantage of the method is that the effects 
of shaping networks in improving the stability of the 
system may be determined by the direct superposition 
of the flow field of the shaping network on the flow field 
of the open loop, since potential fields are directly addi- 
tive. In particular, the fields representing various 
combinations of one and two singularities may be drawn 
for standard spacings of the singularities, thus enabling 
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Fic. 7. Determination of oscillatory mode 

the operator to know at a glance the effects of the 
standard networks, such as low-pass filters, bridged-T 
filters, etc. The logarithmic abscissa is particularly 
useful in this respect, since such fields may be positioned 
at any frequency without changing the shape of the 
potential or stream lines. The concept of regarding the 
zeros and poles of the open-loop transfer function as the 
singularities of a flow field is a useful one in this respect,’ 
in that a physical feeling may be developed for the 
proper positioning of zeros and poles of the shaping net- 
works to reduce the effects from undesirable singu- 
larities in the original function or to provide desired 
effects on the final field. 
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far the pressure and the momentum. These expressions are given (7) Application of the method to a typical multispar thick-skin recent! 
in the Russian literature by Gurevich, but I have not seen them sweptback wing is made, including a tab form solution suitable the lift 
published in this country. for desk calculator techniques. of an u 
The pressure expression is derived as follows: iad Compeions hee senha ag Seances reaper _* 
Agreement is good in most respects. General analytical results fortun: 
—dP/p = (Uy + u) du + vdv (1) appear consistent with those from the large number of less ap- integra 
According to the linear theory p = po[1 — (Uou/ao?)] sothat by | Proximate methods! applying to rather different structure. summ 
substitution of the expression for p and integration we obtain: Basically, the analysis method is one of consistent deformations the cor 
P,— P 4? a in which an elementary stress system is selected which satisfies ; The 
= Uu + i = Bey - (2) | statical requirements and which provides calculable elementary tions 4 
po 2 2 deformations, including warpage of sections out of their planes tributi 


These deformations are not necessarily compatible initially in 
all parts of the structure. Additional self-equilibrating, or re Ap = 
dundant, stresses are then applied which produce corrective de 


where only terms up to second order are retained, My = U/dy is 
the Mach Number of the free stream and Pp is the free-stream 
static pressure. 


The rate of change of momentum (J) is of the form formations so that the final deformations are consistent. Several 


simplifying assumptions are utilized in evaluating the redundant 








dI/dr = Sf pu(uy + u) dy + Sf puv dx (3) etecesieuuiees 
. cer: : : The ell 
making the substitution for p given above yields Sn nae ee 
REFERENCES Wx) /al 
(1 ‘po) (dI, dt) = Sil ou dy - (1 = Ie") Su? dy + S ui dx (4) 1 Williams, M. L., A Review of Certain Analysis Methods for Swe pt-Wing face it 
“ Structures, Journal of the Aeronautical Sciences, Vol. 19, No. 9, p. 615 elliptic 
where the terms up to second order are retained September, 1952. 
2 Hall, D. B., Method of Stress Analysis of a Swept Back Wing, The Glenn scales. 
L. Martin Company, Report No. SDR-37, October, 1947 The st 
- . 7 Bowman, J. C., Method of Root Area Analysis for a Swept Back Wing a(x) ar 
Northrop Aircraft, Inc., Report No. MSRD-6, February, 1952 scale f: 
4 Cole, D. M., Jr., Static Test of Wing, Northrop Aircraft, Inc., Report No 
Note on a Method of Root Area Analysis for GM-409, June, 1949 
Swept Wings *’ Simms, M S., and Benson, A S Evaluation of Wing Static Test Result 
Northrop Aircraft, Inc., Report No. GM325, March, 1950 
which 
John C. Bowman 4 leeth 
Research Engineer, Special Weapons Division, Northrop Aircraft, 
Inc., Hawthorne, Calif. Li . . , 
ift of Highly Swept Wings 
October 27, 1952 gnly P 9g } aVV 
Harold Mirels 
ia A RECENT SUMMARY of swept-wing analysis methods,' men- — Corne// University, Ithaca, N. Y. 
tion was made of an approximate method by Hall? with ap- December 17, 1952 
parently limited applicability. Mr. Hall’s method has been 
the basis of an analysis* conducted at Northrop in which the fol- ies TWO-DIMENSIONAL CROSS-FLOW CONCEPTS of Munk! and 
lowing extensions and refinements have been made in an effort to Jones? have been widely used to obtain the aerodynamic ae 
: rl ae ~ i . . . “q. (5 
increase its usefulness: coefficients of low aspect ratio pointed wings, bodies, and wing isi a 
Sin th 


(1) The analysis is extended from a four-spar system to include body combinations. Relatively few investigations have been A 
-™.: ° : ; - ison s 
any number of spars. This was accomplished by means of a made, however, when the cross flow contains a free vortex sheet 
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This arises, for example, in studies of highly swept wings (Fig. 1) 
a shed vortex sheet whose 


Thus the 


For x > , the cross flow contains 
strength depends on the spanwise variation of lift 
cross flow is no longer independent of upstream conditions, and 
special 

More 
recently, Robinson‘ has presented an elegant general solution for 
the lift of highly swept wings 
of an unknown function, A(x), whose value depends on the wing 


the solution is thereby considerably complicated. A 
family of swept wings was considered in reference 3 


The solution is expressed in terms 


plan form and is determined from an integral equation. Un 
fortunately, Robinson’s work contains an error that makes his 
integral equation incorrect. It may therefore be of interest to 
summarize the lift properties of a highly swept wing, indicating 
the corrected integral equation and some limiting solutions. 

The notation of reference 4 will be used. Only chordwise sta 


tions x > co are considered. The spanwise lifting pressure dis 


tribution is 


a(x) — y? 


i 2oV ae RE’ 
didi a(x) 7 \ atx) — bx) 


. 4 
F’ a(x) — ¥ E y y? — bx) l 
; = ( 
N a%X(x) — b%(x)] K’ a(x) YaXx ) — 
The elliptic integrals are of modulus k’ = VY 1 — k?, where k = 


h(x)/a(x). The complete elliptic integrals are indicated by bold- 
The sine of the amplitude of the incomplete 
elliptic integrals is indicated inside the parentheses. Eq. (1) 
(16) in reference 4 


face italic type 


corresponds to the integrated form of Eq 
The shape of the pressure distribution curve depends only on 
a(x) and b(x). The unknown function A(x) appears only as a 


The chordwise variation of net lift is 


dL y (: =) A) 
= 2rpVa(: 1 — 2) 
dx seid K’ a(x) 


The unknown 


scale factor. 


Which corresponds to Eq. (18) of reference 4. 


) is obtained from the following integral equation: 
“(PE\K *) 
F | = = > 4 
B a(t) ' K 
a(t) a(t) i. 
r(“ )-K|- E ( )-e } 
a(x) a(x) 


a(t) ax) — b>(t)| 
tf 
a(x) \ ax) — axXt)$ 
Eq. (3) contains the complete elliptic integrals K and E, and it 


is in this respect that it differs from Eq. (15) in reference 4. Rob 
(12 He 


function A(. 


— a" = 


aV¥ ‘ark x ) 


lt (3 


mson’s error arises from his treatment of his Eq 
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properly states that a function of x [say f(x)] can be added to Eq 
(12) but claims that this function may be taken to be zero. Ac- 
tually, this function is not arbitrary, and it is not equal to zero 
If the solution of the swept wing is considered as a superposition 
of a simple delta wing and a cancellation wing,’ then f(x) can be 
shown to equal the value of the cross-flow stream function on the 
w = 0 portion of the cancellation wing. 

Let the leading and trailing edges be straight lines with slopes 
tan y; and tan ye, respectively. Define the slope ratio r 
tan y2/tan y:. Then, limiting solutions of Eq. (3) are: 


For r— « 


’ 


A(x) K’ aV tan jy 
; ag (4a) 
a(x) V1 — [b%(x)/a%(x)] 
for r 0, 
A(x)K’/a(x) = aV tan y, (4b) 
. « 
forr=1x~o, 
A(x)K’'/a(x) = aV tan y; (4c) 
for b(x)/a(x) << 1, 
A(x) K’ : { b(x)/4a(x) ge 
= aV tan ¥1 1 + 
a(x) r Jo In s 
1 / d(x) \? 
(4d) 
2\ a(x) 
Eq. (4a) is the same as Robinson’s limiting solution to his in- 


tegral equation. The correction to the integral equation does 


not effect the r > © case. Eq. (4c) corresponds to a two-dimen 


sional swept wing. Eq. (4d) was obtained by an approximation 
in which it is difficult to estimate the error involved, but it is an 
excellent agreement with numerical solutions of the integral 
equation. 

When the leading and trailing edges are parallel (ry = 1), the 


parameter 


A(x) K’ 


a(x)aV tan y; 


equals one for x = ¢, then decreases (with increasing x) to a mini 


mum value of 0.94 at b(x)/a(x) = 0.5, and then increases gradu 


ally toa value of 1. Thus, use of the constant value 1 seems jus 
tified for this case 

The writer has developed a general procedure for calculating 
the aerodynamics of highly swept wings. Expressions similar 
to Eqs. (1) to (4) have been obtained for the rolling and pitching 
case. Early publication is intended. 
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Effects of Angle of Attack on Local Skin- 
Friction Coefficients at M, = 
John Joseph Martin 
Research Engineer, North American Aviation, Inc., Downey, Calif 


November 11, 1952 


-PSarenagenae SKIN-FRICTION COEFFICIENTS are based on prop~ 
erties of air immediately outside the boundary layer, and 
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Number for parameters of Fig. 1. 


this leads to no particular difficulty in supersonic flow, since, 
when there is a change in the properties of air through a shock, 
results of this change may be calculated or obtained from tabu- 
lations for this purpose. It is the intent of this note, however, to 
show how local skin friction, and therefore in some measure 
heat transfer, may vary for a particular case with angle of 
attack. 


Specifically, the case is that of a flat plate with* Jf) = 4 
Ro = | x 107 
Tw/ Ta at i he /of : eed l )M,? 


(i.e., 7, = aconstant = insulated plate temperature) with shock 
wave attached. The local coefficients were computed on the 


basis of reference 1 in which 
Cja = Chall Ma, | Tx / ra) 


Then, given values of J) and Ro, these need be converted to Mg, 
Ra, and T/T, for the plate at an angle of attack; further, from 
these, Cyq may be found. 


* Subscripts denote angle of attack; Mo is free-stream Mach Number for 
plate at a = 0; Mais free-stream Mach Number for plate at angle of attack 


a 
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Now, 


Ma da Pa B ) 


| = Miata a a = R, 
— F a 


Po Ka 
and if 
Ma/wo = [1.505(Ta/To)'/?]/[1 + 0.505/(Te/T 
as is so for air when 7) = 393 °R. (see reference 1), then 
Tw/Ta = 1 + (1/2) (vy — 1)M,? 
as indicated. Also, 
Ma Pa 1 + 0.505/(T2/T) 
Ra = Matava Ha = R =e Paes 
Mo p 1.505 


and M/,asa function of a may be found in tabulations? 
lated. Then, since Ca is a function of Me, Re, and Ty/T,, com 
putationst may be made, and the results are plotted in Fig. | gs 


or caley 


functions of angle of attack. In this figure, for a > 0, the upper 
branch of ¢ya/¢yo and the lower branch of Ra/R> are for a strong 
shock and the remainder of the curves are for a weak shock. Fig 
2 is a graph of the calculations on the usual plot of cya vs. Rg, with 
Ma as a parameter, and shows the relationships among Ca, Ry, 


and M, that obtain in Fig. 1. 
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+ These were very kindly done by Mrs. R. B. Davis 


Comments on ‘‘Remarks on Dynamic Loads in 
Landing’”’ 


Frederic M. Hoblit 
Lockheed Aircraft Corporation, Burbank, Calif. 
December 10, 1952 


gy A PAPER, “‘Remarks on Dynamic Loads in Landing,” in the 
April, 1952, JouRNAL, Zahorski gives an example of calculated 
dynamic loads. In this example, dynamic magnification factors 
are calculated for four cases, with maximum values at the wing 


tip as follows: 


Maximum landing weight, limit drop 1.4 
Maximum landing weight, ultimate drop 2.2 
Normal landing weight, limit drop 3.3 

6.3 


Normal landing weight, ultimate drop. 


This variation is rather startling, and several conclusions might 
be drawn therefrom. 

Zahorski’s conclusion is that, in a dynamic landing analysis, 
it is important to use the correct forcing function and, in par 
ticular, to base the dynamic analysis on the Jimit loading rather 
than the ultimate loading. The desired factor of safety is then 
applied to the resulting accelerations after the calculations are 
completed. The writer heartily concurs in this conclusion. 

A second conclusion that some readers may have drawn is that, 
since the four forcing functions appeared to be of roughly the 
same shape, the results of a dynamic analysis are too sensitive 
to small changes in parameters to be of use in establishing de 
sign load factors. 

This conclusion is entirely too pessimistic 


A major part o 


the variation in dynamic factor in the present case is not due to 


small differences in the forcing functions. Instead, what appeat 
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to be small differences in the forcing function are actually in some 
cases substantial differences, and in other cases the mode shapes 
change considerably with changes in airplane gross weight. Fur 
thermore, whatever degree of sensitivity does actually exist can 
be detected and estimated by examining the mode shapes and the 
response time histories in the various modes 

The change in mode shape between the normal and maximum 
landing weight cases is due to the fact that not only the total 
mass changes but also the distribution of mass. Unfortunately, 
Zahorski did not include either the mass distributions or the mode 
shapes, but one can deduce from the data given that the differ 
ences in mode shape between the two cases are rather drastic 
For example, the signs of the accelerations for the first mode indi 
cate a shift of the node from inboard of the gear for normal land- 
ing weight to outboard of the gear for maximum landing weight, 
possibly as a result of a large addition of fuel in wing-tip tanks 
Naturally, a significant difference in dynamic factor should be 
expected, even if the forcing functions were identical 

The shape of the forcing function for any given one-time load- 
ing is characterized by (1) the time required to reach maximum 
load or “loading period,’’ (2) the period of time the maximum 
load is sustained, and (3) the time required for removal of load 
Each of these can have a great effect on the dynamic response 
In comparing the limit-drop and ultimate-drop cases, it is seen 
that item (2) above is completely different in the two cases, since 
the forcing function in one case is a triangle and in the other is 
a long trapezoid. A considerable difference in dynamic factor 
is therefore to be expected. 

Some of the variation between the limit and ultimate cases 
cannot be accounted for by the above considerations and prob 
ibly must be attributed to’rather small changes in forcing func 
tion (although some question might be raised whether a 40 to 50 
per cent change in loading period is a ‘‘small change’’). In any 
event, this variation could have been anticipated by examining 
the time histories of the responses in the various modes 

First, it is noted that, in this example, the higher frequency 
the second and third—contribute a large part of the 
The time at which the maximum total ac- 


modes 
total acceleration 
celeration occurs, however, is governed largely by the rigid-air- 
plane and first-mode accelerations. If at this time the second 
ind third modes should be in phase, they would add a large 
amount to the total; whereas if they are exactly out of phase, 
their contribution would be negligible. Comparing Zahorski’s 
limit-drop and ultimate-drop cases at maximum weight (his Figs 
jand 6), it is seen that the second and third modes both add at 
the critical time in the ultimate-drop case but are of opposite sign 
in the limit-drop case. If in the limit-drop case both of these 
modes had added, the dynamic magnification factor would have 
increased from 1.4 to 1.8, which is much closer to the ultimate 
lrop value of 2.2 

In this instance the difference in phasing is due to a difference 
in loading period. For the phasing to change as much as in the 
example would appear to be possible only when the loading period 
hanges by an amount that is about the same as, or greater than, 
half the natural period of one of the two modes. In other words 
the modes involved would have to have a relatively high natural 
lrequency. 


Incidentally, significant changes in phasing can also result from 
small changes in natural frequency, especially when the maxi 
mum total acceleration occurs a considerable time after the initial 
disturbance. If, for example, the critical time occurs after two 
cycles of the second mode and three cycles of the third mode, it is 
clear that only a 16 per cent change in the frequency of the higher 
mode will completely reverse the phase relation at the critical 
ume, so that if these two modes originally added, they would 
Considerable uncertainty as to natural frequency 
natural frequencies 


then subtract 
iS to be expected, since 
depend upon approximations to the stiffness properties of the 


calculated 
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structure and upon a mass distribution that may depart con- 
siderably in service from that chosen as a basis for the analysis 

Consequently, consideration should always be given to the 
effect of possible phase variation between the various modes at 
the critical time. The effect of phase variation on total acceler 
ation is, of course, limited to the sum of the acceleration ampli 
tudes of the two modes for which the phase relation is in doubt 

The shape of forcing function also affects the maximum accel 
eration independently of the relative phasing between modes 
Zahorski notes from his Eq. (14) that, for a trapezoidal forcing 
function, if the natural period of a mode is contained exactly 1, 
2, 3, etc., times in the loading period, the ensuing amplitude of ac- 
celeration in the mode is zero. If this number of times is in the 
neighborhood of 2, as for Zahorski’s third mode, it is seen that 
only about a 20 per cent change in natural frequency or loading 
period would make the difference between zero amplitude and a 
peak amplitude. So again the maximum response is seen to be 
sensitive to both the shape of the forcing function and the natural 
period of the mode. This sensitivity can be anticipated by look 
ing at the time history for a single case, and the possible range of 
response can readily be determined. 

It should be emphasized that, when the dynamic analysis indi- 
cates sensitivity of accelerations to shape of forcing function or 
to natural frequency, this sensitivity is a real phenomenon that 
would show up as scatter of test data if tests were made and will 
occur in service, especially as the airplane mass distribution 
varies. 

In summary, a dynamic analysis can be of real value. Ac- 
celerations obtained from a calculated time history must, how 
ever, be interpreted with judgment rather than accepted blindly; 
due consideration should—and readily can—be given to the pos 
sible range of results which might occur with reasonable variation 


of conditions 


A Reply—Single Degree Flutter of an Aileron 


H. L. Runyan, H. J. Cunningham, and C. E. Watkins 
Langley Aeronautical Laboratory, N.A.C.A., Langley Field, Va 
October 3, 1952 


—_ SUBJECT OF single degree of freedom flutter of a control 
surface has been discussed in reference 1, and some of the ma 
terial was presented in reference 2. In these references it was 
analytically demonstrated that a control surface having no aero 
dynamic balance (hinged at its leading edge) could flutter in a 
single degree of freedom. At about the time these papers were 
published, another treatment of this problem was published, 
with conclusions that appeared to contradict those of references 
1 and 2—that is, the conclusions of reference 3 stated that only 
ailerons with large amounts of aerodynamic balance could flutter 
Accordingly, an attempt was made to point out the basic reason 
for the apparent discrepancy between the two conclusions in a 
Readers’ Forum note. This note was then followed by another 
Readers’ Forum item® by the authors of reference 3 in which they 
strongly reiterated their previously made deductions. 

Inasmuch as aileron flutter in a single degree of freedom is theo 
retically possible without aerodynamic balance and, under cer 
tain conditions, may be dangerously significant, we are induced 
to present this brief reply to reference 5 with hopes that the mat 
ter will be laid to rest without false impressions of the theoretical 
possibilities of this type of flutter. This can be simply achieved 
with the equations and notation of reference 3 itself 

In the absence of structural damping, the conditions of aileron 
stability or instability are directly determined by the algebraic 
sign of 7,, the coefficient of the imaginary part of the complex 
Positive values of 77, indicate 


moment term 7, of re‘erence 3 
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Fig. 1. Values of-2 and v/bw for which 7; vanishes. (e= 0.4) 


dynamic instability, negative values indicate stable conditions, 
and the vanishing of 7; indicates the existence of a borderline 
state of equilibrium separating damped and undamped aileron 
motions. We are thus behooved to consider the nature of 7; as 
a function of the various parameters considered in reference 3 

namely, aileron leading-edge position e, hinge position c, aero- 
dynamic balance / = ¢ — e, and reduced speed parameter v/bw. 
If for a given value of e in the range 0 < e < 1 corresponding 
values of / and v/bw for which 7 vanishes are plotted, it will be 
found that such plots invariably consist of two distinct loops or 
branches, one loop corresponding to small values of v/bw and the 
Fig. 1 of reference 3 consists of par- 
that is, for 


other to large values of v/bw. 
tial plots of the first loop of 7; = O just described 
small values of v/bw—for several values of e in the range 0 < 
e < 0.7. A more complete plot of 7; = 0 for one representative 
value of e (e = 0.4) is shown here in Fig. 1. The only portion of 
the curves in this plot which is presented in reference 3 is indi 
cated by the solid curve, and results of more extended calcula 
tions are indicated by dashed curves. Values of 7; corresponding 
to points inside the loops are positive, indicating instability, and 
values corresponding to points outside the loops are negative, indi- 
cating stability. It will be noted that for ] = 0 (no aerodynamic 
balance) that 7; vanishes for v/bw = 33 (approximately) and for 
some value of v/bw > 33 for positive values of / in the range 0 < 
1 < 0.204. This corroborates the statement here and the conclu- 
sions of references 1 and 2 that aileron flutter is possible without 
aerodynamic balance. (Negative values of ] correspond to hinge 
positions ahead of the aileron leading edge and are excluded from 
the present discussion. ) 

As a final remark, it is pointed out that the mass moment of 
inertia parameter J/g is an important parameter that strongly 
affects the degree of stability or instability of an aileron. For any 
combination of values of v/bw and ] which causes the value of 7 
to vanish, there exists real positive value of 7g for which real 
values of flutter frequency are given by the equation of reals in 
Eq. (5) of reference 3. Furthermore for values of / corresponding 
to hinge position ahead of the aerodynamic center of pressure on 
the aileron (about 34 per cent aileron chord), the ratio 
of flutter frequency to aileron frequency will invariably be 
greater than unity, a conclusion apparently at variance with refer- 


ence 3. 
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October 30, 1952 


— COMMENTS—it is hoped that they constitute a closure 
are written without the collaboration of my coauthor in the 
paper ‘‘Aileron Flutter in a Single Degree of Freedom’’! because 
the distance between us makes such a cooperative venture diff 
cult. In spite of the omission of K. P. Abichandani’s name from 
this contribution to the Readers’ Forum, its content is based 
largely on his thesis? for the M.S. degree in Aeronautical Engi 
neering, completed and entered into the records of the University 
of Washington in December, 1950 

It is anomalous that a paper, having had but little circulation 
should become the subject of discussions** that are now entering 
into the third issue of the JoURNAL OF THE AERONAUTICAL 
SCIENCES where the paper under discussion was not published 
However, Runyan, Cunningham, and Watkins did not have the 
opportunity to discuss this paper in its originating publication 
because the Bulletin of the University of Washington Engineering 
Experiment Station is not a periodical. As a result, readers have 
been presented with the discussion of a paper that most of them 
have had no opportunity to examine; they probably enjoy it no 
more than hearing the punchline of an unknown joke 

The above Readers’ Forum contribution by Runyan, Cunning 
ham, and Watkins is substantially in agreement with some recent 
findings by Abichandani, reported to me in a personal communi 
cation,® and it largely resolves the apparent discrepancies that 
have formed the subject of the earlier discussions. 

It appears that the earlier work by Runyan® and by Runyan, 
Cunningham, and Watkins,’ concerned with single degree of free 
dom flutter of ailerons that are not aerodynamically balanced, 
was necessarily restricted to deal with that loop of their Fig. | 
We shall call it the “outer 
1 as the solution 


above which lies in v/bw > 33. 
loop.’”’ If one considers the curves of their Fig. 


to the equation 
1 = f(v/bw) (1) 


then one may regard the problem of references 6 and 7 as that of 
finding the zeros of f(v/bw). But, f(v/bw) has apparently only one 
zero in the outer loop; therefore, the authors were actually con 
cerned with locating that value of v/bw at which the outer loop 
crosses the line 7 = 0, 

We posed the more general problem of solving Eq. (1) 
ever, by using the tables of flutter coefficients tabulated by Smilg 
and Wasserman,’ we were necessarily restricted to deal with the 


how 


“inner loop” of Fig. 1 above because we considered the solution 
of this equation only in the interval 0 < v/bw < 16.67, the range 
of v/bw in reference 8. 

Although the upper branch of the inner loop was not plotted in 
our original paper,! it was not overlooked. In this paper, We 


wrote out the above Eq. (1). It is Eq. (8) of reference 1, 


l b; + be (db; + be)? b; (9 
i 2h, 4b,2 bs 


and the b; are defined by 
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Pp =a, + 1h 
T, = a2 + tbe 
Tg = a3 + 1b; 
P, =a, + 1h 


all functions of v/bw 


It was remarked! that / is a double-valued function of v/bw, and, 
in the earlier thesis? by Abichandani, both values of 1 have been 
computed for a variety of values of ¢e; they are listed in reference 
9 Tables 2.01 to 2.08. However, the upper branch of the inner 
loop was not plotted for the following reason: From the geometric 
relations of the illustration accompanying Fig. 1 of the above 
reply, it will be seen that the aileron chord is 1 — e. Thus, 
1 — e imply that the overhanging portion of the 


values of / 
aileron is larger than the entire aileron chord or that the aileron 
is hinged aft of its trailing edge. We may conclude that not 
only are values of 1 < 0 excluded, as observed in the above reply, 
but so are the values of / > 1 — e or that / must liein0 </<1— 

This excludes the major portion of the upper branch of the 
inner loop of the above illustration, because, in the example 
chosen, 1 — e = 0.6. In fairness to Runyan, Cunningham, and 
Watkins, it should be pointed out that we stated! that the posi 
tive sign in Eq. (2) always implies that the aileron is hinged aft 
of its trailing edge. Actually, this statement is true only for cer- 
tain values of e; but in those cases where it is not strictly true, it 
holds nevertheless for most values of v/bw. For instance, in the 
example of the Readers’ Forum item above, the lower branch only 
of the inner loop has significance in the interval 0 < v/bw < 4.8 
(approximately), while both branches are of interest in 4.8 
v/bw < 5.5. 

It remains to consider the parting remark of the above reply 
This is a repetition of an earlier statement by the same authors.* 
We find that we can no more agree with it now than we could 
In our original paper! there is an equation that is re 
It is Eq. (18) of reference 1 


earlier 
peated here for convenience 


(v,/v4)? = [1 — (w/we)?|/n (3) 


The same paper also contains the statement that the aileron 
divergence speed v, is imaginary if the aileron is hinged foreward 
of the aerodynamic center of pressure of the aileron. It will be 
observed that for this hinge location (v,;/v,)? 0 if the flutter 
speed v, is real, which is the case discussed by Runyan, Cunning- 
ham, and Watkins. But if (v,/v,)? is negative, then Eq. (3) above 
demands that (w/wg)? be larger than unity (because 7 is a real non- 
negative quantity). Thus, not only is the last statement of the 
above reply in accord, but any other statement would not be in 
accord, with our results 


In spite of some minor discrepancies that appear to exist still, 
we have reasons to note with pleasure the results of this discussion 
in the Readers’ Forum. In investigations that obviously were 
independently conducted, the problem of single degree of free- 
dom flutter of ailerons has been investigated. Prior, it was not 
even known that ailerons are capable of this type of flutter; asa 
result, the phenomenon was only imperfectly understood. In 
particular, only one régime of this flutter (the inner loop) was dis 
cussed in reference 1, and only one point on the outer loop was 
discussed in references 6 and 7. Now, both observations are 
drawn into a single fold, and, as a consequence, the phenomenon 
is better understood. This pleasant result appears to have made 


these discussions well worth while. 
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Approximate Calculation of Axisymmetric 
Flow Near the Vertex of a Body Traveling at 
High Supersonic Speeds 

A. J. Eggers, Jr., and Samuel Kraus E 

Ames Aeronautical Laboratory, N.A.C.A., Moffett Field, Calif 


November 3, 1952 


INTRODUCTION 


S TEADY, INVISCID FLOW in the region of the vertex of a curved 
noninclined body of revolution with attached shock wave 
has been studied by Shen and Lin! and Cabannes? by perturbing 
the conical flow solution of Taylor and Maccoll.’ It was found 
that fluid properties in this region may be obtained by performing 
rather lengthy numerical integrations of the perturbed conical 
flow equations. It has recently been indicated,** however, 
that flow about such a body traveling at high supersonic speeds, 
corresponding to values of the hypersonic similarity parameter 
in excess of 1, may be calculated with relative simplicity using a 
generalized shock-expansion method. This note describes the 
application of the method to the determination of flow in the re 


gion of the vertex 
SHOCK-WAVE CURVATURE AT THE VERTEX 


A curved axisymmetric body creates pressure disturbances 
that alter the inclination of the vertex shock wave. Fig. 1 illus 
trates the geometry of this phenomenon. The pressure disturb- 
ances (expansion waves for a convex body) travel along the first 
family Mach lines, F;, and interact with the shock wave. The 
reflected disturbances have been shown®’? to be so weak as to 
create relatively small pressure changes along these lines, pro- 
vided that neither the Mach Number nor the ratio of specific 
heats of the disturbed fluid is too close to 1. Now, at the high 
supersonic speeds under consideration, the variation of Mach 
Number from the surface to the shock wave is extremely small in 














y 
—s 
H 
Mo, % a M,P C4 
3 6 ie 
C F I 
AAW 
HLH /// 
L MM Ls HL MLL HLL x 
N 
Fic. 1. Geometry of the vertex region (a meridian plane) 








216 JOURNAL OF THE AERONAUTICAL SCIENCES—MARCH, 1953 





20 = 
Equation 4 —— 


Shen and Lin (Ret. 1) 


} 


(Ks 
Ky 


Curvature ratio, lf 
=~ 
is) 


Semi-vertex angle, Su, | 30° 




















-_ 4 


o ee ie 10 15 @ 
Free-stream Mach number, M, 





Fic. 2. Variation of curvature ratio with free-stream Mach 
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the unperturbed conical flow; therefore, the variation of pressure 
is also small. Thus, in the perturbed conical flow field created 
by the curved body, it is suggested that the variation of pressure 
along first family Mach lines is small. With this point in mind, it 
is clear, referring to Fig. 1, that the equation 


C Op c DOop .. *B Op 7 D Op ' 
5 = 
f oF, dF, + c 2s d I, aw + B OF, dF, (1) 


relating pressures on-the surface w to those just downstream of 
shock wave s reduces to 


Dod Bo 
f cas = ff P dw (2) 
C OS A OW 


Co eDoO 
f P aF, - J P dF, 
A OF, B OF, 


is indicated to be small indeed. But O0p/OW, the gradient of 
pressure with distance along the body, is, to the accuracy of this 
analysis,‘ the Prandtl-Meyer value. Noting that 00/d0s = K,, 
the shock-wave curvature, and 06,/O0W = Ky, the body curva- 
ture, we have, then, from Eq. (2), in the limit as points 4, B, C, 
and D approach the vertex, the following expression for the ratio 
of shock-wave to body curvature: 


- B , 
(<=) _ (Op /2%w) tien Si? aii 
’ 


| = n (0p/00)n A, B,C, D_, N Scrds 


since the quantity 


c 
The pressure gradient along the shock wave (0p/00@)y can be 
evaluated using the oblique-shock-wave equations, while the pres 
sure gradient along the body (0p/06,)y is given by the Prandtl- 
Meyer equation. With the knowledge that the local Mach 


Number is essentially constant and the difference between floy 
inclination angle and ray angle is small compared to ] in the 
conical flow at the vertex,‘ the ratio of integrals in this relatioy 
is readily evaluated. Thus, according to the generalized shock. 
expansion method, the expression for the curvature ratio may be 
written 


- ) 27M)? sin? vy — (y — 1) 2 tan (0v—5, ; 
: = ; : = ‘ 
KeJ Nn My)? sin 2@y sin 28 y tan By 
} 


where Sy, the Mach angle just downstream of the shock, and 
#y may be determined from existing conical flow calculations’ 4 
for given Wy and byy. 

Results of calculations of curvature ratio using Eq. (4) for 
several free-stream Mach Numbers and semivertex angles are 
compared in Fig. 2 with those presented by Shen and Lin. The 
agreement at higher Mach Numbers would seem to justify appli 
cation of the generalized shock-expansion method to the solution 
of this problem. 


FLow ConpbITIons ALONG A LINE A SHORT DISTANCE Downy 
STREAM OF THE VERTEX 


The variation of a fluid property, ¢, along the vertical line 
GH (see Fig. 1) can be expressed by a Taylor expansion as follows 


& = or + (0¢/Oy)R(y — yr) + (0°6/Oy*) ely — vR)?/2 + 


7) 


where R is a reference point on the line. Since the shape of the 
shock wave has been related to the shape of the body, we may 
proceed with this relation to calculate flow conditions along GH 
To illustrate, at the high supersonic speeds of interest here, the 
rate of change of 6 with y at point G for axisymmetric flow can 
easily be shown‘ to be 


(06/Ov)g = —sin 6¢ cos 6g/y¥Gq 6 
The second derivative (075/Oy")¢ is evaluated using the boundary 
conditions at Gand H. Thus the variation of 6 with y along the 
line GH is 
6 = 5g — (sin 6¢ cos 56G/y¥G) (9 — ya) + 
[654 — 5G + (sin 5g cos 6a/ya) (yu — ya) 


((y — yoe)/(yn — ya)]*? (i 


Similarly, knowing that the variation of static pressure p with 
F, at G is zero, the following relation for p obtains: 


2ypakw, cos (B + 5) 
p= pe - : (y — ¥e) + 


sin Bg sin 28¢ 


2y pckw, cos (B + d)e 
pu — pa + (yu — ¥G) |X 


sin Bg sin 2B¢ 


[(y — ye)/(yn — ye)l* (8) 


The variations of 6 and p/pu with (y — ye)/(yaw — ye) [Eas 
(7) and (8)], are compared in Fig. 3, in the limit of purely conical 
flow, with the variations obtained from M.I.T. tables.’ Good 
agreement is observed 

It is convenient to expand the total pressure ~; about the point 
H, since the gradient of p; with y is readily calculated at this point 
In an analogous manner, then, the variation of total-pressure 


with height is given by the expression 


Pt = Pin + [Opt 00)1Ks, cos 6#/sin (@ — 5)x] (y — yu) + 
[br — Pu — | (Op, /00) 1K sy cos 6x/sin (0 — 5)z)} (ya — YH)|X 
\(y — yx) /(ye — yu)? 
Knowing 4, p, and py, other pertinent fluid properties can be de- 
termined along the line GH. Thus, this line may be used as 4 
starting line for generalized shock-expansion solutions of flow 
about noninclined bodies of revolution traveling at high super 
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The method may also prove useful in starting char- 


sonic speeds 
teristics solutions for flow about a body. 
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Statement on References: On Air Forces on 
Oscillating Wings of Low Aspect Ratio 


|, E. Garrick 

Chief, Dynamic Loads Division, Langley Aeronautical Laboratory, 
N.A.C.A., Langley Field, Va 

November 21, 1952 


I THE INTERESTING PAPER by Lawrence and Gerber,! reference 
is made to an extension of some low aspect ratio wing theory 
of R. T. Jones to nonstationary flows attributed to me.? For 
the sake of the record, I should like to make it clear that the par 
ticular work referred to was done jointly by Dr. Eric Reissner 
and myself, as indicated in reference 3, during June, 1948, while 
Dr. Reissner was at the Langley Aeronautical Laboratory for a 
few weeks’ stay. It is pertinent to note that Miles has also con 
tributed to the same subject.‘ 
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Unsteady Aerodynamics of Potential Flows 


Aeronautical 


Compressible Flow Past a Cascade of Thin 
Airfoils 


G.V.R. Rao 
Veneral Electric Company, Lynn, Mass. 
November 5, 1959 


li WIDER USE of high-velocity flows through axial-flow com- 
pressors has shown the need for a theory of compressible flow 


Past a cascade at high velocities and small turning angles. The 
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to 
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well-known thin airfoil theory has been extended to an airfoil in 
By use of affine transformation, a linearized com- 
As shown 


cascade.» 2 
pressible flow can be obtained past a cascade of airfoils. 
in reference 2 for the incompressible flow, the method adopted 
here can be easily extended to the case of a staggered cascade 

The perturbation velocities caused by the airfoil profiles in 
cascade are assumed to be small compared with the velocity of the 
mean uniform flow. Each airfoil in the cascade can be replaced 
by a distribution of vortices situated on the chord lines of the pro 
files. Further, if the camber is small, these vortices can be lo 
cated on lines parallel to the mean flow direction, as shown in the 
Consider the row of vortices located at a distance Xo from 
the cascade axis. These all have the same strength y(xo) and have 
the coordinates given by 2, = Xo + int, where n takes all values 
from — ~ to + © including zero (Fig. 1) 

The following affine transformation is applied: 


figure 


x =€ 
By =n 


B =V1-— M 
where .V/ is the Mach Number of the mean flow 

It can be easily shown’ that, by choosing a singularity strength 
y(fo) = By(xo) in ¢-plane, the perturbation potential will satisfy 
Laplace’s equation. Thus, from the incompressible flow in the 
¢-plane, the velocity at a point (£, 0), induced by a row of vortices 
8t apart and located at distance & from the point (£, 0), is given by 


-z 


The velocity at the corresponding point in the linearized com 


¥(£o) 
wt, O) = 


2r(é — & — inBt) 
pressible flow is given by 
wx, O) = v(&, 0) 


~ Zed 


2r(—E — k — inst) 


¥(éo) 





Il 


In terms of the coordinates of the physical plane in which the 


compressible flow occurs, 
bey 


The above summation can be easily identified as 


By(Xo) 


2Q2r(x = Xo — inpt) 


— ¥( Xo) mx — Xo) 
wx, 0) = coth 
2t Bt 


Adding up the effect of all the rows of vortices given by —c/2 < 
2 <c/2 
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(x, 0) —1 ies ( tI a(x — Xo) j 
WX, = Xo) coth aXo 
a J. Phas Bt 


= —¢/2Z 


*o 
As t— o the above reduces to the perturbation velocity for an 
isolated airfoil, 
—1 


Visol. (x, 0) = ¥( Xo) fa dx» 


4 Xo c/2 wx — Xo) 


Hence, the effect of the rest of the airfoils in cascade can be given 
as 
c/2 
—1 (xX — Xo) 


( h~ 
— Visol. = y(X%o) | coth — 
i 2t ‘ Bt 


Veasc 


Bl 
dxo 
W(x — Xo) 
Spurr and Allen in their investigation‘ using a different approach 
give this effect as 


_ (Xo): (x — Xo) dxo 

6Bt? J _ 7/2 
after neglecting powers of [(x — xo)/nt] higher than second in the 
summation. These two expressions for the additional effect in 
the vertical velocity can be shown to agree if coth r(x — x9)+ 
8t were expanded and powers of (x — x»)/8t higher than second 
neglected. It can be readily seen that such an approximation 
overestimates the vertical induced velocity, thereby giving a 
higher incidence at the airfoil. Hence, their investigation showed 
a decrease in the lift-curve slope for an airfoil in cascade of high 
solidity at high Mach Numbers. Based on this they concluded: 
“The interference effects of the airfoils in cascade increase so 
rapidly with Mach Number that, contrary to the usual expecta- 
tion, a decrease in lift-curve slope occurs for high solidities.”’ 

Using the more complete expression for the additional induced 
velocity, the calculations by this author showed no such adverse 
effect. 
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The Pressure Distribution on Low Aspect Ratio 
Wings in Steady or Unsteady Incompressible 
Flow 


H.R. Lawrence 
Cornell Aeronautical Laboratory, Inc., Buffalo, N. Y. 
November 17, 1952 


— FOR COMPUTING the chordwise lift and rolling mo- 
ment distribution on low aspect ratio wings with straight 
trailing edges in steady and unsteady linearized incompressible 
flow have appeared recently in this JouRNAL.'» 2 In this note, 
the results of references 1 and 2 will be extended in two directions 
A physical interpretation will be supplied for certain mathemati- 
cal procedures employed in these papers, and a method for com- 
puting the distribution of lifting pressure will be derived. 

The notation of reference 2, which applies to the unsteady state, 
will be employed, but the steady-state results may be found by 
allowing the reduced frequency & to vanish. 


The lift per unit area may be expressed in terms of the Der. 
turbation velocity potential as follows: 


dL/dS = 2pU?[$x(x,y) + ikd(x,y)] 


A major result of references 1 and 2 is a method for computing 
the spanwise integrals g(x) and m(x) defined below 
P — B(x) —" 
g(x) = S P92) dy $(x,y) , 
(x) = J), dy yolxy) 
mx) = Jai. dy yo(x,y) 


In the steady case, the functions 2p U*g’(x) and 2pl ?m'(x) repr 
sent the lift and rolling moment per unit chord. Although 
knowledge of the above spanwise integrals of the potential js 
sufficient to compute a number of important aerodynamic prop 
erties of a wing, the distribution of lift pressures can only be cal 
culated when the velocity potential ¢(x,y) is known. 

In order to compute the velocity potential ¢(x,y) froma know! 
edge of the functions g and m, an additional physical assump 
tion must be introduced. Following reference 2, the integral 
equation for a lifting surface may be written in the form 


l B(x) dn 
WX, VY) — Wil(X,Y) = Waxy) = o(x,n) (4 
2 —B(x) y — 9 


where w;(x,y), the induced angle of attack, represents the addi 
tional terms in the integral equation for the lifting surface and 
Wa(x,y), the aerodynamic angle of attack, is defined by Eq. (4 
The potential ¢(x,y) may be expressed in terms of the aerody 
namic angle of attack wg(x,y) by inverting the integral Eq. (4 
as follows: 
—1 *B(x) dn Wal(x,n) Vv ‘BAx) — 9? 

(x,y) = ; f : 
TY ‘BAx) — re Bla I~” 
For a wing of very low aspect ratio, the induced angle of attack 
wi(x,y) becomes negligible compared with w(x,y), and the integral 
Eq. (4) becomes a form of the Jones* low aspect ratio theory for 
steady flow or the Garrick‘ extension of this theory for unsteady 
flow. For wings of higher aspect ratio, where the induced angle 
of attack w;(x,y) can no longer be neglected relative to the geo 
metric angle of attack w(x,y), a procedure similar to that em 
ployed in developing the Prandtl lifting-line theory may be used 
to obtain an approximate expression for the induced angle of 
attack w;(x,y). In the case where the spanwise lift distribution 
is symmetric, the induced angle of attack will be assumed to be 
constant over the span or w;(x,y) = wj's’(x) This assumption 
is seen to be directly analogous to the assumption of constant 
downwash angle over the chord in Prandtl lifting-line theory. In 
the case where the spanwise lift distribution is antisymmetric 
the induced angle of attack will be assumed to be a linear function 
of the spanwise coordinate or for wj(x,y) = ywi'?(x)/B(x). Ita 
clear that the general unsymmetric case is directly analogous to 
the Weissinger® assumption that the induced angle of attack is 
a linear function of the chordwise coordinate 

The induced angle of attack may be calculated in the follow- 
ing manner. In the symmetric case, Eq. (4) is multiplied by 
v/ 8% x) — y® and is integrated over the span. There results the 
formula 

F w(x) = [f(x) — g(x)] /(«/2)B%(x) 


where the function f is defined as follows: 


: B(x) i ‘ a 

f(x) = dy w(x,y) V B(x) — x? ‘ 
. - B(x) . 

‘ . P oe Witt 

It will be noted that, in the steady case, 2pl/?f'(x) 1s the lift 
per unit chord of a Jones wing. 


An expression for the induced angle of attack in the ant 


symmetric case may be obtained by multiplying Eq. (1 by 
yV 8%x) — y® and integrating over the span 
There results the formula 
(8 


w(x) = [e(x) — m(x)]/(KB?/16) . 








where th 


In the 
per unit ¢ 
It will 
comparet 


sistent ne 


Lawre! 
Subsonic - 
p 683, Oct 

Lawre 
4spect Ra 
\eronautit 

3 Jones 
Below and 

‘ Garric 
Third Ans 
B, 1951 

5 Prand 

t der 
18 

® Weiss! 
T.M. No 


Buserr 


D. Naylo 
Research 
Englar 


Novemb« 


ee 
of tl 
superson 
usually i 
easily ch 
in terms 
avoided 

The c 
values fe 
shock art 


Stead) 
sidered i 
indepen 
tively, t 


the per. 


Mputing 


t) repre 
nough 4 
ntial js 
ic prop 


be ¢ al 


| know! 
iSSuUMp 


integral 


e addi 
ce and 
rq. (4 
erody 
tq. (4 


attack 
itegral 
ry for 
teady 
angle 
e geo 
t em- 
» used 
sle of 
ution 
to be 
ption 
stant 
In 
Ptric, 
ction 
It is 
1s to 


ck is 


jlow- 


| by 


. the 


lift 


nti- 
by 














READERS 


where the function e is defined as follows 


l B(x) 
2 — B(x) 


In the steady case, 2pU%e’(x) represents the rolling moment 


dy yV B(x) — y*w(x,y) (9) 


ver unit chord of a Jones* wing. 
It will be noted that the sign in Eq. (9) has been reversed when 
ompared with references 1 and 2 in order to obtain a more con- 


stent notation 
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Busemann Coefficients for a Supersonic Airfoil 


D. Naylor 

Research Student, University of London, King’s College, London, 
England 

November 11, 19592 


INTRODUCTION 


a HAS BEEN CONSIDERABLE DEBATE on the correct values 
of the Busemann coefficients for the pressure on an airfoil in 
supersonic flew. This is entirely due to the lengthy analysis 
usually involved.1~* Computation is liable to errors and is not 
easily checked, but by resorting to an equation for the pressure 
in terms of the local Mach Number much of the algebra can be 
avoided 

The coefficients are first found when there is no shock; 


values for the combined case when a simple wave follows the 


their 


shock are easily deduced. 


SIMPLE WAVE 


Steady isentropic supersonic flow in two dimensions is con- 
sidered in which the velocity vector reduces to a function of one 
independent variable w, say. If p, p, g,'and a denote, respec- 
tively, the pressure, density, velocity, and acoustic speed of the 
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gas and if @ is the angle the velocity makes with some initial line, 
then the equations of continuity and motion along a stream line 
are 

= () 

0 


(pq)s + pqgOn 
Ps + PqIs 


where suffixes s, » denote differentiation along and normal to a 


stream line. 

Using dp = a*dp and eliminating p,/p yields the equation for 
the velocity in the form 

(a? — q*)gs + qa*O, = 0 
or, since g = g(w) and @ = @(w), 
(a? — @*)quwts + qa*Oytn = O 
Similarly, the equation of irrotational flow 
g8s = gn 
becomes 
Vous = Quen 


Therefore, eliminating wy 


an» 


(a? — q?)qw? + q?a°Oe? = O 
or 
dq/q d@ = +1/B 
where B= (7? — 1)'/?and M = q/a 
Now, g dq + (dp/p) = 0 and a®? = yp/p; therefore, the re- 
quired equation for the pressure is 
dp/p dd = y\M?/B (1) 


the positive sign being taken if the flow is one of expansion in the 
direction 6 decreasing. 
Logarithmic differentiation of ./ 


q/a gives 


d\M/M d0 = (dq/q a0) — (da/a dé) 
= [1 + (1/2) (y — 1)M?Idq/q do 
since 
q dq + |2ada/(y — 1)| = 0 
therefore, using dqg/q d@ = —1/B, 
dM/do = —[1 + (1/2) (vy — 1)M*]M/B (2) 


Differentiating Eq. (1) with respect to 6 and using Eq. (2) to sub- 
stitute for d.//dé@ yields 


d*p/p do? = yM*[(y + 1)M* — 4B?]/2B5 (3) 
Differentiating again gives, after short calculation, 
d°p/p de® = yM?[(y + 1)M8 + (27? — Ty — 5)M* + 
10(y + 1)M* — 12M? + 8]/2B" (4) 


In the case of an initial undisturbed stream parallel to the axis 
6 = 0, with Mach Number /, and pressure p;, Eqs. (1), (3), and 
(4) give at once the first three coefficients of the pressure series ex- 
pansion appropriate to a small expansive deflection @, actually 


negative for the clockwise expansion considered, 


p/p: = 1 + 6,0 + €26? + €:6% + 06%) (5) 
where 
G = yM,?/B,, B, = (M,? — 1)? 
Co = yM2[(y + 1)Mi4 — 4B,?] /4B;4 
Cs = yM,2 [(y + 1) MS + (27? — Ty — 5)ME + 
10/7 + 1)M\4 — 124M,2 + 8] /12B,7 
The following expansions are readily obtained from Eq. (5): 
po/p = & + (2c2 — 67)0 + (8¢3; — Beco + 3/0? + 0(8%) (6) 
poo/2p = C2 + (3¢3 — C¢2)0 + 0(8?) (7) 
poe /6p = cs + O18) (8) 
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By virtue of Eqs. (1), (3), and (4), these are just the expansions 
for the three Busemann coefficients, evaluated on the Mach line 
where the deflection is 0, as a power series whose coefficients are 
functions of the Mach Number ./;. 

If ¥ denotes the actual magnitude of the expansive deflection, 
then Eq. (5) can be written 

p/p: = 1 — ay + ey? — oy? + OY) (9) 
OBLIQUE SHOCK FOLLOWED BY SIMPLE WAVE 

Let MJ, and pp» denote the Mach Number and pressure in front 
of the shock and y denote the deflection through it (Fig. 1) 
From shock theory the pressure ~; immediately behind the shock 
is given by 

Pi/po = 1 + dio + dopo? + (b3 — D)Yo? + O(yo') (10) 
where },, bo, b; are the same functions of A/y as Cc, C2, C2 are of MW, 
and 

D = y¥7v¥ + 1)MA(5 — 8y) Mt — 418 — vy)? + 8] /96B 

Bo = (M.? — 1)? 

Because the formulas for the pressure, density, and velocity 
immediately behind the shock differ from the corresponding 
formulas if the change had been a small isentropic compression 
through the same deflection Yp only in the third order terms, it 
follows that in expressing each of the functions ¢, @, ¢; as power 
series in Wo, the formulas (6), (7), and (8) may be used if 4, C2, ¢s, 
and @are replaced by 4), be, b;, and Yo. Therefore, 

Cy = by + (2b. — d,?) Yo + (3b3 — B3byb2 + d,3) Wo? + Of Yo?) 
co = bo + (3b; — dybo)Wo + OC Yo?) 
= b; + O( Yo) 


Substituting these values in Eq. (9) and using Eq. (10) to elimi 
nate p, by simple multiplication gives the final expansion, 


| 


p/po = 1 + did + bop? + b3p3 — Do? + 0(¢3) 


where @ = yw — y, the resultant compressive deflection 

These results agree with those derived by Jones in reference 4, 
while the derivation of Eq. (10) is given by Laitone in reference 2. 
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A Note on the Solution for the Rectangular 
Membrane 


S. F. Borg* 
Stevens Institute of Technology, Hoboken, N. J 
October 27, 1952 


_— EXACT SOLUTION for the uniformly loaded rectangular 
membrane stresses and deflections has not been obtained 
Hencky! described a rather lengthy numerical finite difference 
solution but applied it only to a square membrane. The Hencky 
computations for the rectangular membrane have not as yet been 
carried through. A solution may be obtained based on a Fourier 
series analysis similar to the Levy Solution? for the nonlinear 


* This work was done while the writer was a Research Engineer for Gru 


man Aircraft Engineering Corporation 


large deflections problem of plates, but this also is a lengthy py 
merical procedure and has not been done. 


An approximate solution obtained by Féppl* assumes values fo; 


the membrane deformations which satisfy the boundary and syn 


metry conditions and evaluates the unknown constant coefficients 
by minimizing the total energy of the membrane. He applied 
this solution to the square membrane only and obtained the stress 


and deflection at the center, which values differ slightly froy 


Hencky’s. Neubert and Sommer‘ carried through the Fé) 


computations for the rectangular case and drew curves for 


stresses and deflections. In addition, these authors obtaing, 


reasonably satisfactory experimental verification for the Foppl 
and Hencky solutions for the square membrane, as did also Head 


and Sechler.6 The latter authors suggest a discrepancy between 
experimental and Féppl stresses for higher aspect ratio mem 
branes. Neubert and Sommer did not test the rectangular 
membrane 


Thus the status of the membrane problem appears to be a 
follows: (a) No complete exact solution has been obtained 
(b) the approximate square membrane solutions of Féppl and 
Hencky have been verified.experimentally; and (c) only the ap 
proximate F6ppl solution is available for the rectangular mem 


brane 


There is, however, a serious discrepancy in the Fépp! solution 
which casts doubt on its validity for a large range of rectangular 
membranes.! It is possible to obtain an exact solution for the 
semi-infinite membrane. This may be done, for example, by tak 
ing the limit of the semi-infinite tied-plate solution® as the plate 
stiffness approaches zero. Doing this, the value for the center 


stress in a direction parallel to the short side is 
=> 3 9.9 , 
or = 0.357 © park h? 


+ Since writing this note, the writer has been advised of a recent N.A.C.A 
Technical Note by Fliigge (No. 2612, ‘Stress Problems in Pressurize 
Cabins,”’ February, 1952), in which a brief statement is made concerning 


this weakness of the Féppl-type solution 
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in which 


p = transverse uniform load 
= short side of membrane 
E = modulus of elasticity 
h = membrane thickness 


If the Féppl solution is carried through to its limit, the semi 
infinite membrane stress oz is given by 
» 


or = 0.425 W p%a2E/h? (2) 


q difference of 19.1 per cent. (A similar comparison and dis- 
crepancy may be obtained for the center deflections. ) 

The reason for this discrepancy is not difficult to trace. The 
trigonometric function assumed by Féppl becomes less accurate 
as the aspect ratio increases. In the limit, b/a = «, the deflec- 
tion function should be a constant but instead maintains the 
sinusoidal form. 

Experiments indicate® that for rectangular membranes of as 
pect ration b/a> 3, the center portion of the membrane is, for all 
practical purposes, in the same stress and deflection condition as 
the semi-infinite membrane. 

The experimental and theoretical requirements outlined above 
are embodied in the curves of Fig. 1. The dotted lines are for 
o, and w (the center deflection) as determined by the Fdéppl 
method. The full lines drawn between the known limits b/a = 1 
Féppl) and a/b = O (exact) in such a way as to satisfy the known 
experimental requirements may be taken to represent actual con 
ditions closely enough for all engineering and design purposes. 

The values of the stresses and deflections in the rectangular 
membrane are of importance not only for problems involving pure 
membrane design but also in conjunction with approximate split 


rigidity designs of plates subjected to large deflections.* 
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A Measurement of Turbulent Boundary-Layer 
Profiles and Heat- Transfer Coefficient at M = 7 


Peter P. Wegener, Eva M. Winkler, and M. Sibulkin 
U.S. Naval Ordnance Laboratory, White Oak, Silver Spring, Md 
October 98, 1959 


geo geammnage IN A naturally turbulent boundary layer 
with heat transfer to the surface were made in the con 
tinuous 12- by 12-cm. N.O.L. hypersonic tunnel.! 
station was located on the centerline of one of the diverging walls 


The survey 
of a straight-walled two-dimensional nozzle. The free-stream 
Mach Number (M 
Stream velocity was 3,400 ft 
Number gradient inherent in this type of nozzle was approxi 
mately 0.05. per in. 
temperature (7)) were kept constant, 7) being 50°C. higher 


) at the survey station was 7.0, and the free 
per sec.; the free-stream Mach 


Wind-tunnel supply pressure (po) and 


than the temperature needed to avoid air condensation during 
the expansion. Wall surface temperature (7,,) was kept constant 
at Ty/T., = 5.7 by a water cooling system in the nozzle wall 


The static pressures measured on the wall and at the edge of 
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Measured turbulent boundary-layer profiles with heat 
transfer at a free-stream Mach Number of 7. 





Fic. 1. 


boundary layer agreed, and this value of static pressure was 
assumed to exist through the boundary layer. Pitot pressure 
was surveyed from wall to free stream, with a probe large enough 
The Mach Number de 
(The 


use of the perfect gas law and y = 1.4 implied in this procedure is 


to avoid errors due to slip flow effects.? 
termined from the Rayleigh formula is shown on Fig. la 


Total temperature ( 79’) distribu 
was surveyed total 
probe. This probe was calibrated in the free stream at about 
the same M and Re encountered in the boundary layer. The 
7,’ profile obtained after applying these calibration factors is 


valid for our test conditions. ) 


tion with a double-shielded température 


shown in Fig. lb. Velocity was calculated using the’ data 
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points and values of 7,’ obtained from a curve faired through the 
measured J)’ points. A straight line was drawn between the 
last To’ point and the measured value of Ty. This velocity pro- 
file is shown in Fig. lc, which combines the results of four inde- 
pendent runs at identical operation conditions. A '/;-power 
profile is also shown in Fig. lc for comparison. Since all the flow 
properties in the boundary layer are known, displacement thick- 
ness (6*) and momentum thickness (@) could be computed and 
are indicated in Fig. le. The Reynolds Number based upon 
free-stream conditions and displacement thickness is 438,900, 
while an equivalent flat plate Re would be of the order of 10’. 

Other measurements show that tunnel supply temperature (7)) 
is constant across the nozzle inlet. Therefore, the rate of heat 
transfer per unit nozzle width (Q) from the nozzle inlet to the 
boundary-layer station may be computed from 


Q = Cypatta(lw — Tod (1) 
o= = (} — 7) dy (2) 
% Pata Nie = Lof ~ 

Eq. (2) represents a definition of an ‘energy thickness” that per- 

mits the calculation of Q from the mass flow weighted defect 

of 7)’ with respect to Ty across the boundary layer at the survey 

station. On the other hand, QO may also be determined by meas- 

uring the discharge rate and temperature increase of the nozzle 

cooling water after a steady state of heat transfer has been estab- 

lished. It was found that these two values of Q were in agreement 

suggesting the reliability of the 7)’ measurement in the boundary 
layer. 

Measurements in the steel nozzle wall showed that the tem- 
perature dropped linearly with distance from the surface at the 
survey station after a, steady state of heat transfer had been es- 
tablished. Knowing the temperature gradient (d7/dy) and the 
thermal conductivity (&) in the steel wall, the rate of heat transfer 
per unit area, 

q = k(dT/dy) = (Te. — Tw) (3) 


could be computed, where h is the heat-transfer coefficient. _How- 
ever, to determine the Stanton Number (.St,,) (also called the 
dimensionless heat-transfer coefficient Cy, ), 
h k(dT/dy) 
St. = - — (4) 
CDP ot wo CrP otal Te — Tw) 

we need to know the insulated flat-plate temperature (7.), a 
quantity that we were not able to measure. Assuming the rela- 
tion between our recovery factor (7) and Prandtl Number (Pr) 
to be r = Pr'/* and taking Pr = 0.73, we can calculate 7, and 
in turn obtain St,, = 0.00075 from Eq. (4) and the measured 
data. Using Pr = 1 (and thus 7, = 7)), we find St,, = 0.00061 
However, after applying the relation® cy,,/2 = St,,Pr’/*, we ob 
tain cy. = 0.0012 regardless of which of these Pr we choose 
Knowing the local friction coefficient (cy, ), we calculate the shear 
stress at the wall (7,). While the heat-transfer measurement 
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Fic. 2. The u*, y* representation of the measured velocity 
profile at 17. = 7 and comparison with theory. 


may uot be too reliable, we are now able to tentatively plot our 
velocity profile in the u*, y* parameters‘ as shown on Fig. 2 
Since these parameters are of most value near the wall, the 
density and kinematic viscosity were based upon wall proper- 
ties. The portion of the velocity profile derived from measured 
M but only interpolated 7)’ (note actual distances from wall) is 
represented by a dashed line. For comparison, von K4rman’s 
incompressible flow semilogarithmic law‘ and assumed linear 
velocity distribution (w* = y*) in the laminar sublayer are shown 
on Fig. 2. Finally, a compressible flow profile calculated from 
Eq. (72) of Van Driest’s analysis® for our M,, and 7,,/T,, is shown 
(Aside from other assumptions inherent in this analysis, the pro- 
file has the boundary condition u = u,, at y = 6.) 

It appears from Figs. lc and 2 that the measured velocity pro 
file is similar to that found in low-speed flow except in the region 
extremely close to the wall, and this difference may be due to the 
lack of 7)’ measurements in that region. 
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Since its incorporation in October, 1932, as a scientific membership society to advance the art and science 
of aeronautics, the Institute of the Aeronautical Sciences, in over twenty years of service to aviation, has 


become the representative technical society for the aviation industry and profession. 


It is the only organiza- 


tion in the United States which brings together all of the sciences and branches of engineering and tech- 


nology which are applied to aeronautics. 


Publications 


Papers presented at meetings or submitted to the Edi- 
torial Board are published in the JourNAL oF THE AERO- 
NAUTICAL SciENcEsS Or the AERONAUTICAL ENGINEERING 
Review. The Journav prints in each issue full-length 
scientific papers on new research and developments 
in various fields applied to aviation. The Review keeps 
members and subscribers up to date on aeronautical news 
and literature through reviews of new books and periodical 
articles, government publications and trade literature and 
publishes papers on applied engineering and aircraft pro- 
duction. 

The AERONAUTICAL ENGINEERING CaTALoG, published 
yearly, is a guide to sources and specifications of materials, 
parts, and accessories used in the design and production of 
aircraft, aircraft engines, and parts. 


Résearch Facilities 


Through the Aeronautical Archives of the Institute, 
members have for their use the most complete aeronautical 
library and reference research facilities maintained by any 
technical society. At the New York building of the Insti- 
tute, The W. A. M. Burden Library of books, reports and 
periodicals, the Acronautical Index of subject files, bio- 
graphic files and bibliography, and the historical and art 
collections in the Archives are available for study. The 
Pacific Aeronautical Library at 7660 Beverly Boulevard, 
Los Angeles, maintains a reference collection and reading 
room for members and loans books to company libraries 
in the area. Through the Paul Kollsman Lending Li- 
brary, acronautical books are lent by mail without charge 
to members anywhere in the continental United States. 
In addition to a complete collection of standard aero- 
nautical reference works, the Kollsman Lending Library 
has available for loan new aviation books as they are re- 
viewed in the AERONAUTICAL ENGINEERING REvIEW. 


Membership 


In addition to membership in professional societies 
representing their particular fields, specialists applying 
their knowledge and experience to aeronautics find mem- 
bership in the Institute of great value in the contacts made 
with others engaged in aeronautical work and in the broad 
coverage of all aeronautical problems furnished through its 
meetings and publications. All applications.for member- 
ship are carefully reviewed by the Membership Committee, 
which recommends the grade of membership for which the 


INSTITUTE OF THE 


AERONAUTICAL 


applicant is found eligible. There are advanced grades of 
membership—Fellow and Honorary Fellow—to which 
members may be elected as their experience warrants. 
Thus, a member receives due recognition—more valuable 
because it is given by his colleagues—for his contributions 
to the progress of aeronautics. New applicants may be 
elected by the Membership Committee to the grade of 
Associate Fellow if they have been engaged for not less 
than ten years in the practice of the aeronautical sciences 
and shall have been in responsible charge of important 
scientific or engineering work, or shall have otherwise 
made outstanding contributions to the aeronautical sci- 
ences; to the grade of MEMBER if they have been en 
gaged in professional work for at least eight years and have 
acquired a recognized standing in engineering, design, or 
other special work applied to aeronautics; to the grade of 
Associate Member if they have acquired a recognized 
standing in an administrative capacity in the aviation in- 
dustry; to the grade of Technical Member if they are recent 
engineering school graduates or are engaged in technical 
aeronautical work. Student Members are admitted by 
application through Student Branches organized at their 
schools. Aeronautical companies, as well as individuals, 
are affiliated with, and participate in, the support of the 
Institute through Corporate Membership. 


Dues 


An entrance fee of $10 is required of all new members 
who apply for admission to a grade of membership to be 
specified by the Membership Committee (Corporate Mem- 
bers, Foreign Members, and Student Members being ex- 
cepted). 

Annual membership dues include a subscription to the 
AERONAUTICAL ENGINEERING Review. Dues for the vari- 
ous grades of membership are as follows: Fellows 
$16.50; Associate Fellows—$16.50; MEMBERS—$13.50; 
Associate Members—$13.50; Technical Members—(26 
years of age and over) $13.50, (under 26) $6.50. 

The regular subscription price for the JouRNAL OF THE 
AERONAUTICAL SciENcEs and the AERONAUTICAL ENGINEER- 
1NG Review is $12 and $3.00 per year, respectively. Mem- 
bers may subscribe to the JourNat at $4.00 per year 

Those who are engaged in aeronautical work or who 
have an interest in any technical phase of aviation will find 
the services of the Institute to be a necessary aid in keeping 
abreast of developments in these times of such rapid prog- 
ress. An application form and further information about 
membership can be obtained from an Institute member or 
by writing to the Secretary. 


SCIENCES, INC. 


2 EAST 64TH STREET 


NEW YORK 2], N.Y. 
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SUGGESTIONS for CONTRIBUTORS 


to the Publications of the 


INSTITUTE of che AERONAUTICAL SCIENCES 


The Institute of the Aeronautical Sciences invites both members and nonmembers from any 
country to submit papers for publication in the JOURNAL OF THE AERONAUTICAL SCIENCES and the 


AERONAUTICAL ENGINEERING REVIEW. 


not pay for contributions. 


The Institute, following the practice of other societies, does 


The following directions for the preparation of papers, if followed by authors, will save corre- 
spondence, avoid the return of papers for changes, minimize the work of preparation for the printer, 
and save the expense due to the charges made for “author's corrections.”’ 





Manuscripts: The original typewritten copy of the paper is 
desired, double or triple spaced on one side of white paper sheets, 
consecutively numbered. There should be wide margins to 
allow for the marking of directions to the printer. Correcting, 
changing, or adding to papers after they are in type is costly. It 
is, therefore, imperative that papers submitted be in final form. 
Typographical errors may be corrected on proofs, but if authors 
wish to add material, they may do so at their own expense. In 
mailing, drawings may be rolled but manuscripts should be sent 
flat. Send by first-class mail (register if you wish for your own 
protection) to the Editorial Office, Institute of the Aeronautical 
Sciences, 2 East 64th Street, New York 21, N.Y. All manu- 
scripts will be examined by the Editorial Committee and by the 
Editor. Authors will be advised as promptly as possible whether 
the paper is acceptable for publication. 


TitLes: The title of the paper should be brief. The name 
and initials of the author should be written as he prefers. The 
use of the full name of an author is advocated because of the fre- 
quent duplication of initials and surnames which sometimes 
makes it difficult to establish the identity of the author. This is 
particularly important for large annual indexing and abstracting 
services. All titles and degrees or honors are omitted. The 
name of the organization with which the author is associated 
should be placed after his name on a separate line. The date on 
which the paper is received will be inserted by the Editor. The 
author’s title should be indicated in a footnote. 


SUMMARIES OR ABsTRACTS: An abstract to be printed at the 
beginning should accompany each article. It should be adequate 
as an index and as a summary. It should contain a statement of 
major conclusions reached, since summaries in many cases con- 
stitute the only source of information used in compiling scieatific 
reference indexes. Abstracts printed in other journals, espe 
cially foreign, in most cases, consist of summaries from printed 
papers. The summary should explain as adequately as possible 
the major conclusions to a nonspecialist in the subject and should 
contain from 100 to 300 words, depending on the length of the 
paper. 


SUBHEADINGS: Subheadings should be inserted by the author 
at frequent intervals. The work of editorial preparation will be 
simplified by the author providing many subheadings 


MATTER USUALLY DELETED: Photographs or illustrations 
of little technical interest and not showing advances in general 
practice. Too detailed tabular matter (general results of such 
tables may be included in the text). Lengthy descriptions of 
materials or processes or of preliminary experiments or theories 
that preceded final results; salient features only are of interest. 


REFERENCES AND Footnotes: References should appear as 
footnotes only, numbered consecutively and grouped together at 
the end of the manuscript. The arrangement should be as fol 
lows: (For books)—! Durand, W. F. Aerodynamic Theory, 1st 
Ed., Vol. 1, p. 23; Julius Springer, Berlin, 1984. (For Maga- 
zines )—! England, C. R., Crawford, A. B., and Mumford, W. W., 
Some Results of a Study of Ultra-Short-Wave Transmission 
Phenomenon, Proc. I.R.E., Vol. 20, No. 12, pp. 481-482, March, 
1933. Please give author, title, edition, volume, page, publisher, 


and date of publicition as indicated. Omission of one required 
fact causes much extra editorial work and possible inaccuracies, 


ILLUSTRATIONS: Illustrations should accompany manuscripts, 
and each should always be referred to in the text, by number, 
Drawings or graphs should not be larger than 12 X 16 inches, 
and must be made with jet black India ink on-white paper or 
tracing cloth, the latter being preferred. Do not use typewriter 
for lettering. The smallest lettering on 8- & 10-inch figures should 
be no less than !/; inch high. Cross-section paper (white with 
black lines) may be used but it should not have more than 4 lines 
per inch. If finer ruled paper is used, the major division lines 
should be drawn in with black ink, omitting the finer divisions. 
In the case of finely ruled paper, only blue-lined paper can be 
accepted. Tracing paper and blueprints are not acceptable. 
Lettering and all markings must be large enough to be readable 
after reduction. Mail rolled or flat; never fold. Drawings 
that cannot be reproduced (including pencil drawings) will be 
returned to the author for redrawing thus delaying publication of 
the paper. Photographs should be distinct and show clear black 
and white contrasts. They must be on glossy white paper. 
Avoid round and oval photographs. 


CAPTIONS AND LEGENDS: Legends or captions must accompany 
each drawing or photograph submitted. If written on the draw- 
ing or photograph, they should be placed below and well outside 
the part to be reproduced. Each table should have a caption 
such as Table 1, Table 2, Table 3, etc. Captions should be com- 
plete in themselves so as to make the data intelligible to the 
reader without reference to the text. A duplicate list of captions 
for figures should be included as the last page of tie manuscript. 
Use ‘Fig. 1’ (not Figure 1), Figs. 3 and 4, etc., in both the text 
and the numbering of illustrations. In the text, ‘Eq. (1)” or 
“Eqs. (1) and (2)” is used, not ‘“‘Equation (1).’’ In captions 
and legends, except for ‘‘Fig.’’ and ‘‘Eq.,’’ and in table headings 
write all words in full; do not abbreviate. 


MATHEMATICAL WorRK: Only the simplest formulas should be 
typewritten; all others should be carefully written in pen and 
ink, the writing to be large enough so that ample room is provided 
to mark mathematical matter for the printer. A considerable 
spacing for marking should be allowed above and below all equa- 
tions. All complicated equations should be repeated on separate 
sheets with plenty of space left for marking. The solidus should 
be used for simple fractions appearing within the text. Make 
all expressions clear to the typesetter. Greek letters used in 
formulas should be clearly designated by name on the margin of 
the manuscript. All symbols should be clearly written and care- 
fully checked. The difference between capital and lower-case 
letters should be clearly distinguished and care taken to avoid 
confusion between zero (0) and the letter (0), between the 
numeral (one) and the letter (ell) and the prime (‘), between 
alpha and a, kappa and k, u and mu, v and nu, n and eta. All 
subscripts and exponents should be clearly marked, and dots and 
bars over letters or mathematical expressions should be avoided. 
Avoid complicated exponents and subscripts. When it is neces- 
sary to repeat a complicated expression, it should be represented 
by some convenient symbol. 

NOMENCLATURE AND ABBREVIATIONS: Standard abbreviations 


should be used, and it should be noted that most abbreviations 
are lower case such as m.p.h., b.m.e.p., i-hp., b.hp., hp., ete. 
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